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Two-Terminal Dielectric Measurements Up to 6 x 10° Hz 


Martin G. Broadhurst and Anthony J. Bur 
(February 23, 1965 


\ two-terminal dielectric specimen holder has been constructed and used to make 
dielectric constant and loss measurements on a single disk specimen at room temperature 
over a frequency range from 10 to 6108 Hz. The measurement procedures are outlined 


and a detailed analysis of the working equations and measurement errors is presented, 


ie Introduction of the measurement techniques used here are related 
to techniques described in a previous publication 
It is often important for the study of dielectric from this laboratory [3]. 
properties that measurements be made on a single 
specimen using a minimum amount of the material | 2. Design and Construction of Specimen 
over as wide a range of temperatures and frequencies HB ld. 
as re: Low-frequency measurements Can con- seed 
veniently be made on a single disk-shaped spec imen, , 3 a a ; 
but at higher frequencies conventional equipment A schematic diagram of the specimen holder is 
requires the use of evlindrieal or ree ti angul: ar speci- shown In figure 1. The cas », Cap, and elect rodes are 
mens of ari us sizes to fit diffe rent microw: ive made ol brass, gold-plated to reduce surtace losses 
cuides lt Cavities. If tl he dielectric properties of a at high lrequency. lhe lower ele ctrode, which IS 
material are sensitive to differences in specimen operate dat high pote ntial, is Suppor ted and insulated 
preparation, then the use of different specimens to from the case by a fused silica rmg, which in turn is 
cover the desired frequeney range could lead to supported by the bottom rim on the holder case. 
results which seriously misrepresent the frequency he lower electrode and silica ring are held fixed by 
dependence of the dielectric constant and loss even | 5@™ all beads of epoxy cement. The upper on ctrode 
though the individual measurements may themselves | 8 Connected securely to a brass pin, finely threaded 
be quite accurate. on the upper end to accommodate the at ing nut. 
The two-terminal parallel-plate dielectric speci- 
men holder described in this paper was designed to 
make measurements on single disk-shaped specimens 
(about 1 to 2 ¢ of material) over as wide a range of 
frequencies as possible. Minimum specimen and 
holder size together with simplicity of construction —_ 
and operation have been emphasized in order that Sie on 
the cell could easily be fitted with an insulating eae 
jacket and operated at temperatures from —200 to MICROME 
oon “<;. The variable tempe rature oper: ation of MOUNT 
the holder will be reported in a later paper. 
In this paper, we will discuss the construction of 
the holder currently in operation at room tempera- 
ture, and develop the equations for its use with the 
low-frequency Scheiber bridge [1],) the General é CASE 
Radio Models 716C Schering Bridge and 1615A 
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Transformer Bridge, and the Boonton Radio Corp- | compression S SHORTING 
oration Models 260A, 280A, and 190A Q-meters. SPRING STRIP 
Dieleetrie constant and loss data obtained with the 

holder for a disk of commercial poly(methyl metha- UPPER 
crvlate) over a frequency range of 10>? to 610° Hz ELECTRODE 
will be presented to illustrate the use of the holder. FUSED SILICA 
The experimental errors will be analyzed and dis- CONNECTORS — ae 

















cussed in detail. Use of the holder to measure liquid LOWER ELECTRODE 
specimens is described in a separate paper [2]. Some | th 
FIGURE e A sche matic diagram of the two-tern 


1 Figures in brackets indicate the literature references at the end of this paper. specimen holder 





The’compression spring holds the spacing nut against 
the shoulder on the holder cap and rotation of the 
spacing nut raises or lowers the electrode and hence 
regulates the spacing. The spacing ‘between the 
electrodes is determined by measuring the position 
of a steel ball, mounted at the top of the upper 
electrode assembly, with a micrometer barrel mounted 
rigidly to the holder cap. Electrical connection to 
the upper (ground potential) electrode is made by 
soldering four flexible gold foil strips between the 
upper electrode and case. Such an arrangement was 
preferred to conventional metal bellows mounting 
because of the former’s shorter electrical path. 
Electrical connection to the holder is made through 
two connector caps soldered to the case and lower 
electrode, which fit snugly over pins suitably attached 
to the measuring equipment. The holder electrodes 
are 2.54 cm in diameter. The micrometer measures 
10 » per division and readings are estimated to 1 y 
(1 w=10™ m). 


3. Theory of Measurements 


Each dielectric constant and loss determination 
involves a specimen in and specimen out measure- 
ment of the holder impedance at the measurement 
frequency, and two calibration measurements of the 
empty bolder capacitance at some convenient audio 
frequency. The calibration measurements can either 
be done separately and compiled in the form of a 
curve of measured capacitance versus electrode 
spacing which is then applicable to all measurement 
frequencies, or they can be performed at the time of 
each determination with increased accuracy. 

Figure 2 shows a schematic diagram of the two- 
terminal holder for the two measurement frequency 
conditions, together with the assumed equivalent 
circuit and the effective series resistance and capaci- 
tance of the cell. The symbols used in figure 2 are 
explained in table 1. The two audio-frequency 
calibration measurements can be represented in a 


Figure 2. A _ schematic representation of the 
holder with (a) and without (b) a specimen 
sponding assumed equivalent circuits. 


two-terminal 


and the corre- 


The symbols are defined in table 1. 


way similar to figure 2b except that the effects of 
lead inductance can be ignored. When the measure- 
ment frequency is itself in the audio range then the 
general measurement procedure outlined here can 
be simplified as is described in the following section. 


TaRLE 1. Description of symbols used in text 


A ; ; 
= €'€y ; equivalent parallel capacitance of the 


x 


specimen. 


A, : . 
€05 =equivalent parallel vacuum capacitance 


of the specimen. 


equivalent parallel resistance of the specimen. 


= €05 - 7, = Capacitance of the airgap between the 
specimen and electrode. 

€0 a =capacitance berween adjacent areas A, of 
the electrodes. C; corresponds to t. 

the capacitance reading of the standard capacitor 
in a bridge when the cell is connected. ( 
corresponds to t; and C,, to tr. 

capacitance between the remainder of the elec- 
trodes including that portion outside the area 
A, and the edges. C,, corresponds to spacing 
t; and C,, to tr. 

capacitance between the leads when the holder is 
in the measurement circuit. 

error capacitance due to distortion of the electric 
field (see discussion of errors). 
equivalent resistance and inductance 
leads when the holder is in the 
circuit. 

effective resistance 
holder with specimen in. 
capacitance between the leads when the holder is 
in the calibration circuit. 

the difference in frequency between the upper and 
lower half-power points of an L C-R series resonant 
circuit. 
effective 
holder with specimen out. 

thickness and area of specimen. 

separation of electrodes at settings 1 and 2. 
dielectric constant. 

permittiv ity of free space. 

dielectric loss. 


cl 


of holder 


measurement 


series and of 


capacitance 


series resistance and capacitance of 


9 
Ce 


g Re \Ce+Cc, 


FieureE 3. Two simple transformations used to derive the 


working equations. 
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In order to relate the equivalent circuits in figure 
2 to their effective measured circuits, the approxim: ite 
transformations shown in figure 3 were used. The 
approximations will be suitable for most low loss 
materials such as many polymers, glasses, and 
organic compounds. The validity criteria are shown 
in figure 3, and in general if ¢’’/e’<(0.1 the trans- 
formation approximations will introduce less than 
1 percent error in the results. 

At the measurement frequency, the effective 
measured series capacitance and resistance of the 
holder with specimen in, Cy; and Ry;, and _ 
specimen out, Cy, and Py, can be expressed a 
follows: 


Cn 


wo 


R, 


Rag 
tv (wl 'R.)' 


( Ho 


wD aril Ho 


ois Ears 


where 


and 


One can write an expression for C, which is inde- 
pendent of frequency or the type of electrical equip- 
ment used. If Cy, is experimentally set equal to 
Cyr:, eqs (1) and (3) can be equated with the result, 


-C.,)—C"] , 
C)—0"} (7) 


From the calibration measurements (circuits corre- 


spond to that of figure 2b ignoring the lead induct- 
ance), 


one finds 


which can be substituted into eq (7) to give 


CI[(C..—C, 


C,—C") 
C, ibe 


~~ C,-[(C,— C,) FG OY (0) 


Writing C, and C; in the numerator in terms of 


4 


measured quantities and applying the definition of 


| recorded. 


e’ one finds 


t 
fi 


The last term brackets is the airgap correction 
term whose value is close to unity. It should be 
noted that this airgap term will be used in this form 
throughout the text. The effect of the error capac- 
itance ©” is also small (see discussion of errors) and 
hence the value of ¢’ is primarily dependent on the 
difference between two empty holder capacitances 
measured at audio frequencies and on the sample 
dimensions and initial holder spacing. 

Since the holder is a part of the measurement 
circuit when the specimen is both in and out of the 
holder, the lead losses Ry; can be eliminated by 
combining eqs (2) and (4), giving, 


CoRR 


ARy Ru - Rigo 


ss G. w® Cy (Cr+ C,) - 
€ ( DD: CO ARx k C. : ] : 


Unfortunately, there is not a convenient way to 
make resistance substitutions for the specimen loss 
so that the calculation of R, will depend on the par- 
ticular methods used at the measurement frequency. 


(13) 


4. Measurement Procedure 


4.1. @-Meter Procedure of Measurement 


For Q-meter measurements in the frequency range 


10° Hz to 2.5 10° Hz, two Q-meters were used, the 
Boonton Radio 260A Q-meter (100 kHz to 20 MHz) 
and the Boonton Radio 190A Q-meter (20 MHz to 
250 MHz). The measurement procedure used is as 
follows: The holder is mounted on top of the Q-meter 
and connected in parallel with the Q-meter capacitor, 
C,, with the shortest possible leads. With the 
specimen in the holder, the circuit is brought to 
resonance by tuning the (J-meter capacitor, C,, to 
©). The values of C/, the elec sold spacing t,, and 
the Q of the circuit with the specimen in the holder, 
Q;, are recorded. The specimen is then removed 
from the holder and the electrode spacing is de- 
ree ased until resonance is achieved while holding 
, fixed at Ci. The electrode spi icing tf and the Q 
of the circuit with the specimen out, Q,, are 


The two resonant situations, with the specimen 
in and with the specimen out, occur at the same 
frequency and the total capacitance in the resonant 
circuits must be the same in both cases. The value 
of the total capacitance, Cr, is taken as the value of 
CO, when the Q-meter is tuned to resonance with the 
holder removed. 
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For the Q-meter measurement, the circuits of 
figure 2 may be viewed in parallel with the Q-meter 
capacitance C Considering figure 2a, one finds 
that if (,,; and R,,; are placed in parallel with C; 
then the transformation of figure 3b can be used 
to reduce the circuit to a simple L—C-R series circuit, 
for which the total series capacitance, the total series 
resistance and the Q of the resonant circuit with the 
specimen in the holder, C;, R,;, and Q,;, are given as 


(14) 


(15) 


and 
] @ wh ( yi (16) 
where (’y,, and P,,; are defined by eqs (1) and (2). 
Likewise, with the specimen out and with the 
holder at spacing fy, the resonant circuit can be 
represented by figure 2b in parallel with C;. The 
total capacitance, the total series resistance and the 
@ of the resonant circuit with the specimen out of 


the holder, Cr, F,,, and Q,, are given as 


ana 


19 


where Cy, and R,,, are civen by eqs (5 and (4). 
Equating the expression for Cy; m eqs (14) and (17 
will yield eq (11 
measured value of ¢’. 
An expression for “3h is obtained from the difference, 


of the previous section for the 


Q>?—Q,', using eqs (16) and (19 Thus, 


2()) 


dividing through and usine 


, one obtains, 


Rearranging, 
( H ( H 


) 
¢ 


In the above equation, the term in brackets is the 
dominant term and it alone would be present if there 
were no airgap or lead inductance. The term 
(Cyo/Co)* results from the effect of lead inductance, 
which becomes significant above 30 MHz (see eq 
The third term, [(C,+C,)/C,}*, is the airgap 
correction. 

Equation (21) can be reduced to following working 
equation: 


”) 
())- 


In the above equation Cy, bas been replaced by 
Cy—C, in xecordance with eq (17). The capaci- 
tances, Cy and C”, are obtained from the calibrated 
(-meter capacitor. The capacitance Cy, is obtained 
from the audio frequency calibration. Omitted 
from the above equations are the high-frequency 
corrections which must be applied to the values of 
@ and C, as read from the 190A Q-meter. The 
manner in which these corrections are carried out is 
explicitly stated in the 190A @-meter manual. 

A more detailed description of Q-meter dielectric 
measurements may be obtained from a paper by 
Hazen [4] in which he shows how the Q-meter circuit 
may be adapted for the measurement of very low 


losses. 


4.2. Two-Terminal Measurements at Audio 
Frequencies 


The General Radio (GR) 716C Schering bridge 
was used for audio frequencies from 10° Hz to 10° 
Hz, and the GR 1615A transformer ratio arm bridge 
was used from 510° Hz to 2*10* Hz. 

For audio frequencies the effects of the lead in- 
ductance can be ignored, and the circuits of figure 2 
without the inductance are considered here. With 
these bridges one measures the dissipation factor 
and the effective series capacitance of the unknown, 
However, when ¢’’/e’<0.1 the series capacitance 
and the parallel capacitance of the unknown are 
equal (within 1 percent and this is treated 
here. ‘/e’ >0.1 the GR bridge manuals should 


Case 
For e’’/e 
be consulted. 

If the holder and specimen are connected as a 
two-terminal device in an arm of the bridge and if 
the procedure of measurement which was used with 
the @-meter is also used with the bridge, then the 
11) and 


result for e’ is the same as eq 


where D; and D, are the dissipation factors with the 
specimen in and out of the holder, and Cr is the 
total capacitance in the arm of the bridge which 
contains the holder. with the (-meter 
measurement, Cy, must be the same for the D; and 
D, measurements. 

A variation on the above experimental procedure 
which utilizes the standard capacitor in the bridge 
is preferable when using a bridge circuit and null 
detector for a two-terminal measurement. First, 
a bridge balance is obtained with the specimen and 
holder in the bridge circuit. Then ¢,, D;, and the 
standard bridge capacitor reading (,; are recorded. 


Just as 
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The sample is taken out of the holder. Keeping 
the spacing fixed at ¢,, the capacitance of the 
standard capacitor, (,, is varied to rebalance the 
bridge. The specimen-out values, and D, are 
recorded. 

If «a bridge measurement is made by changing 
the standard capacitor while holding the electrode 
spacing at ¢,, then a distinction must be made 
between using the substitution method of measure- 
ment and using the bridge as a direct reading 
instrument. In the substitution method, the un- 
known is connected in parallel with the standard 
capacitor, whereas for the direct method the un- 
known is connected in the bridge arm electrically 
opposite the standard capacitor. The GR 161: 5A 
conveniently used to measure C and D 
directly. The substitution method is used with the 
GR 716C bridge by connecting a precision air 
ballast capacitor to the “unknown direct” terminals. 
When the electrode spacing is ¢, for both ‘in’ and 
out”? measurements and a direct measurement is 
made, the total capacitance in the unknown bridge 
arm will be different for the “in” and ‘out’ balances. 
Let Cy; and Cr, be the total capacitance in the 
unknown bridge arm with the sample in and out of 
the holder respectively. For a direct measurement 
we have, 


IS most 


“c 


(24) 


(25) 
from which €’ can be calculated in terms of measurable 
quantities. Using Cy;—Cy,=C.,.—C, we have for a 
direct measurement 


with the GR 


conveniently true 


1615A 
that 


For a direct measurement 
transformer bridge it is 
Crn=C,, and Cr:=C,:. 

On the other hand, for the substitution measure- 
ment, C,,—( Cr and CC C For ¢ and &’’ 
we have, 


where the bracketed term in eqs (26) to (29) is the 
airgap correction term. To obtain the total capac- 
itance, Cy, in the unknown arm of the GR 716C 
Schering bridge a separate bridge measurement is 
necessary. During the “in” and ‘‘out”? measurement 
the ballast capacitance, C,, remains fixed at C 
After the specimen bolder has been disconnected, 
Cz is measured by switching the bridge to “‘direct’ 
and balancing the bridge by changing the standard 
capacitor. The dial reading of the standard capac- 
itor is calibrated by measuring a known (';. The 
calibration should also include the 1 pF capacitance 
of the bridge terminals. 


4.3. Measurements at Ultra-Low Frequencies 


The Scheiber ultra-low-frequency bridge was used 
in the frequency range 10-?Hz to 210°Hz [I]. 
Although the bridge was designed for use with a 
three-terminal dielectric cell, we have found that it 
may be easily adapted for measurements with a 
two-terminal cell. Ordinarily the leads from the 
bridge to the electrodes are shielded at ground 
potential and the shield encompasses the entire cell. 
For the two-terminal holder the electrode leads are 
shielded but the shield is abandoned the base of 
the holder where the electrodes connect. A shield 
around the entire cell was improvised by lowering an 
aluminum box over the holder and connecting the 
the box to cround. 

The substitution method of measurement with the 
electrode spacing fixed at ¢, for the ‘in’ and ‘‘out”’ 
bridge balance is used. Equations (22) and (23) of 
Scheiber’s paper [1] in conjunction with the airgap 
correction term of this paper were used to calculate 
the capacitance and conductance of the unknown. 

For measurements below 10 Hz the surface con- 
ductivity of the sample can become a problem. It 
can be eliminated by wiping and cleaning the edges 
of the sample and thereafter handling the sample 
with tweezers and by introducing dry nitrogen 
through one of the ports of the cell. 


4.4. Self-Resonant Measurements 


The two-terminal holder can be used up to 600 
MHz by shorting the terminals with an appropriate 
length conductor and exciting the resulting circuit 
at its self-resonant frequency. The Boonton model 
280A UHF @Q-meter was used a combined oscil- 
lator and detector. The oscillator was coupled to 
the holder through a variable attenuator by 
of a single-loop inductive probe. The signal was 
picked up with a capactive probe and rectified to 
give a d-c voltage which was applied to a voltemeter 
in the model 280A. The voltmeter is calibrated 


means 
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full and half-power levels which can be conveniently 
checked by switching 3 dB into the inductive probe. 
The oscillator frequency was measured with a 
Hewlett Packard model 524C frequency meter and 
540A transfer oscillator. 

Measurements were made by placing the specimen 
in the holder and tuning the circuit by adjusting the 
oscillator frequency. The frequency corresponding 
to the half-power points is then determined and 
recorded along with the electrode spacing. The 
specimen is then removed and the circuit retuned 
by adjus ting the holder electrode spacing. The 
frequencies corresponding to the half-power points 
are again measured and recorded. One or more 
additional measurements of electrode spacing and 
corresponding resonance frequency must be made in 
order to determine the change in capacitance with 
resonance frequency. 

Calculations are based on eqs (1) through (6), 
where the terminals in figures 2a and 2b are shorted. 
Thus we treat the problem as a simple series 
L,CrR, cireuit for which the dissipation factor D; 
and total capacitance Cy are 


D wl’, ie. 


w*L, 


One may obtain e’ from eq (11 
remembering that here Cy 


and "bed from eq 
Cy, and 


ADxy Dp. -D yy wl; Ew 


_ R uv 


Substituting (32) into eq (13) gives 


The ratio of 


the frequency separation of the 
half-power 


points to the resonance frequency is 
measured to obtain the dissipation factor. The 
total capacitance which cannot be measured directly 
can be obtained by differentiating eq (31) to obtain 


dC, 


Af yy (34) 


Hence we can write 


_ dy Af Ajo (Get Ce 
dfn 2¢ ( is 


g 


Since the only part of C7 which ec a with spacing, 
t, is the calibration capacitance C, of the holder 
can write 


one 


e dt dl } Af —Afo 
dfn dt ?/ OF 
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Finally 


— dt d¢ ' 
df dt 


Afi—Afho 
26,A,/t, 


t_—? Ces Ci, Od 2 a 
a ge CE 


5. Sample Data 


In order to illustrate the 
obtained with the dielectric specimen holder de- 
scribed in this paper, room temperature measure- 
ments were made on two specimens machined from 
the same 2.54-cm-diam poly(methyl methacrylate) 
rod into disks of average thicknesses 0.3150 em and 
0.1577 em. The disk faces were made flat and 
parallel to better than 10 uw and no contact electrodes 
were applied. Both specimens were kept at about 
23 °C at a relative humidity of 50 percent for 
several months before being measured. The measur- 
ing techniques used are those described in detail in 
the previous sections, and the measurements on the 
two specimens over the whole frequency range took 
about 10 hr. Calculations of e’ and e’’ were done 
by hand. The results of the measurements are 
shown in figure 4 which is a plot of the values of 
e’ and e’’ versus log frequency for the 0.3150-cm- 
thick specimen of poly(methyl methacrylate). 


results which can be 


The data for the two specimens measured were 
found to agree to within 1 percent for ¢’ and 5 
percent for e’’ at all frequencies, and with just a few 
exceptions, the agreement was within 0.3 percent 
for e’ and 2 percent for e’’ over the frequency range 
covered. In addition, the highest frequency values, 
which are the most crucial as far as the magnitude of 
the circuit correction is concerned, agree with each 
other and with measurements made using a re- 
entrant cavity [3] within 0.3 percent in ¢’ and 
percent in e’’. Additional measurements made on 
specimens of polychlorotrifluoroethylene and fused 
silica gave results which were also accurate to 
within the limits claimed in this paper except for the 
loss in fused silica which could not be detected at 
high frequency. 


6. Discussion of Measurement Errors 


The two-terminal method described in this paper 
was intended to provide values of ¢’ and e’’ over the 
entire range of applicable frequencies without im- 
posing rigid requirements on specimen preparation 
and cell construction. In addition to uncertainties 
in the quantities explicitly appearing in the working 
equations of the previous section, attention must be 
given to the assumptions and approximations upon 
which these equations are based, the errors due to 
which will appear in an error capacitance term C’. 

The error capacitance (in fig. 2 and eq (5)) in- 
cludes corrections to the ideal assumed capacitance 





FIGURE 4. The dielectric 


mercial poly(methyl methacrylate) as a 


constant (e’) and loss index (e’’ 
function ot 


of a 0.3150-cm-thick disk 
measured with the 


of com- 


fre quency as two- 


terminal dielectric spec zimen holder using various measuring equipment. 


resulting from distortions of the electric field be- 
tween the measuring electrodes because of (a) tilted, 
misalined or nonflat electrode surfaces, (b) the 
presence of an airgap, (c) departure of the specimen 
geometry from an ideal right circular cylinder, 
(d) electrode edge effects, ni (e) the presence of 
air between the electrodes. 

(a) The errors due to nonideal electrodes are all 
second order effects with the actual capacitance 


given by 
d ay 
€’ €o ak +0(5 ) | 


where A and ¢ are the electrode area and separation 


and A represents the variation in ¢ in the case of 
tilted or uneven electrodes or the distance between 
axes of the electrodes in the case of misalined elec- 
trodes. It is not difficult to keep mechanical errors 
to within 25 u (25 4=0.001 in.). With typical sepa- 
rations, ft, of 1000 uw these errors amount to roughly 

0.1 percent. In addition, the tilt error increases 
the capacitance whereas the misalinement error 
decreases it. Also since each determination of C 
involves a difference between two measurements, the 
effects tend to cancel. Thus the mechanical con- 
struction need not require special attention. 

(b) The airgap causes some distortion of the field 
at the edge of the specimen (which increases with 
gap width and dielectric constant) and makes the 
actual capacitance higher than its uncorrected value. 
This in turn leads to a high value for the specimen 
capacitance C,. An upper limit to this error (e’ = ) 
for a typical measurement was calculated to be 0.2 
percent [6], and one can reasonably conclude that in 
most cases the error from this source is <0.1 percent. 

(c) Geometric imperfections in the specimen also 
lead to distortions. Rounded corners would cause 
an error similar to the airgap error. Uneven or non- 
parallel specimen surfaces lead to an actual capaci- 
tance greater than that assumed. The airgap capaci- 
tance is particularly affected by nonideal specimen 
surfaces. The specimen surface may well be tilted 


with respect to the adjacent electrode surface so that 
the airgap is a wedge of average thickness f[=15 yp 
and variation A=10 pu. In this case the assumed 
airgap capacitance may be in error by 20 percent. 

(This calculation is based on the worst the 
pheno ea of two concentric but mutually tilted 
circular disks in the limit of minimum separation.) 
Fortunately, the airgap capacitance enters the calcu- 
lations in a very insensitive way SO that an error of 
20 percent in C ‘ would result in an error of perhaps 
0.5 percent in C,, and the error would be such as to 
increase the value of C, above its true value. The 
geometry of the specimen does not appear to be 
particularly critical and special molding and machin- 
ing techniques are not required. 

(d) Electrode edge effects have been studied 
extensively [5]. With the specimen in place the 
edge capacitance is higher than one assumes and 
hence the calculated value of C, is too high. In 
general, this error increases with e’ and also with the 
proximity of the specimen to the edge of the electrode. 
If the distance between the specimen and the elec- 
trode edges is twice the separation of the electrodes 
then this error is <.0.1 percent [5]. Thus it is ad- 
visable to design the holder accordingly keeping in 
mind that measurement sensitivity decreases as the 
ratio of electrode area to specimen area increases, 
i.e., as C becomes a smaller part of the total 
capacitance. 


case 


(e) Since ¢’ for air is about 0.06 percent higher 
than assumed, then the actual capacitance with 
specimen removed is higher than calculated, thus 
lowering «’ by <0.1 percent. The effect on the 
airgap is negligible. 

A final approximation of the accuracy of the meas- 
urements can be obtained from an estimate of the 
accuracy of the measured quantities in the working 
equations. Generally, the calibration capacitances 
can be measured to the nearest 0.1 percent, the 
linear dimensions to the nearest micron and the area 
of the sample using a traveling microscope to 0.02 
percent. Considering eq (11) for e’ the factors and 
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their percent uncertainty are C, 
—C.+4 6A, /t,—C l, 
tion term, 0.5 percent. 


+0.3 percent; [(C., 
+0.3 percent; the airgap correc- 


In eqs (22), (23), 


(27), and 


(29) for e’’, the factors and their percent uncertainty 


I 
ara i 9) rea } 2 rea . ° 
are C';, +0.2 percent; eA,/f;, 40.3 percent; 4(5) 
D,), 


} 9¢ 
t (22 


+3 percent; Cy./(Cr—C.), 40.3 percent; (D 
+ (0.5% +0.00005) for GR 1615A; (D;—D,), 
+0.0005) for GR 716C; and the airgap correct- 
ion 1 percent. For eq (37) the percent uncer- 
tainty of the terms are dt d fas, 0.5 percent; dC./dt, 

1 percent; and Af; and Af,, + (0.2%+10kHz). 

The total percent uncertainty for ¢’ and ¢’’ is 
obtained from the square root of the sum of the 
squares of the uncertainties in the working equations. 
For «’ the percent uncertainty is +0.7 percent over 
the entire frequency range 10-? to 610° Hz. For 
e’’ we bave (1.5%+0.0005) for the ultra low 
frequency bridge, +(1.8°>+0.0005) for the GR 
716C bridge, 0.6°>+0.0001) for the GR 1615A 
bridge, (4°>+0.0005) for the 260A and 190A 
(-meters, and + (1.69%+0.0005) for the self-resonant 
measurements. 
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Some defects in a set of fused silica dielectric capacitors constructed in 1961 are listed, 


and methods for their 
set of capacit 


presented indicate that 


ors completed in 1964 is 


ts in 10° in five 


was less than tv Dp: 


1. Introduction 


as the dielectric 
investigation at NBS and else- 
where [1] * for many years. Early indications were 
that suitably constructed capacitors of this type 
would be at least as stable as air capacitors, and addi- 
tionally would be sufficiently rugged to withstand 
shipment between laboratories. 


Capacitors utilizing fused silica 


have been under 


A set of four fused silica dielectric 10-pF capacitors 
was constructed at NBS in 1961 [2]. The capacitance 
elements of these capacitors consisted of fused silica 
disks 7 cm in diameter and 1 em thick. Silver elec- 
trodes were fired on the two faces of the disks and a 
third electrode was fired on the cylindrical surface. 
The cell in which the element was mounted shielded 
the face electrodes from each other in the region 
external to the fused silica element, and was con- 
nected electrically to the cylindrical electrode of the 
element. The housing contained a four-terminal 
resistance thermometer for temperature measure- 
ment and was hermetically sealed. 

One of the 1961 capacitors was successfully hand- 
carried to several distant laboratories and returned. 
These laboratories included NRC (Canada), NPL 
(England), and NBS Boulder. The largest change 
in capacitance observed upon completion of the round 
trip was 0.2 parts in 10° (ppm). A shipment of this 
capacitor by commercial carrier was not successful 
and resulted in a change in its capacitance of 36 
ppm. This change was traced to a sensitivity to 
severe mechanical shock. 

The 1961 capacitors had a number of undesirable 
features, the shock sensitivity being by far the most 
important. Further investigation indicated that all 
of the capacitors in the set suffered from this flaw, 
and that the change in capacitance was caused by 
motion of the fused silica element in its cell. The 
element was supported in the cell by means of 
phosphor bronze springs which permitted motion of 
the element along the cell axis. The element could 
not return to its equilibrium position after displace- 


‘This work was supported in part by the 
Bureau of Naval Weapons, Pomona, Calif 
? Figures in brackets indicate the literature references at the end of this paper. 


Metrology Engineering Center, 


reduction or elimination 
deseribed in 
the typical drift of the 
months. 


are ce scribed. The 
detail. 


construction of a ne 
Preliminary stability 
1964 capacitors with respect to their m«¢ 


data 


ment because of frictional forces between the cylin- 
drical surface of the element and the spring contacts 
attached to the cell. This mechanism could not 
have caused a capacitance change if the direct elec- 
trical field between the face electrodes had been 
entirely within the fused silica. Unfortunately, the 
direct capacitance included a small contribution due 
to fields from the back of a face electrode which 
entered the fused silica element at the insulation 
gap between this face electrode and the evlindrical 
shield electrode, and terminated at the opposite 
face electrode. Motion of the element within its 
cell altered the leakage field and consequently 
changed the total direct capac itance. 

A second tiaw in the 1961 capacitors was a pro- 
nounced dependence of capacitance upon voltage. 
It was necessary to restrict the voltage applied 
these capacitors to 10 V rms to reduce the measure- 
ment uncertainty from this cause below 0.1 ppm. 
Since sufficient sensitivity is available at this voltage 
to see one. part in 10°, the voltage de ~pendence did not 
seriously limit measurement accuracy. An investi- 
gation into the mechanism of the voltage dependence 
was made, and it was found that the dependence 
could be almost eliminated by increasing the thick- 
ness of the fired silver electrodes. The evidence 
indicates that the existence of islands of silver on 
the dielectric surface which are not in good electrical 
contact with the body of the electrode is in some way 
responsible for the large voltage dependence. A 
similar effect has been observed in mica capacitors 
with fired or deposited silver electrodes [3]. 

Measurements of capacitance and dissipation fac- 
tor as a function of frequency indicated a capacitance 
difference from 50 Hz to 20 kHz of about 25 ppm, 
with a maximum slope and a maximum dissipation 
factor at about 1600 Hz. The shape of the curve 
suggested a relaxation mechanism involving polar 
impurities in the fused silica. An investigation of 
several — ‘al grades of fused silica disclosed that the 
use of a grade possessing exceptionally good ultra- 
violet pa ansmission characteristics resulted in capaci- 
tors with very much smaller frequency dependencies 
and with dissipation factors not greater than 3 107° 
throughout the audio-frequency region. 

Examination of the published optical characteris- 
tics for the various grades tested disclosed no correla- 
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tion between optical absorption and audio-frequency 
behavior. It is felt that the presence or absence of 
the audio-frequency relaxation phenomenon is an 
accidental result of the technique for producing 
fused silica. It is possible that changes in manu- 
facturing technique which have no effect on the 
optical properties of the fused silica may greatly 
affect the audio-frequency characteristics. 

The temperature coefficients of the 1961 capacitors 
were about 14 ppm/°C. The large values were be- 
lieved to be the result of a large temperature coeffi- 
cient of dielectric constant in fused silica. Resist- 
ance thermometers were built into the capacitor 
housings to eliminate the measurement uncertainties 
caused by unknown temperatures. 

A set of twelve 10-pF fused silica dielectric capaci- 
tors was built in 1964. Care was taken to reduce 
greatly all of the flaws listed above. The capaci- 
tance elements were constructed of type II Suprasil,® 
and were supported in their cells by means of tightly 
fitting polytetrafluoroethylene (PTFE) rings. Care 
was taken to keep the gaps between the face elec- 
trodes and the cylindrical shield electrode as small 
as seemed consistent with operation up to 200 V rms. 
The voltage dependence of each capacitor was meas- 
ured before final sealing, and the electrodes were 
replated when necessary. The 1964 capacitors are 
physically smaller and much lighter than the 1961 
capacitors, but were not designed with this as a 
principal objective. 

The construction of these new capacitors is de- 
scribed in some detail to allow their reproduction 
elsewhere. Some characteristics of the completed 
set are then presented. 


2. Construction 
2.1. Housing 


A cross section of the capacitance standard is 
shown in figure 1. The exterior of the assembly is 
constructed entirely of stainless steel. The hermeti- 
cally sealed region at the bottom containing the 
fused silica element is welded shut upon completion, 
and after baking, the container is filled with dry N, 
through a copper sealoff tube. Electrical connec- 
tions are made to the capacitor element and to the 
four-terminal resistance thermometer by means of 
glass-Kovar seals soft soldered to the 
compartment. 

Electrostatic shields not shown in figure 1 isolate 
the two face electrode leads from each other and from 
the resistance thermometer leads. The face elec- 
trode leads pass through stainless steel tubes as 
shown to a pair of coaxial connectors located above 
the liquid level of a constant-temperature oil bath. 
The resistance thermometer leads are attached to 
four binding posts located below the liquid level. 
Details of the construction can be seen in figure 2. 


sealed 


3 Available from Amersil Quartz Cumpany. 
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Figure 1. Cross section of the capacitance standard. 


Ficure 2. Details of the capacitor housing and superstructure. 

The dimensions of the fused silica element and its 
mounting are shown in figure 3. The PTFE radial 
support is inserted into the slot of the brass cell be- 
fore the inside diameter of the PTFE is finish 
machined. Each fused silica element is individually 
fitted to its cell by machining this radial support so 
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| that the element with its silver electrodes in place 
| can just be removed from the cell with a rubber 
suction cup. The axial (vertical) constraint on the 
element is controlled by the use of shims or by 
machining, so that an axial interference of about 
0.001 in. is obtained. This results in a slight com- 
| pression of the PTFE axial supports. 

Electrical connections to the face electrodes on the 
fused silica element are made by means of wires 
soldered to the phosphor bronze disk springs shown 
in figure 3. The disks are cut to the shape of an X 
and bent as shown in figure 4. Electrical connection 
between the cylindrical electrode and the cell is 
made by inserting three 0.030-in. diam phosphor 
bronze wires about 2 in. long into the circular gap 
between the cell and the element. The ends of 
these wires lie in the slot at the sides of the cell with 
the PTFE radial support. The three wires are 
positioned uniformly around the circumference of 
the cell and provide three points of electrical contact. 


2.2. Resistance Thermometer 


The four-terminal thermometer is located on the 
outside of the cell containing the fused silica element, 
as shown in figures 1 and 4. The cell is prepared 
by cementing a layer of linen cloth over it with 
shellac and baking at 100 °C. A leneth of B& S 

No. 36 Formvar insulated copper wire with a re- 
sistance of 25 2 is doubled over and wound bi- 
Ficgure3. Fused silica capaci- filarly over the linen. More shellac is then applied 
tanceelementand cell dimensions and baked. Leads to the current and _ potential 


terminals are connected to the resistance wire with 
soft solder. The junctions are tied down securely 


) 


and shellac is applied. The completed thermometer 


is baked for 24 hr at 75 °C. 


2.3. Fused Silica Element 


The fused silica element is ground to the dimensions 
shown in figure 3 with conventional toolroom grind- 
ing equipment. Some care is necessary to prevent 
chipping of the edges. After inspection for chips, 
the element is cleaned with alcohol, and then with 
soap and water. It is next immersed in an ultrasonic 
cleaning bath containing a wetting agent, and rinsed 
with water. Finally, it is immersed in an ultrasonic 
cleaning water bath and rinsed with distilled water. 
After drying at 100 °C, Dupont Silver Paint #4666 is 
sprayed onto one face and the cylindrical surface, 
and allowed to air dry. The other face is then 
sprayed and dried, and the silver is fired onto the 
element at 480 °C. It is found that single silver 
deposits thicker than 0.0005 in. are often defective, 
and that a single application is not sufficient to 
prevent a dependence of capacitance upon voltage. 
A satisfactory remedy is to buff the electrode surface 
with fine steel wool and apply two additional coats 
of paint to the electrodes. The preferred final 
electrode thickness is between 0.0010 and 0.0015 in. 
If the electrode thickness is increased to 0.002 in., 
the large thermal expansion coefficient of the silver 

Figure 4. Partially assembled capacitance element cell. often causes the fused silica element to chip when it is 


175 





cooled from the firing temperature. No noticeable 
chipping occurs when the silver thickness is less than 
0.0015 in. 

After the third electrode coat has been fired, 
edges of the element are very lightly 
separate the evlindrical electrode from the 
electrodes. This is done by hand with a diamond 
wheel and a special jig, using water as a lubricant. 
Care must “a taken at this point to avoid making the 
gap excessively wide and to remove all traces of 
silver from the beveled portions of the fused silica. 
If the diamond wheel is allowed to become con- 
taminated with silver, a thin layer of silver will 
remain on the fused silica, and a large dependence of 
capacitance upon voltage will result. 

The element is next fitted to its cell by machining 
the PTFE radial support. The element is then 
baked 24 hr in vacuum at 175 °C to eliminate 
adsorbed water, placed in its cell, and the cell is 
evacuated. The baking and evacuation are es- 
sential to produce a small dissipation factor and a 
stable c: ipacitance. 

After the capacitor has reached room temperature, 
it is measured and the adjustment necessary to pro- 
duce exactly 10 pF at 25 °C is estimated. Using the 
element and cell dimensions shown in figure 3, the 
capacitance will probably not depart from 10 pF by 
more than +0.05 percent. An adjustment as small 
as this can be made quite easily by cutting a ca\ ity 
into either a face or the cylindrical surface of the 
element with a small diamond wheel. The cavity is 
then sprayed with silver paint and refired. The 
adjustment sensitivity depends upon the diamond 
wheel diameter and thickness, the depth of cut, and 
the position of the cavity on the element. A jig may 
be readily set up to regulate depth of cut, and a 
calibration of capacitance change versus depth of cut 
measured with a trial element. Usually two or three 
successive adjustments suffice to produce a capaci- 
tance of 10 pF + 50 ppm. A typic al fused silica 
element after adjustment is shown in figure 5. 

After final adjustment, a test is made for voltage 
dependence, and if necessary, the entire element is 
replated. The stainless steel housing is then welded 
shut, and the system is baked at 65 °C with a vacuum 
pump attached to the sealoff tube. After 


the 
beveled to 
face 


e 


cooling, 


Ficure 5. Completed fused silica element with an adjustment 


of minus 50 ppm. 


the housing is filled with dry N» at atmospheric pres- 
sure, and the sealoff tube is pinched off and soldered. 
A series of temperature cyclings is made between 50 
and 0°C, and the unit is pli aced into operation. 


Figure 6 shows a completed capacitor. 


3. Performance 
3.1. Capacitance Stability 


All measurements of capacitance stability reported 
here were made at 1592 Hz using a 10:1 bridge, with 
all capacitors in an oil bath maintained at 25.0 °C 
The measurements seemed to indicate that these 
capacitors were at least as stable as any other 
capacitors in our laboratory. Accordingly, the mean 
of five of these ec: apac ‘itors, 109, 110, 112, 
and 113 was taken to be the basis for the comparisons. 
No temperature corrections were applied since the 

capac itors were at nearly the same temperature and 
had nearly identic ral temperature coefficients. The 
application of individual corrections for temperature 
slightly reduces the scatter in the measurement series 
reported below. Initially, all measurements 
made with 10 V rms applied. to the 
low ready comparison with the 1961 set of capacitors 
which had large voltage dependencies. Bridge 
readings were recorded under these conditions to the 


numbers 108, 


were 
capacitors to al- 


FIGURE 6. Complete d capacitance standard. 
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nearest part in 10°. Beginning November 9, 1964, 
the new capacitors were intercompared at 100 V 
rms. This allowed considerably greater precision, 
but prevented simple comparison with the 1961 fused 
silica capacitors. 

Capacitance differences between each capacitor in 
the set and the mean of the select five are listed in 


TABLE 1. 


table 1. The table demonstrates that most of the 
capacitors are stable with respect to each other, but 
they could possibly all be drifting at the same rate. 
Comparisons of the 1964 set of capacitors with the 
1961 set indicate no appreciable relative drifts, but 
conclusive proof of stability must await the comple- 
tion of a new calculable capacitor. 


Relative stability of eleven fused silica capacitors 


Decimal points indicate parts in 10° 





nh Nett 
FFFHFFSF 


non 
wmryeyyeyty 
F 
bo bo bo bo te 


1 


wets wwyow vs ty 


FFFEFFEFFFFFFFFF 


o—t 


otennwre 
! 
tober 
OO 


! 
t 


t 


Capacitor number 117 exhibits a fairly steady in- 
crease of capacitance with time. The electrodes of 
this capacitor were about 0.002 in. thick. This 
capacitor also has a large voltage dependence of 
capacitance, as will be seen later. It is believed that 
the diamond wheel used to bevel the edges of the 
element was contaminated with silver. 

3.2. Resistance Stability 
Measurements of each resistance thermometer 
are plotted in figure 7. The effect of variations in the 
mean oil bath temperature was eliminated from the 
data by plotting the differences between each resist- 
ance and the mean of the resistances of capacitors 
108, 109, 110, 112, and 113. 

Figure 7 indicates that the resistance thermometers 
are quite stable relative to each other, but shows a 
flaw in our oil bath. Between September 20 and 
October 10, a blast of air from a fan was directed at 
the southwest corner of the oil bath, where capacitor 
number 115 was located. This blast of air introduced 
a temperature gradient of about 0.02 °C into the 
bath and is believed responsible for the relative varia- 
tions observed on September 22 and October 5. 

The stabilities of the copper resistance thermome- 
ters were determined from a series of measurements 
made with a Mueller bridge. A calibrated platinum 
resistance thermometer in conjunction with measured 
data for the dependence of the copper resistances 
upon temperature was used to refer all measurements 








Capacitor 





Figure 7. Relative resistance thermometer stabil 


Insert shows relative positions in the oil bath. The southwest corr 


normally cold between September 20 and October 10. A resistans 
0.001 2 corresponds to a temperature change of about 0.01 °C, or ¢ 
change of about 0.1 ppm, 
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to a standard temperature. The limited information 
currently available indicates that the mean of the 
five selected resistances at a constant temperature 
has not changed more than 0.0001 Q in 2 months. 
This corresponds to about 0.001 °C; or, referred to a 
capacitance correction, about 0.01 ppm. 


3.3. Temperature Coefficients 


The temperatures of the 
may be measured by means of the internal 25-0 
resistance thermometers. The capacitances at a 
specific temperature may then be determined if the 
capacitance temperature coefficients are known. 

In practice, the temperature is usually treated as 
an implicit variable. A report of capacitance cali- 
bration then contains values for the capacitance and 
for the thermometer resistance at a temperature 
near 25 °C. A later measurement vielding a different 
thermometer resistance is corrected to the standard 
resistance thermometer value by reference to a plot 
of capacitance versus resistance for the standard. 

Measurements of the capacitance-resistance func- 
tion have been made over the temperature range 
20 °C to 30 °C. The curves deviate from a straight 
line by about 2 ppm at the extremes and are all con- 
cave upwards. The slopes at 25 °C range from 100.4 
ppm/2 to 110.0 ppm/Q. actual capacitance 


capacitance standards 


The 
temperature coefficients are about 10 ppm/°C which 
is substantially smaller than the 14 ppm/°C tempera- 
ture coefficients of the 1961 fused silica capacitors. 
No explanation for this difference has been found. 


3.4. Hysteresis 


The capacitances of most air or solid dielectric 
standards at a given temperature depend upon the 
previous temperatures to which the standards have 
been subjected. Table 2 shows the magnitude of 
this hysteresis effect in the 1964 set of capacitors. 
Prior to September 8 all of the capacitors were 
subjected to a temperature of 50 °C, and then 
replaced in the oil bath at 25 °C. On September 8, 
their capacitances were measured with respect to a 
capacitor maintained at 25 °C throughout the experi- 
ment. Following these measurements, the capaci- 
tors were cooled to 0 °C, and then placed back into 
the oil bath. They were measured a second time 
on September 9. Table 2 shows that the capaci- 
tances increased an average of 0.28 ppm following 
the 0 °C treatment. A second subjection to 50 °C 
and a return to 25 °C resulted in the September 11 
data in table 2. The capacitances returned most of 
the way to their original September 8 values, but an 
average capacitance increase of 0.06 ppm remained. 
Since the data were only recorded to the nearest part 
in 10’, this overall change may not be real. The 
effect illustrated by table 2 represents one of the 
most serious faults in the capacitors, and one for 
which the mechanism is not completely understood. 
It is believed that the large differential expansion 
coefficient between the fused silica element and the 
silver electrodes may be partially responsible. 





TABLE 2. Dependence of capacitance upon prior thermal 


history 

September 8 and September 11 data were taken at 25 °C after subjection to 50 °C; 
September 9 data were taken at 25 °C after subjection to 0 °C (see text). Capac- 
itance differences from the September 8 data are given in parts per million. 


Capacitor September 8 September 9 


September 11 


108 
109 
110 
112 
113 
114 
115 
116 
117 


C“mmenrkne 


Mean change 


va 


3.5. Voltage Dependence 


Some techniques and apparatus were developed 
recently at NBS for measuring the dependence of 
100 and 1000 pF capacitors upon voltage, with un- 
certainties in the order of 1 part in 10° [4]. The 
equipment was used for measuring the voltage 
dependencies of the 10 pF capacitors with slightly 
larger uncertainties. Table 3 shows the capacitance 
differences observed when the voltage was changed 
from 100 V to 200 V. The computed probable 
error of the data in table 3 is 6 parts in 10°, and 
the systematic error is estimated to be less than 2 
parts in 10°. Since the voltage dependence of a 
standard depends upon which electrode is at ground 
potential, all of the standards are marked with a 
“DPD” near the coaxial connector chosen for the 
ground potential lead. 


TABLE 3. Dependence of capacitance and dissipation factor 


upon voltage 
The increases in capacitance and dissipatior 


factor due to increasing 
from 100 V to 200 V rms at 1592 Hz aic shown 


the voltage 


3.6. Frequency Dependence and Dissipation Factor 


The result of comparing the fused silica 
with two 10 pF air capacitors at three frequencies 
is reported in table 4. Capacitor A was constructed 
with concentric cylindrical electrodes, and capacitor 
B with rigid parallel plates. Identical results were 
obtained at 159 Hz, using the two capacitors, but 
a significant difference was observed at 15900 Hz. 
The cause of the discrepancy is not known, but 
probably involves either mechanical resonances or 
excessive series inductances in one or both of the 


capacitors 


178 





air capacitors. The measurements at 159 Hz were 
extremely difficult to make due to poor bridge 
sensitivity. The computed probable error of the 
159 Hz data in table 4 is 21077. 


TABLE 4. Frequency dependence of the fused silica capacitors 
relative to two air capacitors (see text) 


Capacitor 159 Hz 1592 Hz 


15900 Hz A | 15900 Hz B 


ppm ppm 


—0.6 —1.9 
} —1.5 


—1, 


—2. 5 


9 


aa 
—3! 
—1.9 


Accurate measurements of dissipation factor can- 
not be made at NBS at this time. Comparisons of 
the fused silica dielectric capacitors with several 
types of air capacitors suggest that the dissipation 
factors of the fused silica capacitors are probably 
between 0 and +3107° at 1592 Hz. 


3.7. Shock Sensitivity 


The capacitors were tested upon completion for 
shock sensitivity by dropping them onto a hardwood 
table. Various angles of impact and various heights 
up to 8 in. were tried. The largest capacitance 


change observed was 1 part in 10° which occurred 
when a capacitor fell about 8 in. in an upright posi- 


tion. Falls through smaller distances usually pro- 
duced no observable change in capacitance. 


3.8. Transportation Experiences 


*oA fused silica dielectric capacitor designated num- 
ber 107 and which was identical with those described 
in this paper was shipped to the National Research 
Council of Canada (NRC) on July 28, 1964, via air 
parcel post. Prior to shipment, it was compared 
with capacitor number 108. Number 108 was hand- 





carried to NRC and returned between August 4 and 
August 17, 1964, and again between December 8 and 
December 11, 1964. Comparisons of the two ca- 
pacitors at NRC indicated no relative change in the 
two capacitances larger than 0.2 ppm; and as can 
be seen from table 1, no change in number 108 oc- 
curred during either of the two round trips. 


4. Conclusion 


Much more data is needed to evaluate conclusively 
the reliability of the 1964 set of fused silica dielectric 
capacitors. Existing evidence suggests that sta- 
bilities in the order of a few parts in 10’ may be ex- 
pected, and that the capacitors can probably with- 
stand normal handling during shipment between 
laboratories. 

Some further modifications of the design may re- 
sult in an improved standard. A study of the 
mechanism which produces the capacitance hysteresis 
effect with large temperature changes would be of 
value, and consideration might be given to improving 
the support system for the fused silica element. 

It is felt that the choice of 10 pF for a fused silica 
dielectric capacitor results in a design of nearly 
optimum stability. Other values would be useful 
as secondary standards and should not be difficult to 
construct. 
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Errors in the Series-Parallel Buildup of Four-Terminal 
Resistors 


Chester H. Page 
(April 20, 1965) 


The use of n equal resistors (a) in series and (b) in parallel provides an n?:1 ratio of 


potentially high accuracy. 
one-ohm standard to the hundred-ohm, 
Formulas are derived for the error in the 
and (b 
order error of the ratio. 
It seems feasible to construct a 1 
less than 1 in 108. 


[100 o} 


The use of n equal resistors (a) in series and (b) in 
parallel provides an n?:1 ratio of potentially high 
accuracy. Such devices are important for extending 
the use of the national 1-Q standard to the 100-Q, and 
thence to the 10 000-2 level. 

The first stage of this step-up involves the use of 
ten 10-ohm resistors, which must be four-terminal 
resistors to avoid large errors due to the connecting 
networks. 

Hamon! has described an arrangement 
terminal resistors permanently connected in series, 
and convertible to a parallel connection by adding 
jumpers. In that paper, he shows the use of com- 
pensating resistors in the potential leads to eliminate 
errors introduced by the added connections. 

Compensated lead ‘fans’? may be used 
potential terminals of the paralleled resistors, for the 
current terminals, for both. The aim of the 
present paper is to present a complete analysis of the 
general case, and formulas for the errors introduced 
by imperfectly compensated fans. 


of four- 


for the 


or 


lild-up resistor for the ¢ 
Dee. 1954 


ilibration of standard resistors, 
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Such devices are important for extending the 
and thence 
ratio, expressed (a 
in terms of first-order residual misadjustments which combine to yield the second- 


m buildup device with a 


use of the national 


level. 


in terms of design tolerance, 


to the 10 000-ohm, 


ratio uncertainty of 


We consider four-terminal resistors connected 
series by means of “tetrahedral” junctions, having 
the equivalent circuit shown in figure 1. Each june- 
tion supplies current and potential leads. This array 
can be converted to a parallel connection by adding 
four “terminal fans,” in figure 2. For analysis, 
the junction resistances can be considered as ab- 
sorbed in the fan-conductor resistances. The prob- 
lem is to make the four-terminal resistance of the 
combination precisely equal to R/n. If the resist- 
ances of the various arms of the current fan are 
adjusted to make each main resistor (2? Carry iden- 
tical current, the corresponding voltage drops will 
be identical and there will be no circulating current 
in the potential fans, and the potential across V;—V, 
will be independent of the resistances of the potential 
fans, and equal to JR/n. This requires that each 
arm of the current fan on the left have the same 
resistance, say r, except for the top and bottom arms, 
each of which feeds only one main resistor. 
end arms must have the 
right-hand current fan, each arm must have the 
same resistance, say 7’. The reciprocity theorem 
leads to the conclusion that if we use these com- 
pensated fans as potential fans, the four-terminal 


as 


These 
resistance 27. For the 
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resistance would be ?/n, independent of the current- 
fan resistances.:_;Since the potential fans can tolerate 
larger resistance than can the current fans, the 
compensation is usually made in the potential fans; 
resistances of several tenths of an ohm can be added 
to the arms, allowing adjustment to reasonable 
accuracy. 

Making either set of fans perfect yields zero error. 
This suggests that the overall error is in the nature 
of a product of fan errors, and suggests the possibility 
of reducing the effect of residual potential fan errors 
by making at least a rough adjustment of the current 
fan. For analysis, we consider four perfectly com- 
pensated fans as a nominal condition, with arbitrary 
maladjustments allowed in each arm. In addition 
each main resistor is allowed a departure from 
nominal. 

The complete circuit to be analyzed is shown in 
figure 3, where the labels indicate conductance rather 
than resistance. Capital letters indicate ‘‘average”’ 
values, 1.e., 


= =a, a 
24= 25 b,=Doe: 


Considering terminal G as the “ground” terminal, 
we have a network possessing three external nodes, 
and (n+1) internal nodes. A complete description 
of the network requires 3+(n+1) simultaneous 
equations, and we must eliminate the last n+1. 


If we impress currents J;, Jy;, Jy into the external 
nodes, and J;, J; . . . Zs; into the internal nodes, 
the voltages on these nodes are implied by the 
simultaneous equations: 


I; YirV1 T YinVu J 
Tn YuriV1 r 


YuVi+ + Yi niV n41 


+Yuaes Vins 


I, +1 


=VrnaiyVit VntingiV n+1 


(1) 


These equations are conveniently symbolized by the 
matric equation 


I=YV (2) 


where J and V are column vectors (each having n+4 
components), and Y is an (n+4) X(n+4) square 
matrix. 

Consider J;, Jy, and Jy; as components of an 
“external” current vector, J,, and J; . .. Ins, as 
components of an “internal” current vector, J;; 
similarly consider ‘‘external”’ and “internal” unknown 
voltage vectors, V, and V;. Grouping the terms 
appropriately: 


# = 
I; 
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FIGURE 3 


where a is a 33 matrix, B is a 3 (n+1) matrix, 6 
is its transpose, and y is an (n+1)*(n+1) matrix. 
In other words, the matrix is partitioned: 


a Bp 
[=| - ) V. 
B ¥ 


Since we are interested in knowing the external 
voltage with no impressed internal current, we must 
eliminate V, from 


(4) 


I,=aV.+BV; 
0=pV.+W, 
finding ; 
I,=(a—By"'8)V. 
V.=(a—By718)-1T,. 


This last equation can be written as 


Vi 
Vu 
Vin 


The required four-terminal resistance is the transfer 
resistance (V},;—Vy)/J; under the condition Jy= 


Inn=0, hence 


I; 
In | 
Tun 


The problem is to carry out all the indicated 
algebraic manipulations, with all resistors subject 
to arbitrary tolerances. This will be done in an 
appendix, using the following strategy. 

First, we carry out the manipulations required 
for (a—By~'8)- assuming the nominal values of the 
network elements; then, assuming that the de- 
partures from nominal are small, find the resulting 


, : 1 : : 
correction to the nominal Ro ( 7) as a series in 
na 


powers of the perturbations. This yields a formula 
for the fractional error as a sum of terms of the 
type gigs; gi, gibi, ete.; aib;, ed; The terms in- 
volving the g; (tolerances in the resistors of the 
series-parallel set) are ordinarily negligible com- 
pared with the pure measurement-network errors. 
These latter give (from eq (A38)): 


4 G Gy Dy. Ge Ds 
a ‘  — tes 
Tae aBtAB 


Error 
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Since Ya;—0, ete., this error arises from accidental 
correlation between errors of fan connections to 
common points. The worst case is where half the 
a; have maximum positive error, the other half 
having maximum negative error—and where these 
pair up with the same situation among the 6,. 
Normal design would have D=A, B=C, and 
A B. Now let the maximum fractional error 
in the voltage fans be e,; in the current fans, €;. 
he oi 
The resulting extreme error is 4 — e, 


A 


mQ, e 10, <,=10 


Reasonable 


l ] 
—=10 Q, l 
G A 
vielding a maximum error in the series-parallel 
ratio of 4107. Fractional errors of 10~° in the 
main resistors make the g; contributions completely 
negligible. 


estimates are 


The above computation is for maximum error. 
Random pairing of the fan errors gives a much 
smaller expected error. 

The next problem is to adjust the device after it 
has been constructed, and to evaluate the residual 
error. We note that if the current fans are perfect, 
there will be no voltage across any arm of the 
potential fans, hence any voltage across an arm of a 
potential fan is a measure of current-fan errors. 

The potential fans have relatively high resistance, 
so are conveniently ‘‘trimmed” on a resistance 
bridge. To adjust the current fans, we can connect 
a microvoltmeter across each potential fan resistor 
in turn. This essentially the deviation 
of the potentials on the internal nodes, from their 
average. The current-fan arms can be filed to 
reduce these voltages to approximately zero. 


measures 


In the analysis, we have treated the tetrahedral 
junctions as though their centers were available 
terminals. Any adjustment procedure must take 
cognizance of the fact that these terminals are 
fictitious; any physical connection is separated from 
the tetrahedron center by a low resistance. This 
could cause errors in direct measurement of current- 
fan resistance. The proposed scheme of measuring 
potential-fan resistances and voltage drops avoids 
this difficulty, since the potential fans are loaded 
with additional resistance. 

In appendix A, we evaluate these residual adjust- 
ment errors quantitatively. We find to first order 
that the voltage across the arm B-+-b; is given by 


(10) 


and similarly the voltage across A 
into terminal 


as, for current J 
Il) and out of terminal (ITT), is 


Comparison with eq (9) yields 





. if ( 
epoca ' (A+B) 


et cat XnsiWoss 
GT XWeAX Wet. . pet st) 


LX,.W,) }. 


(12) 


+(C+D)(X2W2+X W4 


For a test current of 1 A, the fan tolerances 
assumed before yield voltages of the order 


IX~2X10-V 
IW~2x10-7V 


with 
nGA~10 s 


Note that the larger voltage is generated by the 
current-fan error (under the assumed construction 
tolerances); this should make it feasible to adjust 
the current fan to the next order of magnitude. 

The most sensitive set of measurements for evalu- 
ating the residual error is probably the set of residual 
adjustment errors; ie., the residual potential-fan 
resistance differences and voltage drops. The ap- 
propriate formula is (from eq (A34 
Residual error 


2G} .X,6,/B+X3)3/B.. . 


a ee: Oe 


$ 


where positive XY; is associated with node 7 positive 
with respect to the appropriate potential terminal. 

Note that the measured error quantities are first- 
order effects in the network tolerances; the computed 
residual error is a second-order product of these 
terms, so can be reasonably accurately evaluated by 
this procedure. : 

There are, of course, additional errors not attrib- 
utable to the current and potential fans. The first 
of these arises from the fact that the series resistance 
is not nN, G, but is 


Seer 9 G) 


vielding a fractional error of only 107° for g,/@~10~ 
Another source of error lies in the imperfections of 
the tetrahedral junctions. Under conditions of 
suitable design and proper adjustment procedure, 
these junctions can be balanced to transfer resist- 
ances of less than 107° Q In a string of 10-Q 
resistors, this is then a potential source of error of 
the order of 107°. 

It seems feasible, therefore, to construct a 1:100 Q 
buildup device with a ratio uncertainty of less than 
1 in 10°. 
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Appendix A 


The submatrices of equation (4 











It is convenient to decompose these matrices into 


sums of nominal values and perturbations: 


B; 


where 
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We write the 
nominal value 


admittance matrix as the sum of a 


and a perturbation: 


on ae 
P+) B NTY 


YK 


(A11) | 


(A12) | 


where 
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This yields 











These matrices are readily verified by inspection of 
the circuit when the elements have their nominal 


| values, 


The matrix Z, involves the term (['+)~' standing 


| by itself; we shall eventually need this term and 
shall use a series expansion for it. 


The four-terminal resistance of interest is the 


| transfer resistance expressed as Vy;— Vi for J;=1, 


(A12a) 


(A12b) | 


+0) ~"BoZ1. 
Now r t 


combination 6, 


not 
' is, from appendix B: 


+| BE 0 Be oO ok 
a CF D4 sk 
where 


| the partitioning of Y. 
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readily found, but the | 


all other input currents being zero: 


? [(O 


1—1)(0 














Lod] 
In (A16), we have partitioned the vectors into 
“external” and “‘internal’’ vectors, corresponding to 
For the nominal case we find 


1 
1—-1)Z,1 0 
0 


(A17) 


| We are also interested in the transfer resistances that 


are measured as voltage drops across the potential- 


fan arms: 


l 
0 
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(oO—1 OQ O O.. 











All X,=0 for the nominal case. 

The perturbation from nominal is found by ex- 
panding Y~ in a series (as shown in appendix B for 
(T+70)7): 


| y-. Yi'-Ya"(: ie ¥: 
By € 
re A (; *) rs? (; *) Y5'+ 
Bi € By € 


to second order in the perturbation. 
X¥,-[0-1 0)(0 0 1 0 0))¥— The algebra will be alleviated by noting some 
—. special properties of the matrices 8, y;, and y2: 


X»=([(0 0—-1)(0 1 0...0))¥= 
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Note that premultiplication of either y, ory, by either | These simplify the multiplications required by eq 
(1 01 0...) or (0 1 0 1...), followed by | (A19): 
postmultiplication by either 

Zo jo (: ~ ( 0 246:+Zie 
@ (0 ~ > , - we ¢ , 
7 4 Za By € ZB 2181 +Z2e 
0 ] 
1 0 ; 
of], results in a null vector. 


| Note that 
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| so that 


Thus from (A12b) and (A13) in conjunction with 0 B 0 MZ, 
(A20) and (A21), we have the very useful relations: ¥:*6., % eee is Mee 
By € Z.M Z2B, 2, +Z18:Z.+Z.€Z> 
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and finally 
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Neglecting g/G relative to 6/B and ¢/C, we have 


l 2 : P ‘ 
i. a Oe & >. 
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nG 4 


A 19 1 
for evaluation of error in terms of measured residuals. 


To express /? in terms of circuit element tolerances, 
we ! eed to ¢ xpress A: 


b 1 0 O)M(T+4 
1 O 0O){ 48, +7119 


Wwihuch becomes 








2K) 


and the successive terms of (A386) diminish as powers 
of GE and GF, which are small compared to unity. 
Hence the first term of (A36) is sufficient for our 


purposes, for use in (A385 
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A power series expansion yields then 
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Centerable Rotator for Measuring Properties of Crystals 
Charles P. Saylor and Herbert B. Lowey 


(February 


23, 1965) 


By allowing the unsupported portion of a spindle to be relatively long so as to extend from 
the rim of the microscope stage all the way to the axis, simple controls can be provided for 


centering a crystal with respect to a horizontal axis of rotation. 
described, centration can be assured within about a micron. 


In the device which is 
The angle of rotation with 


respect to a fixed position can be read with an error not greater than 0.05 degree. 


1. Necessary Adjustments With the Single 
Axis Method 


In the single axis of rotation method of measuring 
optical properties {1],! a crystal is held upon a spindle, 
it is rotated about an axis parallel to the stage of the 
microscope until one after another of the three char- 
acteristic vibration directions is made parallel to the 
stage, and each corresponding index of refraction of 
the crystal is then matched with that of an immersing 
liquid of known or subsequently determined refrac- 
tivity. A crucial activity is the coordinated control 
of crystal orientation and polarizer vibration direc- 
tion so as to achieve and recognize the three sig- 
nificant vibration directions. In this action, while 
the spindle is rotated slowly with one hand, the other 
hand, on the of information derived from 
changes in the interference figure, executes a con- 
trolled rotation of crossed polarizer and analyzer. 
The hand on the nicols sustains the condition by 
which an isogyre * will always cross the center of the 
interference figure. 

Devices for rotating the spindles have been 
described by Wood and Ayliffe [2], Bernal and 
Carlyle [3], Wilcox [4], Hartshorne and Swift [5], 
Hartshorne [6], and others. Well planned as these 
rotators are for relatively large crystals, they 
imperfectly cover one condition, important with 
crystals which cover less than the full microscopical 
field. If the interference figure of a small crystal is 
to be sufficiently contrasty, its image must be isolated 
from other light which might pass to the observer. 
This is best accomplished with a field diaphragm 
placed above the Bertrand lens at the level where the 
image of the crystal itself is formed or, if a Klein 
magnifier is used, by a diaphragm in the image plane 
of the eyepiece. Under these conditions, if the 
crystal is uncentered with respect to the axis of 
rotation, it is necessary continuously to adjust the 
crystal vertically or horizontally as its image strays 
from the position of the field diaphragm. Both 
hands and the mind, however, are already occupied 
with the control of two rotations. While foot 


basis 


1 Figures in brackets indicate the literature references at the end of this paper. 

2 A locus of directions in a crystal (as seen in orthographic projection) for which 
the projected vibration directions, whether trivial or characteristic, are each paral- 
lel to the vibration directions of a polarizer or analyzer 
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regulated motors might be installed, these would not 
overcome the problems that arise from division of 
attention. The effective solution is a rotator that 
permits centration of the crystal with respect to the 
axis of its rotation. 


2. Design of a Centerable Rotator 


In designing the apparatus which is described here, 
we found it expedient to depart from preceding prac- 
tice to the extent of placing the entire rotator beyond 
the stage of the microscope. The actual support for 
the crystal, in the form of a small glass tube, thus 
extends unsupported across the radius of the stage. 
Because of the elasticity of glass and the smallness of 
the forces, this introduces no difficulty. Sensitive 
centering controls are easily achieved. In the 
device here described these operate at positions 90° 
from each other. The adjusting of one control does 
not upset the centering that has been achieved with 
the other. 

The centerable rotator which we are now using is 
illustrated in figure 1. The stationary tube of the 


ie aca So 


Shite 


FiGuRE 1. The centerable rotator. 

The back and bottom nobs drive two jaws by reason of the difference between 
28 and 32 thread screws. The others are spring loaded. The front and back 
jaws, cut away over half their diameters, lap over the top and bottom jaws. 
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rotator, normally held in the chuck of a micro- 
manipulator, is attached to an end plate with a 


20-interval vernier extending across 38°. At the 


right end, a thumb screw compresses a 2 mm long 


annulus of “‘microchemical” rubber tubing, 1.5 mm 
in internal and 3.5 mm in external diameter. This 
serves as a flexible chuck to hold and center the end 
of the tube farthest from the crystal. The hollow 
shaft has an internal diameter of about 3 mm. 
Since the device was planned for use with 1.35+0.15 
mm glass tubing a total excursion 
of about 14; mm at the controlled end is possible. 
The overhang of the tubing (the distance from the 
is not less than length 
of the rotator so that the range of adjustment at the 
crystal is about 3mm. In the drum at the left are 
the controllers, those with jaws driven by differential 
screws at the back and bottom and those with jaws 
pushed by light springs at the front and top. The 
sets which are separated by 90° overlap like the 
parts of a lap joint. 


a standard size), 


controlled end to the ervstal 





FIGURE 2, 


Details of the centering controller are shown in 
figure 2. For movement of two jaws, a }4-in.-28 
machine screw is drilled along the center and tapped 
with a No. 8-32 thread. The actual jaws are cut 
into the central screws, one side of each of these being 
cut away for 11 mm. They thus glide in a channel 
outside the point where they overlap the other set. 
This keeps them from rotating. The other two 
jaws are spring loaded. The differential movement 
is such that the jaw advances a little more than 0.1 
mm for each full turn of the larger screw. This 
corresponds to about 0.2 mm at the crystal, thus 
providing all necessary sensitivity. By comparison, 
one full turn of the fine adjustment of the standard 
microscope produces a movement of 0.1 mim. 

In order to insert the glass tube into the rotator, it 
is convenient to have means of retracting both the 
screw driven and the spring loaded jaws. As we 
have chosen to manage this, the spring-loaded jaws 
are momentarily held back by a flat collar in the back 
of the jaws with two camlike wedges that can press 
the jaws outward. 


De tails of the rotator. 


ifromend. (b) Differential screw drive. 
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The graduated drum provides an easy means for 
recording positions of the spindle. Since the drum 
is only 5 cm in diameter, markings at 1 deg intervals 
would be separated by less than 0.5 mm. This 
would be too close. When 2 deg marks are combined 
with the 20 interval vernier there is adequate 
accuracy. Thus the position of the drum as seen 
in figure 1 clearly lies between 213.5 and 213.6°, 
being the former. Obviously the un- 
certainties are less than 0.05 deg. Since the accuracy 
of settings of the crystal is rarely better than 1 deg, 
this precision of measurement meets all requirements. 


closer to 


3. Supporting the Crystals 


It is well to insert a word about the use that is 
made here of glass tubing for the spindle. Equip- 
ment is available for drawing tips on capillary glass 
tubing with practically perfect control. By having a 
sharply tapered end on such a tubular spindle, 
crystals can be held in either of two satisfactory 
ways. This can be by suction applied to the tubing 
by means of the nipple shown at the right, figure 1. 
Although some immersion fluid leaks between the 
crystal and the end of the tube, the steadily in- 
creasing diameter prevents viscosity from 
destroy ine the force at the attachment. Occasionally 
the leakage is so great that too much of the immersion 
compound will be lost. 

The technique, which is also limited to 
tubes, makes use of drying cements. If a tapered 
very fine rod is inserted into a cement and with- 
drawn, surface tension tends to remove the cement 
down to such thinness that the residue cannot provide 
With a somewhat larger rod, the residual 
the sides, but not to the end. 


lk msses 


second 


adhesion. 


cement adheres Lo 


With a tapered tube, however, the outside is left 
nearly bare whereas the cement flowing into the 
canal forms a nearly plane surface across the end. 
When this plane end is touched to a crystal, the 
crystal is caught up by capillarity. As the cement 
dries in a circle of contact between the crystal and 
the tube, the attachment becomes very strong. 
Of the many different cements which can be used, 
one of the most satisfactory is ordinary corn syrup. 

By either of these methods of attachment, a 
crystal will hold an unvarying position at the end of a 
tubular spindle. The attachment endures through 
the necessary reorientations of the crystal, changes 
in composition of the immersing fluid, and determina- 
tions of the refractive index of the fluid. Since the 
correctness of orientation of the crystals and their 
vibration directions can be found with adequate 
sensitivity, and since the accuracy of a refractive 
index measurement is limited only by the criterion 
of match and the patience of the observer, the 
rotator which is here described facilitates highly 
accurate determinations of the principal optical 
properties of cry stals. 
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Equipment for Single-Crystal Growth From the Melt 
Suitable for Substances With a Low Melting Point 


Avery T. Horton and Augustus R. Glasgow 


(April 27, 1965) 


Glass apparatus is described for the purification and growth of crystals with melting 


points up to about 130 °C by 


perature drop exists 
the interface. 
with the glass container. 
men to solidify 


oxalate, benzoie acid, and benzene. 
radiation stability of pure substances. 


1. Introduction 


The procedures for the growth of a single crystal 
from the melt commonly restrict growth to an iso- 
thermal growth surface held in a horizontal plane and 
perpendicular to the flow of heat from the melt [1]. 
The surface is not permitted to regress (remelt). In 
the Bridgman method the growth of a single crystal 
takes place at the boundary between two tempera- 
ture zones, above and below the melting point of 
the substance. The temperature boundary is nor- 
mally produced by arranging vertically two furnaces 
which are separated by a thin laver of heat-insulating 
material with an opening for passage of the container 
in which the material is held. Crystal growth pro- 
ceeds by relative movement of container and 
furnaces. 

Too often the temperature gradient at the boundary 
of the two zones is neither sharp nor well-defined. 
Convection currents of air in and between the zones 
cause nonuniform temper: atures at the boundary and 
in the two zones. Consequently the uniform growth 
surface is not realized. 

Following an observation on the remarkably low 
transfer of heat across the interface between two 
immiscible liquids, equipment [2] for the growth of 
single crystals from the melt has been developed to 
use the very sharp, uniform, easily controllable tem- 
perature gradient thus produced. ‘This type of equip- 
ment is limited to temperature ranges in which 
liquids are thermally stable, do not react with con- 
tainers, and in which liquid pairs remain immiscible. 
With each liquid heated independently at the wall of 
a container, thermal or mechanical stirring in each 
layer can be expected to maintain two uniform tem- 
perature zones above and below the melting point 


1 Figures in brackets indicate the literature references at the end of this p iper. 


a modified Bridgman technique. The 
melt and the solid are maintained by two immiscible liquids at whose interface 


also be 


temperatures of the 
a sharp tem- 


A sealed Pyrex tube containing the specimen is slowly lowered through 
The specimen is maintained under 


its own vi ‘por pressure in contact only 


Provision is made for the removal of the last portions of the speci- 
from each of several successive 
rhe equipment has been found effective for the purification of trans-stilbene 
It could 


rece rystallizations. 
dimethyl 


used for careful studies of f the ‘rmal and 


of the substance being crystallized. Good crystals 
are produced only when the lines of heat flow from 
the melt are substantially perpendicular to the grow- 
ing crystal surface. This means that the growth proc- 
ess must be slow enough for most of the heat to be 
dissipated through the erystal rather than through 
the container walls near the interface. A convex 
solid surface is indicative of satisfactory heat flow 
and minimizes the effects of unwanted nucleation at 
the wall of the growth tube. 

With all-glass construction, visibility is unimpaired 
to observe nucleation, “lone seed survival” in the 
early stages of growth, and the entire growth process. 
A specially designed growth tube provides for the 
introduction of material by distillation or sublima- 
tion in a vacuum, subsequent sealing, and the growth 
of single-crystal material in the presence of its own 
vapor only. Appended seal-off bulbs allow the iso- 
lation and removal of the last portion of melt (or 
the last portion of crystal to grow) in which im- 
purities have been concentr: se by continuous rejec- 
tion into the molten phase at the crystal-melt 
interface. 


2. Apparatus 
2.1. Growth Tube 


The growth tube, shown in figure 1, is used as 
follows: The material to be crystallized is introduced 
at A by vacuum distillation or sublimation, with 
subsequent sealing at S,. The tube is attached by 
eyelets C and E to cables that control its movement 
through the immiscible liquid bath by means of 

synchronous motor. Crystallization is initiated - 
the tip of the capillary portion D and allowed to 
progress through the length of the capillary with 
the object of promoting the selection of a single 
oriented crystal at its juncture with the main body 
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of the tube; crystallization in the body of the tube 


progresses in the form of one or, at most, a few 
cry stals as the tube passes continuously through the 
interface between the two liquids. The last (im- 
pure) portions V,, V2, V3 of the crystal to form at 
each successive pass are remelted and transferred 
into the attached bulbs (*V,, *V., *V;) by inverting 
the entire growth tube and sealing off at S., 53, or 54. 
The purified material is transferred to evacuated 
containers through the break-tip at B. Typically 
the growth tube contains about 500 cm?’ of liquid 
substance, the attached bulbs *y, etc., have a volume 
of 10 em’ each. Growth tubes of a variety of sizes, 
but similarly proportioned, have been used. Various 
shapes of the capillary portion D have been tried 
with varying success for different substances. 


2.2. Dithermal Jacket 


Suitable immiscible liquid systems which have 
been studied include paraffinic oil (Marcol)-ethylene 
glycol, silicone sil-ethylene glycol, toluene-water, 
and water-carbon tetrachloride. Good boundary 
conditions were obtained between the two zones in 
all cases, provided that the temperature of the 
liquid in the lower layer was 20 to 30 °C below its 
boiling point. Silicone oil (Dow-Corning 500) and 


outh appa atus. 


ethy lene glycol had the widest range. After about 
6 months as bath liquids at approximately 150 and 
100 =. respectively, the silicone oil remained color- 
less, but the ethy lene oly col tended to become cloudy 
and was generally replaced. 

A diagram of the complete apparatus is shown in 
figure 2. The immiscible liquids, L-1 and L-2, are 
heated in the glass column by resistance-wire heaters 
H-1 and H-2. Thermocouples for measuring tem- 
peratures hear the interface, I, are connected to a 
selector switch Ss. They are located in six re-entrant 
wells, TW, placed 0, 1, and 2 cm above and below 
the interface. ‘Typical temperatures in the dithermal 
jacket, observed at these six positions during the 
growth of trans-stilbene, were in order of decreasing 
height 157.3, 155.5, 143.7, 112.7, 100.1, and 99.2°, 
respectively. Hence, the temperature difference 
across the interface was 31 °C. The growth tube and 
weight, GW, are centered and lowered in the column 
by the cable, C, from a drum which is driven by the 
hour shaft of a synchronous clock motor, M-—2, 
which in turn is controlled by a current interrupter. 
Speeds down to!,em per dayhave been used. About 2 
weeks are required to complete one crystallization 
cycle. The growth tube is placed in or removed 
from the column by detaching the inlet tube, IL-1, 
and the water condenser, WC, at the stopper, ST. 
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During crystal growth the upper liquid was generally 
held 25 °C above and the lower liquid 25 °C below 
the melting point of the sample. It is suspected that 
experimentation with smaller temperature differ- 
ences will prove interesting; the lower thermodynamic 
driving foree and the reduced thermal shock may 
tend to lead to the production of more perfect 
crystals. 

Since the descending growth tube displaces first 
one, then the other, of the two immiscible liquids, 
provision is made to maintain the interface in fixed 
position. The decrease in volume in the upper liquid 
chamber is compensated by the addition of fresh 
liquid from the reservoir, L-1, R, through the pump, 
HS, and inlet tube IL-1; excess liquid overflows back 
to the reservoir. The pump is made from a hypo- 
dermic syringe, HS, and glass check valves, V, and is 
driven by a reciprocating device, M-1. The bearing, 
SB, serves to keep the piston rod centered. Liquid 
displaced from the lower laver by the descending 
tube overflows from the leveling bulb, LB, at OL 
and collects in the flask, L-2, R, from which the bulb 
LB can be manually replenished at the end of a run. 


3. Notes on Experimental Procedures 


In experiments in which the aim is to preserve 
the single crystal intact, the entire sample is per- 
mitted to crystallize and the assembly is brought 
slowly to room temperature. Most crystals have a 
higher coefficient of expansion than glass and will 
therefore pull away from the container walls. The 
growth tube is parted by scoring with a glass cutter, 
and the crystal is removed intact. 

In experiments in which the aim is to acheive a 
high degree of purification, crystallization can be 
stopped before the last, most impure portion has 
solidified. The erystallizer is quickly withdrawn 
from the bath, inverted, and the liquid collected in 
one of the small bulbs V (fig. 1) which, after solidi- 
fication by cooling, is sealed off and rejected. The 
sample is not allowed to cool to room temperature 
whenever another crystallization cycle is planned. 
Otherwise on reheating there is danger of cracking 
the glass container. If reheating is necessary after 
cooling, the best procedure is to heat the glass rapidly 
to permit melting of the surface regions of the solid. 
The liquid so formed will yield during subsequent 
expansion of the remaining sample. 

Another procedure which has been employed is to 
allow the entire sample to crystallize. The tube is 
now removed and the temperatures of the two 
immiscible liquids are reversed, the upper liquid is 
lowered to about 25° C below the melting point of 
the ree and the lower liquid is raised to about 

‘above the melting point of the sample. The 
alk tube is now reinserted in an inverted position 
and the guide weight and support cable attached to 
opposite eyelets. The crystal is placed in the upper 
bath, allowed to come to temperature, then lowered 
to a position in the interface determined by the 
amount of material to be rejected. This portion 


is melted and drained into one of the bulbs V (fig. 1), 
and the tube and contents are removed from the 
column. The liquid in V, is frozen, cooled low 
temperature, and the bulb is sealed off at Sy. 

The growth cycle is repeated until purity assays of 
small portions last to crystallize are of a satisfactory 
level of purity. The purified substance may then be 
transferred to other vessels through the break-tip, B. 

For a given sample some indication of the impurity 
is obtained by observing the growth process, particu- 
larly the shape of the solid-liquid interface. When 
the contamination is considerable the solid tends to 
be polverystalline especially in a typically conical 
core at the top of the solidified portion. Several 
crystallization cycles lead to better single-crystal 
growth. In the final growth of a very-high-purity 
sample, the solid-liquid interface tends to have a 
convex shape. 


4. Substances Purified by 
Growth 


Single-Crystal 


4.1. trans-Stilbene 


Commercial trans-stilbene — (scintillation-grade, 
about 99.9 mole percent, mp 124 °C) was melted, 
filtered, vacuum-sublimed twice, and a 200 cm® 
portion transferred to the growth tube. The tube 
was lowered through the interface between upper and 
lower bath liquids at 155 and 100 °C 
Four passes, with removal of 5 em® last portions 
after the first and second pass, were required to 
produce a single crystal which was optically clear. 

The contour of the crystal surface changed as 
impurities were removed in successive recrystalliza- 
tions. The first pass resulted in a shallow V-shaped 
surface containing colored material which remained 
liquid in the bath at 100 °C and was shown to have a 
melting point 24° below that of trans-stilbene. A 
marked improvement in the surface was noted 
during the second cycle. The third growth produced 
a specimen which was still not optically clear, but 
the fourth run resulted in a clear single crystal. An 
analysis of the purity of the final product was not 
made. 


respectiy ely. 


4.2. Dimethyl Oxalate 


Commercial dimethyl oxalate (93.7 mole percent) 


was purified by two vacuum sublimations to a 
purity of 99.6 mole percent measured cryoscopically 
[3]. During sublimation methyl alcohol was removed 
in a volatile fraction, while methyl hydrogen oxalate 
and oxalic acid remained in the residue. 

Two samples of the sublimed dimethyl oxalate 
(mp 54 °C) were grown as single crystals by lowering 
the growth tube through the interface between an 
upper liquid at SO °C and a lower ~— at 30 
°C. The first sample (about 250 em‘), erys- 
tallized twice with one intermittent al 
of a small last portion, gave a single crystal 
which was perfectly clear. The second sample 
(about 400 em*) was crystallized three times 
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with two impure last portions removed. The 
final surface to crystallize was flat and the crystal 
was clear. An analysis of the purity of the final 
product was not made. 


4.3. Benzoic Acid 


The starting material was benzoic acid of purity 
99.997 mole percent (mp 122.37 °C) as determined 
by comparison with a thermometric cell of benzoic 
acid of known purity. The sample was degassed in a 
vacuum to remove water vapor and air, then vacuum 
distilled into the growth tube at 150 °C. The tube 
was lowered through the interface between liquids at 
150 and 100 °C, 15 em? of the last portion to crystal- 
lize was removed, and the crystal was melted and 
regrown. The final specimen was not a_ single 
crystal; half of the material appeared to be a single 
crystal, while the other portion appeared to contain 
two crystals. 

The specimen was taken from the tube and cut 
with a jeweler’s saw to isolate the single-crystal por- 
tion. All imperfect portions were removed by cleav- 
ing. The specimen was sealed in a platinum en- 
velope, gold-soldered to prevent contamination, and 
was used as an acidimetric standard of high purity in 
comparison studies by Bates and Wichers [4]. 

In the latter study an intercomparison by precise 
differential titrimetry corroborated the above purity 
of the starting material (determined independently 
by calorimetric analysis), and showed the single 
cry er benzoic acid and single crystal potassium acid 
phthalate (grown from aqueous solution) to be not 
less than 99.9998 mole percent pure. 


4.4. Benzene 

Benzene (fp 5.5 °C, purity 99.8+0.2 percent) was 
purified by single crystal growth from the melt with 
the upper liquid at room temperature, the lower 
liquid a, to about —28 °C by immersion of the 
lower half the column in a circulating bath of 
cious. “The purified benzene, used in compara- 
tive tests of methods for the determination of purity 
[5], contained impurity believed to be less than could 
be disclosed by thermal analysis (<(10-* mole per- 


cent). Actual cryometric and calorimetric assays 
showed a purity of 99.995 mole percent, but con- 
tamination by water adsorbed on the inner surfaces of 
the measuring apparatus was believed to have 
contributed to this value. 


5. Further Potential Applications 


The immiscible liquid-pair system may operate at 
much higher temperatures than have been employed 
in the present work. Thermal stability and im- 
miscibility of liquids are major requirements. 
Molten salts, metals, and metal alloys could con- 
ceivably be used, although alterations in heating 
methods, insulation, vessel and column construction, 
ete., would be involved. Such developments have 
not been pursued, since present programs have dealt 
with substances melting below 200 °C 

The usefulness of the equipment here described 
could be extended to the study of stability of pure 
substances under controlled conditions of, say, 
temperature or radiation. The purified sample while 
sealed in the tube could be subjected to the relevant 
experimental conditions and the decomposition 
products could be concentrated and isolated in a 
subsequent recrystallization cycle. 


The authors thank C. P. Saylor for his initiation of 
and interest in this project. The type of pump used 
(HS fig. 2) was the idea of G. Ross. The benzoic 
acid starting material was prepared by D. Enagonio. 
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This paper describes the effects obtained in phase contrast imagery with partially coherent illumination 
of phase objects. The objects chosen for study include simple Zernike type of periodic phase gratings of the 
form (i) A(u 1+7-A-cos wu and (ii) A(w)=1+7- > A,:- Cos nwu, With axial or oblique partially coherent 

n=l 
illumination. The concept of “effective source’? has been utilized, and the nonlinearity in the intensity 
transmission due to partial coherence has been analyzed in terms of “fundamental contrast” and “remnant.” 


1. Introduction 


It is known that ordinary methods of observation are not suitable for obtaining information 
about phase objects. Special techniques have to be used to render phase variations visible, 
namely, the dark-ground method of observation in which the direct beam illuminating the 
object is stopped. There is the other method, the Schlieren method, which differs from the 
above in that not only the direct light but also a part of the diffracted beam is also stopped 
from entering the ultimate detector. However, the most powerful method in this respect is 
due to Zernike [1934, 1942] which makes use of a phase modulation of the central or zero order 


of the diffraction produced by the object. In fact, it was shown by Zernike that with this 
method, which he termed as ‘phase-contrast’? method, phase variations in the object are 


practically linearly transformed into changes in intensity and thus made observable. Later, 
Hopkins [1952] pointed out that because of the finite size of the aperture, some details of phase 
structures could be seen in ordinary microscopes as well. 

The Zernike theory of the phase contrast method basically makes use of the Abbé theory 
of microscopic imagery in coherent light. We shall then refer to this as the Abbé-Zernike 
theory. A simple and obvious extension may be made to the theory to include the effect of 
amplitude modulation of the diffraction spectra, thus including the dark-ground and Schlieren 
methods of observation as well. In the original elementary theory of Zernike there were 
certain implicit assumptions which were later treated by Picht [1936] and Kahn [1955]. Oster- 
berg [1944] has also made a thorough study of the phase-contrast method as a problem in 
diffraction, in which he showed that the phase-plate characteristics could simply be treated 
asa coating function” modifying the pupil function of the microscope objective. His basic 
diffraction integral, which he terms as “generalized transport integral’ includes the effect of 
not fully coherent illumination, though in the final analysis he restricts himself to a narrow 
cone, or effectively coherent illumination of the object. 

In the meantime, the concept of partial coherence has been utilized by Hopkins [1953] 
in formulating what he describes as the generalized Abbé theory, the generalization being due 
to his inclusion of the effect of partially coherent illumination of the object. Steel [195s] 
also gave a generalized treatment of image formation in partially coherent light. We followed 
up the obvious next step in generalizing the Abbé-Zernike theory of “phase and amplitude 
contrast microscopy in partially coherent light.’ As a first step in this direction De and 
Som [1963] extended the study of Hopkins [1952] to show that phase structures could be ob- 
served in ordinary microscopes even in partially coherent light, the extent of observation 
depending on the magnitude of phase variations and the degree of partial coherence. 
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In this paper, we now discuss the effect of phase and amplitude contrast technique on 
the images of simple Zernike-type objects represented by A(w)=1+7 A(w), in which the mag- 
nitude of A is small enough to justify the presumption that the object shows no variable 
intensity transmission. In our treatment, we have utilized the theorem of effective source 
as postulated by Hopkins [1953, 1957], and the observation of Osterberg that the phase and 
amplitude contrast plate characteristics may be simply associated with the pupil function of 
the objective. This treatment may be further extended to cover a general type of phase 
erating represented by A(u)=exp |? A(u)], of which the Zernike-type object is just a special 
ease. We may further assert that our method is perfectly amenable to the consideration of 
even nonperiodic objects in a rather straightforward manner, and also of the effects of small 
residual aberrations in the objectives. Slansky [1962] also studied a similar problem, though 
his treatment is restricted to the case of Zernike-type objects only. 


2. General Solution 


In deriving the general solution of the problems of diffraction in phase and amplitude 
contrast microscopy in partially coherent light, we shall make use of the formulation by Hop- 
kins [1953]. In his so-called generalized Abbé theory in partially coherent light, Hopkins 
showed that, in all practical cases, the illumination of the object may be considered as due to 
a self-luminous source placed at the condenser exit pupil. This is termed the “effective source,” 
and insofar as the degree of partial coherence on the object plane is concerned, there seems to 
be no difference of practical significance between the two types of illuminating systems 
viz. the Kohler and the critical illumination. This and the further consideration of the 
progress of mutual intensity between a pair of points through the object to the objective aper- 
ture shows that this ‘effective source’? may be orthoscopically projected on to the objective 
aperture through the object point. 

Earlier Osterberg and Wilkins [1949] showed that with objectives satisfying Abbé’s sine 
condition, the diffraction plate characteristics may be simply treated as modifying the pupil 
function of the objective. Thus, the overall pupil function of the phase- and amplitude-con- 
trast objective may be taken as 


tla, yy=f (x, y) pla, y) 


where p(z, y) is the complex amplitude transmission of the diffraction plate. That this is pos- 
sible is seen from the fact that (x, y) denotes not only the coordinates of the point on the pupil but 
also the angles a ray makes with the axis as it comes from the object point on the axis to the 
point on the pupil. Thus, this also refers to a corresponding point on the back focal plane on 
which the diffraction plate is situated. 


Based on these two premises, the final intensity distribution in the image plane is given by 


B’(u’, 2’) ¥(20, Yo) O( Lo, Yo: U’, 0") *drodyo (1) 


ee 
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, a , +s 
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and where 
reduced coordinates of a point on the image plane, 


coordinates of a point on the image plane, 
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Lo, Yo fractional coordinates of a point on the pupil, 
sin a’—semiangular aperture of the objective referred to the object point, 
ay, by) —coordinates of a point on the pupil or on the effective source, 
radius of the pupil, 
refractive index of the medium between the object and the objective, 
is the effective source function described in terms of the pupil coordinates, 
partial intensity at the point (w’, v’) of the image plane due to an element 
dryly of the effective source. 


The physical significance of the above integral is obvious. To obtain the partial intensity, 
@°, we may proceed as follows: 

Let A(u, v)=complex amplitude transmission of the object, which may be Fourier analyzed 
into 


eet 
a(r, ¥) | A(u, v)e~**™ dudv 


giving the components of the spectrum of A(u, v). 


If we consider the object to be illuminated by an element drydy of the effective source at 
Zo, Yo, the spectrum would be shifted by a corresponding amount. Thus, the amplitude 
distribution of the spectrum would be given by a(z—z», y—Yyo), Which, after passing through 
the objective and the diffraction plate, will be modified to a(r—sp, Y— Yo) fla, y). The partial 
intensity, then, is given by the squared modulus of the inverse transform of the above. In 
other words, 


(a0, Yo: U’, 0’) | a(r—2o, ¥ yo) t te apye tO daly 


gives the complex amplitude distribution in the image plane due to an element of the source 
at ro, Yo It may also be interpreted as the complex amplitude distribution in the image plane 
when the object is illuminated by a plane wave front proceeding in the direction making an 
angle x9, yo with the axis. 

It now remains to evaluate explicitly the function ¢ pertaining to a specific object, which 
is possible if the composite pupil function is known. 


3. Solutions of Cases of Periodic Grating Objects of the Zernike Type 


We are now in a position to obtain solutions of specific cases. First, we shall consider 
those types of diffraction (or Zernike) plates which have a phase and amplitude contrast 
coating over an area which has the same relative size as the effective source. The complex 
amplitude transmission of this area, which we shall call the “conjugate area,” is, say, given 
by he'¥, while that of the remaining area, the so-called “complementary area,” is unity. / in 
the above implies amplitude contrast and y, phase contrast. 

Let us first consider the simple Zernike-type object: 


A(u)=1+7- A(u) (4) 


in which A is so small that A? is neglected. In other words, the intensity transmission of the 
above object is constant while the phase varies as A(w). Assuming A(w) to be periodic, we 
have 
A(u)=>5 A, : exp inwit| 
where 
| 


PAT] 7 


=O, 


A(v) exp [inwn|du 


2uy is the length of the fundamental period. @=(a/up). 
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The Fourier-spectrum of the object function (4) is thus given by 
a(x, y)=46(x, y) +7 >>) A, - 6(x+-nw, y) 


where 4(z, y) is the Dirac delta function. 
Thus, we have for ¢ as defined in eq (3) 


(29, Yo: Uv’) E I—2X0, Y—Yo) +t Dd) An-5(x—Ip+nw, y w) [ie y) exp {i(u’a+v’y) }dardy 


‘eae. f( Lg—N, Yo) exp —inwu’ exp i( 


o|)?=!| f(x, Yo) |? +2-Im { f(z, Yo) D5 A* f*(to—nw, Yo) exp (inwu’) (7) 


in which the asterisk denotes the complex corijugate, and Im | | indicates that the 
imaginary part is to be taken. Terms involving A’ have been neglected. The final intensity 
would then be obtained by an integration of \¢? over the effective source y(x, yo) as indicated 
in eq 1) above. The integration involved is of the form 


('(0, ») ¥ (205 Yo)F (20, Yo) F*(ao—nw, Yo) drydyo 


which is the generalized partially coherent response as described by Hopkins. 

Since we shall be concerned with pupils or Zernike plates which are symmetrical about 
the y-axis, f(2)+ ue, yo) =f(to—Ue, Yo) and if, further, we restrict ourselves to cases where A, 
is real, then eq (7) reduces to the following simple form: 


f(x, Yo) |? +2 >3 Im {f(x0, Yo) f*(to—nw. Yo) (2A,) COS nwu’. 
u=1 


The final intensity is thus given by 


NV I Y 
m; €'(0, n) 
B’(u’) LoS" (2A.) : 


— 


: cos neu’ (9) 
ea | C0, 0) 


in which a normalization is effected by the constant term C(0, 0). The C’s in the above equation 
are as defined in eq (8). The summation has been shown to extend from 1 to N, where 
NwS 1+ ) since C(O, n) shall be 0 for larger values of n. 

The problem now reduces simply to an evaluation of the C’s with given Zernike plate 
and effective source characteristics. 


4. Evaluation of C(O, n) 


In this section we shall describe the mode of evaluation of C(O, 7). We shall consider 
first axial illumination of the object with a circular disk-shaped effective source. The radius 
of this source projected on the aperture of the objective is, say, p. The actual radius of the 
condenser aperture may be obtained from the following [viz: Som, 1963}: 


nh, -sin Qe No? Sin ay No p+ sin a 


where ve and ny are the refractive indices of the media between the condenser and the object, 
and between the object and the objective respectively, ac is the semiangular aperture of the 
condenser from the object point, a is the corresponding semiangular aperture of the effective 
source projected on the objective, and @ is the semiangular aperture of the objective. In what 
follows we shall consider only cases in which p< 1. 
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Second we shall consider annular illumination, i.e. 


, cases of an annular condenser aperture 
with radii 1 and p,<1. 


The objective shall be assumed to be free from aberrations. Further, 
the conjugate area of the Zernike plate shall have the same size and shape as that of the “ef- 
fective source.” 


i) Axial illumination: 


Let us recall the expression for C(O, n) in eq (8) 


C'(0, n) | | ¥(2Xo, Yo) (20, Yo) f*(ap—ne, Yo) dx pd Yo. 


As indicated above, we may write 

y(2o, Yo) =1 within a circle 
0 outside. 

I (xo, Yo) =f (Lo, Yo) P(LosYo) 

A(X, Yo) =1 within a circle 
0 outside 

P(to, Yo) =he™ within a circle 
1 outside. (10) 


Figure 1 shows the geometrical implication of the integral (S). Evidently the region of inte- 
gration is complicated by the various limits imposed on the functions that appear in the inte- 
grand. The values of C(0, 7) have been obtained by subdividing the area of integration into 
different areas of overlap of two circles of same or different radu. These latter areas could be 
easily determined analytically, and with the values of the function within these different 
subdivided areas, one would obtain the final values of C(0, n) in terms of h, WY, n, &, p. 


They are 
eiven as follows: 


l me) } h sin y-mp" for w<l —p 
p 


( 
(“*)-5 i. ( 22 Tee. , D,(2 -nw- 7) } h sin y+ rp? 


2p 


for 1—p<w<yl—p 


{ 
f, 
{ 


fk 
D, 
’ Dy (== )+5 Di(=)+55 Dy (2 -nw4 7) } h sin p- mp’, 


where 


fornw>1—p’, 


l M ° _& 
D(x) 2. cos —sin - ( 2 cos ) . (12) 
T 2 2 


It should be noted that Do(z) gives the response of a perfect lens system imaging a self- 
luminous object. These values are already tabulated, and hence the formulas for C(0, n) given 
as above permit a straightforward evaluation, and also may be analyzed for aberrational effects 
with relative convenience. 
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Figure 1. Evaluation of C(O,n); axial illumination. Figure 2. Evaluation of C(O,n); annular illumination. 


(ii) Annular illumination: In this case, we have 


(20, Yo) =1 within an annulus bounded by circles 25+ ys=1 and 23+ y= pi 
0 outside 

f(xo, Yo) =1 within a circle 23+y5=1 
0 outside 


P(o, Yo) =he™ within an annulus bounded by circles 29+ yo=1 and 29+ y= 2 


1 outside. 


Figure 2 shows a typical geometrical representation of the area of integration in this case. The 
values of C(0, n) are as follows: 


Im {C(0,n) roi| 1 p,(“) |r sin y fornw<1—p,y 
Pa 
lw l l—pi—7w ] l 2+n*w ° 
roi| 1 D, (“* )—5 Do( at )+53 D,(—*&— )| hsiny 


for 1—p,<nw<y1l—p;3 


lw l l—pi—n-w | l—pit+7w ; 
iat) hf bbe Sh : - Do |: sin 
Tp | ( rs ar ( naan )+553 ( ) d v 
fornw> yl—p;.- 

C(0, 0) =2r(1—p?2)h . (13) 


5. Some Results 
(i) Axial illumination: 


a) Singly periodic grating: 


We shall first consider the results for a singly periodic phase-grating. 


Such a grating is 
represented by 


A i) +-9+A+COS wil 





Hence, from (9), we have 
B’(u’) +2a,-A-cos wu’ (14) 

where 
Im [C(0, 1)]. 


(15) 
(0.0) ’ 


With coherent illumination, ie., in the limit of the effective source y(z», yo) tends towards 
the Dirac delta function at thé origin, a, tends towards (sin y/h) for all values of w<1, so that 
the intensity distribution is given by 


; sin 
BAu’) rs: ene A (16) 
h 
Then the contrast in the image with partially coherent phase and amplitude contrast relative 
ayh 
sin y 


narrow-cone axial illumination, the contrast will be zero for w>1, while in the partially coher- 


to that with a limitingly narrow-cone axial illumination is given by a Further, with 


ent case, this limit is w>1+pwhena becomes zero. Thus, there is at least a possibility that 


there would be a gain in resolving power of the phase-contrast objective with a finite width of 
the illuminating beam. However, this is obtainable at the cost of contrast. 

The values of a have been calculated for various values of » and p. Figure 3 shows the 
results graphically. The following summarizes some of the pertinent observations in respect 
of the above. 

(1) So long as 2p <w<1—p, the imagery in partially coherent light is indistinguishable from 
that with narrow-cone illumination in so far as the contrast is concerned. Hence, if one can 
have a prior knowledge of the range of frequencies in the object, it might be possible to choose 


FiGurE 3. Contrast in the images of phase gratings; axial illumination. 
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a finite size of p to obtain a better photometric efficiency of the setup. Figure 4a shows the 
formation of the spectra under these conditions. All the spectra are wholly within the pupil 
and completely separated. That is why the ultimate image shows no difference in contrast 
from what is obtainable with coherent illumination. 

(2) If w<2p<1—p, the spectra are still wholly within the pupil, but they mutually overlap 
to a certain extent. Thus, not only the entire zero-order spectrum, but also parts of the 
first-order spectra are affected by the phase plate. This latter neutralizes the phase modula- 
tion, and hence contributes to a general loss in contrast. Figure 4b represents the state of 
affairs under these conditions. 

(3) If 1+p>w>1—p, parts of the first-order spectra fall outside the pupil and hence do 
not contribute to the image formation. The result is a further decrease in contrast. Figure 4¢ 
represents these conditions. 


ST a iT ai 
a a 8 














Fict RE 4a. Fo mation of the spectra: 2pwoX 1 p. Fiat RE 4b. Formation of the spect a; w 2p= l—p. 














Fiagure 4c. Formation of the spectra; | 








FIGURE 5 Useful frequency-range for partially 
coherent phase-contrast rmages 








The values of a as shown in figure 3 may now be analyzed in the light of the above observa- 
tions. Figure 5 shows (i) the range of frequencies which shall have a contrast same as in the 
coherent case for different values of p, and (11) the range of frequencies for which the relative 
contrast is equal to 0.8. This gives one an idea of the permissible maximum size of source 
that can be used for a satisfactory phase contrast image. 

The curves in figures 3 and 5 have been drawn for cases when h=1, that is, when no 
amplitude modulation is used. Further sin y has been taken as unity as well. This implies 
that the phase of the zero-order spectrum is advanced or retarded by 2/2. Noting that a Is 
proportional to sin y, the maximum phase contrast efficiency is obtained at this value. How- 
ever, d; is seen to be inversely proportional to h, meaning thereby that the contrast may be 
enhanced by using a slight amplitude contrast as well. This, of course, is possible at the cost 
of the mean level of illumination of the final image. 

(b) Square-topped phase-grating: 

This is the case of a simple periodic structure where troughs and crests represent differences 
in optical path only. 

Let 

A() +7 A(w) 


where A(w) is periodic in uw with a fundamental length of period 21. , over one 
fundamental interval, 


A(t) Q) aAUyS ou <u, 


A ul <alo. 
The Fourier spectrum is given by 


A A ae 4. 18 


his 


ain the above gives the ratio between the line width to line separations. Hence, from eq (9) 
we have 


ww C'(0, 0) 


NV j { ) 
. / SIN:a:ir Im ((O.n . 
B'(u ) l 1DA->\2- | COS Tiwi 
1 


where Nw < 
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In the limit of a point source, the expression for intensity will be 


N - \ /as 

‘ : SIN + anr sin . 

Bu’) =1 +2 >) (2. =\ YY cos « new’ (20) 
2=1 nw h 


where Nw <1 since Im [C(0, n)] in the coherent limit will be / sin y for all values of nw<1. 

Figures 6a to d show the image of such gratings (A=0-5) for four different frequencies, and 
for various effective source sizes at each frequency. The dotted curves represent an ideal 
ease of linear unitary transformation of phase variation into amplitude variation. 

The expression for B’(u’) in eq (19) or for Bi(u’) in eq (20) shows that the final image is 
given by a truncated Fourier series, and hence one would expect to observe the so-called Gibbs’ 
phenomenon of large oscillations in the region of the edges of the square-topped grating. In 
the presence of harmonics in the image intensity distributions it is indeed a bit difficult to 
define a unique contrast for such cases. The presence of these harmonics is in fact an indica- 
tion of the effect of nonlinearity in the intensity transmission in partially coherent light, and 
as vet no unique specification has been suggested for a proper measure of such nonlinear effects. 
One suggestion is to express the contrast in terms of the fundamental frequency distribution, 
and lump together the effects of higher harmonics in a so-called “remnant” [Ingelstam, 1964]. 
This “remnant”? may be described as the square root of the sum of the squares of all coefficients 
other than the constant and that of the fundamental frequency component. Figure 7 gives a 
plot of this fundamental contrast and the remnant against source sizes for different frequencies. 

It is of interest to note that the remnant always decreases with increasing frequency, 
falling gradually down to zero. In fact, if N<2, there are no harmonics in the image intensity 
distribution. Thus, for a frequency up to 0.5, harmonic distortion would be present for all 
source sizes, though its magnitude as reflected in the value of the “remnant” may be negligibly 
small. For values of w lying between 0.5 and 1, the “remnant” would be zero only if p is less 
than a certain value given by p<.2w—1. 








I maae s of square -lop pe 


d phase gratings FIGURE 6b Images of sq sare-topped phase gratings; 
arial illumination; @=—0.25. 


arial illumination; w= 0.50. 








Figure 6e. Images of square-topped phase gratings; } Figure 6d, Images of square-topped phase gratings; 
arial illumination: w= 1.0. axial illumination; w= 1.5. 


Figure 7. ‘Fundamental contrast’? and 
“remnants” in the images of square-topped 
phase gratings; axial illumination. 








In the cases for which computation has been done, a has been chosen as equal to 1/2, im- 


plying equal widths for the troughs and crests of the phase grating. Jn these cases, the spectra 


») 


consist of only the odd harmonic components and hence the “remnant” would be zero if N<3; 
i.e., if p<<3w—1, w >0.33. This expresses the range of validity of using the criterion of ‘‘funda- 
mental contrast” for the study of the efficiency of phase-contrast devices in partially coherent 
light. 

Gi) Annular illumination: 

(a) Singly periodic grating. 

The expression for intensity is the same as in the previous case, except that the values of 
Im [C(O0, n)] to be used are as given in (13). 
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If the width of the annular source tends towards an infinitesimally small value, the area 
integral for C(O, n) will tend towards a line integral, and we shall have 


° , . N@ . 
lim Im [C(0, n)| 2-cos!- ( = yh sin yp 
py?! of 


lim [C(0, 0)|=27h?. (21) 
Pa! 


Thus, unlike the case of axial illumination discussed above, we have in this case a contrast 
varying with frequency even with a source of infinitesimally narrow width; we have, in this 
limiting case, 


nw 
COS 


l _ * 

“ sin P 

+DA.- ( ¥) cos wu’. (29) 
T h , 


The difference with the corresponding expression (16) for the axial case is to be noted with 
interest. 

Evidently in this case the limiting resolvable frequency is nw <2. This is just a theoretical 
limit, as the contrast also tends towards zero at the limit. 

If we now consider the case of the annular source with a finite width, we shall find that 
Im [C(0,)] becomes zero for nw >1+ p, which implies that we lose on the count of ultimate 
resolving power. Figure 8 shows graphically the values of a,=Im [C(0, n)]/C(0, 0) as a funetion 
of frequeney » for various values of p,. The dashed curve in figure 8 indicates the range of 
frequencies which are so imaged that the relative contrast is nowhere less than 0.8. 

(b) Square-topped erating: 

As in the case of axial illumination, we have, in the limit of a line source coinciding with 
the margin of the pupil, 


mee et) p 
COS + Th@il 


Wa T h 


where 


> and N<(1-+ p,)/w. 


annular illumination. 











Figure 8. Contrast in the images of phase gratings; 


With a finite width of the source, we have an expression similar to (20) above. Figures 9a to d 


; : : , => ; T , F 
show the images of such gratings (A=0.5) of a number of frequencies w——~’ where 2% is the 


length of the period, and a=} implies equal widths for the troughs and crests of path-differences 


FIGURE 9a I mages of square-topped phase gratings; FIGURE Yb. Images of square-top ped phase gratings; 
0.25. | annular illumination; w=0.50. 


annular illumination; w 





Figure 9c. Images of square-topped phase gratings; Figure 9d. Images of square-topped phase grat 


annular illumination; w= 1.0. annular illumination; w= 1.5. 





Figure 10. “Fundamental contrast’ and ‘“‘remnants”’ 
in the images of square-topped phase gratings; 
annular illumination. 











in the grating. The dotted curves represent an ideal case of linear unitary transformation of 


phase variation into amplitude variation. Here also we may make the same remarks as done 


in the case of axial illumination of square-topped gratings. Figure 10 describes the situation, 
and the conclusions are also similar. Remembering that the coherent limit of resolution in 
this case is 2, we find that practically the upper half of the range of resolvable frequencies is 
still assessable in terms of the “fundamental contrast” only. 


This study has been undertaken under the sponsorship of the National Bureau of Stand- 
ards. The authors are deeply indebted to the stimulating discussion with 5S. C. Som on the 
various aspects of this study. 
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Exposure Time Relations for Kossel Microdiffraction 
Photographs 


H. Yakowitz and D. L. Vieth 


(April 23, 1965) 


The exposure time for a Kossel photograph may vary from a few seconds to a few hours. 
Therefore, it is desirable to be able to estimate the exposure time for various experimental 


conditions. 
pattern have been developed. 
reflection Kossel photographs. 
commercially available x-ray films. 


Hence, semiempirical relations for the exposure time of a Kossel microdiffraction 
Equations are presented for both transmission and back 
These equations are tested for validity using two different 
It is shown that the agreement of actual exposure times 


with predicted exposure times is valid within 10 to 15 percent. 


It has long been known that the exposure time of a 
Kossel photograph has a well-defined optimum [1].! 
However, useful analytical expressions enabling one 
to calculate the exposure time have not been pre- 
sented. Therefore, the purpose of this paper is to 
propose expressions for the exposure time in both the 
transmission and the back reflection Kossel regions. 

In the transmission mode, the Kossel conics 
usually appear light on a darker background and the 
contrast is, at best, not good [2, 3]. In the back 
reflection region, the conics appear darker than the 
background, and the contrast may be somewhat 
better than in transmission [2]. The only practical 
source for Kossel patterns is a finely focused electron 
beam which is allowed to strike either the sample or 
a source foil in close proximity to the sample. The 
latter case will be discussed as it is the more general 
and more useful. It will be assumed that the film is 
in a vacuum and that once the x rays leave the 
sample, having undergone the usual exponential 
reduction in intensity, they travel unimpeded to the 
film. If an x-ray window and airpath intervene, the 
reduction in intensity of the x rays emergent from 
the sample can also be accounted for by the usual 
exponential retardation law which simply appears as 
a multiplier. 

Using a focused beam of electrons, one may vary 
the specimen current, i.e., the number of electrons 
flowing to ground per unit time from specimen or 
foil, and the accelerating potential of the electrons. 
A knowledge of the number of photons per electron 
which strike the film and of the film area is necessary. 
Some provision for the fact that the exposure is not 
constant over the expanse of a flat film is also 
necessary. Bearing each of these requirements in 
mind, we may write for the exposure: 


—_$______— 


1 Figures in brackets indicate the literature references at the end of this paper, 





| p (2,) erg.) 
1, ( ( 
fi (000) G5 G) ( { ) 


rey 


(eq) 7 > —_ = 
6.2510" [exp [—(ayr,+hss) 
E ] t | t M oo Le 


) is the exposure in photons/cm?. 
is the specimen current in microamperes. 
is a factor accounting for variation of 
exposure With position on the flat film. 
, is the solid angle subtended by the spheri- 
cal segment containing the film. 
A, is the exposed film area in em’. 

(nx) ,/4x is the total number of AK photons pro- 
duced per unit solid angle by the foil per 
incident electron. 

6.2510" is the number of electrons 
coulomb of charge. 
is the electron backscatter coefficient of 
the foil [4]. 

wyand g, are the equivalent linear absorption co- 
efficients for the total spectrum produced 
in the foil and the sample respectively in 
em™ {1}. 

z,and x, are the thicknesses in cm of the foil and 
sample respectively. 

t is the exposure time in minutes. 

Values of no=[(nx),/4x] as a function of atomie 
number and the accelerating potential of the elec- 
trons are available in the literature [5, 6, 7). The 
area of the film can be found by reference to figure 
1 as: 


where 


per micro- 


Ap=rl> tan? n, (2) 
where 


D is the source to film distance in em; 
n is the half-cone angle subtended by the film. 
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The values of D and 7 are usually well known. 

The inclusion of @ arises from the fact that the 
exposure will be greatest at the film center and least 
at the edge. This is a combination of two effects, 
the first of which is an increasing absorption path 
within the sample with increasing angle of emergence. 
The other effect is that of the geometry, 1.e., a unit 
of solid angle subtends a larger area on the film at 
the edge than at the center. Jt has been observed 
that the net effect is a slowly decreasing function of 
the angle 7». This can be suitably approximated by 


defining G as: 
G=|! — "| 
2 


Over the usual range of n values, i.e., 20°<7<36°, 
the inclusion of this form of G increases the exposure 
time 3 to 10 percent. It is interesting to note that 


the important central film area in figure 1 defined 
by the angle /2 and the ray, Q, will be nearly 
uniformly blackened. 


The factor — is 
|—? 
the beam current, i.e., the incident electron current. 
It is for this current that np values have been derived. 
The simple value “r’’ can be used since the effect on 
the specimen current, 7,, due to secondary electrons 
which may be generated is quite small. 


RAY SOURCE 


SPECIMEN ——™ 


ACTUAL FILM 


AREA BLACKENED 


FIGURE 1. Geometrical relations between the x-ray source and 


the filn ° 


The solid angle subtended by the spherical segment 
containing the film is given by [S]: 

Q,=27r(1—cos n). (4) 

Thus, (m2) is a measure of the number of photons 
per incident electron which would actually fall on the 
film if there were no exponential absorption. 

The value @ can be approximated by the relation 
uu, w being the linear absorption coefficient for 
Kk, or L,, for excitation over-voltage ratios of 3 
to 4 [3] and for most uw values used in Kossel analysis. 
The value of x, should be kept as small as practical. 
The value of r, should be the approximate optimum 
thickness for the transmission case [3]. In back 
reflection, an effective depth can be taken as an 
approximation to a,. This point will be discussed 
later. 


Equation (1) may now be rewritten as: 


(7,) Leos n (] 
(t) 2r 
(1—r) 2 


COS 7) iio 
rl)” tan® 7 


(3.75 X10") {exp—[32 (uyrp+unr.)]}. (5) 


Solving for ¢ we obtain: 


2.67 X107E (sec? n)D? 


lhots 


t (min) Sng) 


exp—[2 (uyr,+usr J}. (6) 

It is necessary to obtain a value for F empirically 
for each type of film to be used. This F value, 
henceforth called ££, is that exposure density 
yielding the maximum contrast between the Kossel 
conics and background for a given film at a given 
distance, D. Jt is emphasized that the /) value in 
eq (6) represents a compromise between the exposure 
at the extreme edges of the film and the exposure at 
the film center. It is to be expected that the ex- 
posure density / will be essentially constant for a 
given film type, independent of other camera 
parameters. 

Equation (6) as written is strictly speaking only 
applicable to transmission Kossel photographs. In 
the back reflection region, precisely the same re- 
lation holds true except that x, is undefined. How- 
ever, a reasonably good approximation to sz, can be 
made. 

Such an approximation can be made following 
Il’in’s treatment [9]. Il’in proposes an _ effective 
thickness, d,, of the absorbing layer chosen such that 
the function exp—(ud,) would give total attenuation 
of an analytic line such as Ka, or La;, in the working 
volume of the specimen, i.e., at 7»=90°. It is pre- 
sumed that 


pd, H (7) 


in which p is the density of the specimen and // is a 
constant for a given 7 value. As support for eq (7) 
the fact that d, is related to the total thickness of 
the emitting portion of the target is invoked [9]. 
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According to the Whiddington law, this thickness 
is determined for a given accelerating potential by 
the density of the emitter. 

By rearranging II’in’s relations [10] and invoking 
electron probe microanalyzer results presented in the 
literature [11-14], it is possible to plot an /7 versus 
n curve. This plot is shown in figure 2. 

Noting that zs, is to be replaced by d, in eq (6) 
for back reflection cases, one obtains: 

Mgt s 


uss H(u P)s- 


Hence for back reflection, eq (6) becomes 


2.67 


10—-5F (sec? n) D? 


(Vig) (2) 


t (min) (1—r) 


exp —2[u,2, +H (u/p)s]} (9) 


In order to test the validity of eqs (6) and (9), 
the value of / was calculated for two different film 
types commonly used in the preparation of Kossel 
microdiffraction photographs. One of these was 
moderately coarse-grained duplitized film ? while the 
other was a moderately fine-grained single emulsion 
film.’ These shall be designated I and II 
respectively. 

Using type I film for about 200 Fe-K radiation 
exposures of Fe-3 w/o Si alloy, it was found that the 
value of Ay was 1.6 10° photons/em? for high con- 
trast, high resolution transmission photographs [15]. 

For type II film, the exposure time versus iron 
thickness curve given by Morris and Ogilvie [16| 
Was used to calculate /),. The value obtained was 
1.410° photons/em*. Inserting the values of F, 
and /),; into eq (6) we obtain: 

4.1 107'D? (sec? n) 


) (No) 


t; (min) (1- 


r) 


(7 


exp —[3 (uyr p+ ya's) (10a) 
4.0 10-*D? (sec? n) 


t;, (min) . 
it (Ve) (No) 


exp (10b) 

It remained to check eqs (10) with independent 
data. For type | film, Gielen’s exposure of germa- 
nium with Cu-K radiation was used [17]. The re- 
quired parameters were: 


6 ¢m 

20.3 

O.lluA 

2.2107 photons/electron/unit solid angle at 
30 keV [7] 

1.8 X107* cm 

55 em7! 

.o2 X10 

1 


l 
4 
l 


- em 
385 em 
0.29. 


Vou 


Kodak AA fil 
Kodak M filn 


21 


- 


S] 





H x 10°, g/cm? 








) , DEGREES 


») 


FIGt RE <2. The value of H as a function of the half-angle of 


the film, n. 


Gielen used a camera with an aluminum foil light- 
tight cap. The (ur) value of the foil was 0.20. 
The camera operated in vacuum. Using these data, 
eq (12a) predicted an exposure time of 49 min. The 
actual exposure time for a high-contrast photograph 
was 50 min [18]. 

In order to check both eqs (9) and (10b), a back 
reflection photograph of iron using Fe-K radiation 
taken on type II film was used as a model [19]. The 
required parameters were 


12.4 em 

26° 

2.0uA 

6107? photons/electron 
at 40 keV [6] 


unit solid angle 


0:27. 


Using these data, the predicted exposure time is 
0.47 min or 28 sec. The actual time was 25 sec {19}. 





It is perhaps of interest to comment that for a given 
set of conditions a simple equation results. For 
example, using the equipment parameters given for 
the iron exposure in back reflection, one obtains: 


78><10"* ero Ze os 
ty -—— (1—r) exp [5.7 107*(u/p),+0.67y,2,] 
(Mn) (de) oA ee 
(lla) 
14xXi90-* 


(No) (4s) 


ty (l—r) exp [5.7107*(u/p),+0.67y,2,). 


(11b) 


It has been shown that eqs (6) and (9) represent 
the correct form of exposure time relations for Kossel 
photography. Because of the approximations made, 
e.g., for « and d,, the uncertainties in absorption 
coefficients themselves, and possible inherent sample 
limitations [1], the value of an empirical £) is 
probably good to only 10 or 15 percent. Hence, 
agreement with eqs (10) and (11) to only about 
10 to 15 percent is to be expected. Nevertheless, 
this is a significant improvement over a pure trial 
and error method. Furthermore, as more exposures 
are taken /, values can be refined. 

In the transmission method, there appears to be 
a definite upper limit on (u,7,) in order to obtain 
any pattern [3]. This limit is about p,r,=10. It is 
therefore recommended that eq (6) not be employed 
with values of u,r,>10. 


used to 


The photographs 
values were prepared while one of the authors 
(HY) was a guest of the Massachusetts Institute 
Technology assigned to the laboratory of Prof. 
R. E. Ogilvie. 


determine the £, 
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Cartesian Diver as a Density Comparator 


Horace A. Bowman and Randall M. Schoonover 
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\ hydrostatic weighing system is described that provides independent values of apparent 
mass with a standard deviation of about 0.2 microgram. The definition of “independent” 
used here includes the requirement that the sample under observation be removed from the 
water, dried, and recleaned between measurements. This precision is between one and two 
orders better than existing hydrostatic balances, and permits high-quality density measure- 
ments without recourse to large sample sizes. The system is used as a comparator, hence 
the absolute accuracy of results can be no better than that of the standards used. Data 

re presented which were taken during experiments on two 2-gram pieces of single crystal 
silicon, and the standard deviation in density of a single determination was computed to be 
4.41077 g/em® (0.22 ppm). 

The system is particularly well suited to detection of slight density changes in small 
size samples. In this service, it is not subject to the requirement for high accuracy standards. 
Data are reported on density changes of about 50 ppm that occurred in a 250 milligram 
ceramic crystal. Two independent determinations on this change differed by only 7 percent. 


1. Introduction with the meniscus surrounding the point where the 
support wire penetrates the water surface. Under 
The density of an irregularly shaped object. is nonideal conditions, these forces may vary by 10 ug 
usually determined by a combination of air and | Or more. 
water weighing in which the apparent mass of the | .. Thus, the accuracy in a measured value of B is 
object is determined in both media. The apparent limited by our knowledge of py, and the eben: 1 
mass in air, A, is determined on a balance in the | bility of the measured value is limited by the 
usual manner. The apparent mass in water, B, is variability of the nonbuoyant forces mentioned 
determined by suspending the object from one pan | above. ; 7 
of the balance by a fine wire into a body of water. For a given hydrostatic weighing system the ab- 
The mass, 7, and volume, V (the ratio of which is | solute pgs of the variability (due to nonbuoyant 
density), result immediately from a simultaneous | forces) is largely independent of load. Hence, to 
solution of achieve 1 seit. values of percentage re produc ibility 
in B, we must use samples of very large sizes. 
A=M—p,V (air weighing) Unfortunately, objects for which density deter- 
minations have the maximum scientific value are 
B=M—pywV (water weighing) rarely available in samples of unlimited sizes. It is 
seldom that samples above 2 ¢ are obtainable, and 
where p, is the density of the air displaced by the | they are usually of fractional-gram size. Such 
object during the air weighing, and py the density | items as uniformly irradiated crystals, samples of 
of the water displaced by the object during the | maximum purity subst: ances , and balance sensitivity 
hydrostatic weighing. weights, usually weigh 250 mg or less. For a 250- 
The main limitation in the above scheme occurs | mg object, a variability of 10 ug in the nonbuoyant 
in the determination of B. There are two reasons | forces ool result in a variability in B of about 
for this: 40:10°—a highly unsatisfactory measurement for 
(1) Although we know the density of a particular | many technical purposes. 
sample of water to almost two orders better than In the Cartesian diver system to be described, 
we know the density of a particular sample of air, | is possible to determine B of a 250-mg object to a 
water density is about three orders greater than that | reproducibility (standard deviation) of about 1:10". 
of air. Hence, the uncertainty in pyV is greater | The reason for this is that the two major sources of 
than that of anaes variability in ordin: ary hydrosts atic weighing do not 
(2) In water weighing, the variability of non- | exist in the diver system or are minimized: 
buovant forces is of much greater magnitude than (A) No portion of the diver structure penetrates 
those occurring in air weighing. First, there are | a water Ah so uncertainties caused by surface 
variable momentum forces associated with turbu- | tension do not exist. 
lence and convective currents in water which act on (B) The diver system is operated at a nominal 
the object being weighed. Such forces also exist in | temperature of 4 °C; hence, nonideal thermostating 
air weighing but to a much smaller extent. Second, | does not result in variability of water density. In 
there is variability in surface tension forces associated | the absence of density gradients (other than those 
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caused by gravity) convection currents cannot be- 
come established, and turbulence damps out quickly. 
Thus momentum forces are minimized. It is gen- 
erally considered impractical to perform ordinary 
hydrostatic weighing at 4 °C. 


2. Description of the System 


The Cartesian diver is not a new device. A re- 
view of prior art is given by Chiarulli and Chilton 
f1}.!. Recent work by Haller and Calecamuggio [2] 
and Spaepen [3] has raised the state of the art 
greatly. 

Basically the diver, figure 1, is a hollow, compres- 
sible body with a closed surface to which are attached 
various hooks, guides, and reticles to facilitate load- 
ing, manipulation, and observation. When such a 
body is immersed in water and the pressure in the 
water is increased, the positive buoyancy of the 
body will decrease or increase respectively, depend- 
ing upon whether the compressibility of the diver is 
ereater than or less than that of water. If diver 
compressibility is equal to that of water, buoyancy 
is independent of pressure. 


1 Figures in brackets indicate the literature r 
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Our most useful diver is made of fused quartz 
(vitrified silica) and weighs about 5 g. It displaces 
about 10 em*; hence in water it has a positive buoy- 
ancy of about 5g. It has a very thin wall so it is 
quite compressible. In practice it Is loaded by the 
object to be measured and enough tare weight to 
cause it to sink slowly in water at 4 °C. The pres- 
sure system, figure 2, is then closed and the water 
pressure reduced slightly. This causes the diver 
volume to expand, thereby increasing diver buoy- 
ancy. This causes the loaded diver to stop falling 
and begin to rise. 

There is a pressure gradient in the water, due to 
the gravity field, so that as the diver rises, it moves 
into a region of lower pressure. This causes addi- 
tional expansion and an increase in buoyance, so the 
diver is accelerated upward. Conversely, when the 
water pressure is increased, the diver body is com- 
pressed which results in smaller buoyance. The 
diver will then sink. As it sinks into a region of 
higher pressure, additional buoyancy is lost and the 
diver is accelerated downward. Thus, the gravity 
pressure gradient in the fluid always accelerates 
diver motion, and as a result, there is no position 
of stable equilibrium. 

On the other hand, 


had the diver compressibility 
been less than that 


of water, the reverse set of 
circumstances would have held. Under such con- 
ditions, the pressure gradient in the water would 
always impede the diver motion and a position of 
stable equilibrium would exist. It is possible to 
operate a diver system under either condition, but 
practical considerations such as load range and 
sensitivity led to our choice of a diver more com- 
pressible than water. 

Diver motion is observed through a_ rigidly 
mounted microscope and a system of windows 
and mirrors. The optical center of this system is 
carefully adjusted horizontally, and the horizontal 
plane within the diver chamber containing the 
optical center is called the reference plane. Although 
(as explained above) there is no level of stable equi- 
librium for the diver in the system, there is a pressure 
which will result in the loaded diver having a density 
precisely equal to the density of the water it displaces 
when the observed upper tip of the diver is in the 
reference plane. This pressure is called the equi- 
librium pressure, ?,, and it can be determined by 
averaging the two pressures which cause the diver 
to move upward and downward across a vertical 
range defined by an eyepiece reticle and extending 
about 0.013 em above and below the reference plane 
in equal intervals of time (usually about 30 see). 

When the pressure is very close to Pe. we are 
justified in assuming that the average velocity, y, of 
the diver, during its transit of the range, is nearly 
linear with pressure. Thus, when the average 
rising velocity, ype, does not precisely equal the 
average falling velocity, ye, an acceptable value of 
Py, can be caleulated from the formula 
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Alternatively, there is a quick graphic solution, DIVER VELOCITY 
which we used, illustrated in figure 3. 0.4 L (IN EYEPIECE DIVISIONS 
PER SECOND) 


3. Pressure System DIVER MOVES UPWARD ACRO 
DIVISION RA 
OCITY OF 0.334 
Our pressure regulator system is a modification of DIVISIONS / sec. PRESSURE WAS 
the two-bellows barostat developed by D.P. Johnson i. 
{4] for use in the NBS barometer calibration service. CHECK POUT 
It consists of two stainless steel bellows mounted 
one above the other. The upper bellows is evacu- 
ated, and the lower is water-filled and connected 
to the diver chamber. The upper end of the upper 
bellows and the lower end of the lower bellows are 
rigidly attached to the frame. The two bellows are 
connected mechanically together by a freely floating 
separator plate. Under such circumstances, the Las Ps “er 
pressure developed in the water in the lower chamber 
is directly proportional to the deadweight loading Pe TAKEN TO 
of the separator plate. The effective areas of the BE 71.3370 kg 
bellows are such that a 76 kg deadweight load on the 
separator plate results in about 1l-atm pressure 
being developed in the lower chamber. The loading 
of the separator plate is conveniently changed in 
increments of about 100 g, from a minimum of 2 
ke (the tare load) to 85 kg, by a semiautomatic 
loading machine developed by David Sklar, Design 
Development Laboratories, for the purpose. A OBSERVED POINT 
small tray on top of the assembly is used to receive DIVER FELL AT 0.347 —______. 
additional small hand-loaded weights of between 1 pay fay ree 
and 100 ¢ (see fig. 2). - , 

The separator plate is electrically connected as 
the moving plate of a two-plate variable capacitor. 
The fixed plate is rigidly attached inside of the 
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evacuated bellows, about 0.013 em above the mid- | 


point of the 0.025 cm free travel allowed the separator 
plate by its stops. The capacitance between the 
tixed and moving plates is fed into a Wien bridge 
oscillator and discriminator, which provides servo 
power to drive a motor-driven bellows connected in 
the hydraulic system, in such a direction that the 
separator plate is maintained at the midpoint of 
its allowed travel, regardless of loading. 

With this arrangement the two-bellows assembly 
is used as an error signal generator and actual pres- 
sure changes are made by the servobellows. Inas- 
much as the separator plate remains at the same 
vertical level, regardless of pressure, the effects of 
the elastic properties of the stainless steel bellows 
are largely eliminated. 

The elastic properties of one of the divers are 
such that the diver buoyancy changes by 1 ug when 
the separator plate load is changed by 12 g._ Inas- 
much as diver buoyancy on differences is constant 
to about 2/10 yg, it can be inferred that the short 
term stability of drift of the pressure regulator is 
about 210° atm/hr. This is consistent with 
other estimates of stability of similar devices [5]. 

As use of the diver system does not require absolute 
knowledge of pressure, no effort has been made to 
calibrate the system in absolute terms. In_ use, 
pressures developed are stated in terms of the 
kilogram load on the separator plate, which is 
acceptably linear with water pressure. 


4. Use of the Diver—Density Measurements 


If G, is the mass of water displaced by the Carte- 
sian diver body at some pressure, then 


G; Pi?! 


where p, is the density of water at that pressure, and 
vr, the associated diver volume. Both p,; and 7, are 
functions of pressure, and G will increase or decrease 
with increasing pressure depending upon whether the 
compressibility of water is greater than or less than 
that of the diver body. AG is not a perfectly linear 
function of AP. We have carefully measured the 
nonlinear term and found 


A’G 


AP’ 8x10 


So/em Hg/em He.” 


Integrating twice gives us 
G 


where (’; and Cy are the constants of integration in- 
troduced thereby. This term adequately represents 
the mass of water displaced by the diver as a function 
of pressure. 


Although, as previously stated, the pressure system has not been calibrated 
in absolute terms, within 1 percent, a lo 


1 | unge of 1 kg on the separator plate 
results in a pressure change of 1 cm Hg 


The experiment consists of loading the diver, first 
with two standards, 51 and S82, and noting the equi- 
librium pressures, /s; and Pg» associated with each. 
Following this, the unknown object, Y,is placed on 
the diver and P, observed. The masses, Mg, and 
Mego, and the volumes, Vg; and Vgo, of the two stand- 
ards are assumed known. The equilibrium equation 
of the diver is 

M,—G T M,—pwV1 0 

where .V/, is the diver mass, ./, and V; the mass and 
volume, respectively, of the load and py the density 
of water at the equilibrium pressure. As pointed 
out in the introduction M,;—pwV_, is the apparent 
mass of the load in water, B,. For each of the three 
loads placed on the diver during the experiment, we 
may write an equilibrium equation 

(A) Myp—(4 107°) P3, 


-C’":Ps:—C2 Bs; 


(B) My—(4 


(C) Mn—(4> 

subtracting the second equilibrium equation from 
the first, and then the third from the first gives us 
two difference equations 

— (4X 1078) (P2,— P2.)—Ci (Ps: 


Pas) T (Bs, — Bz) 


—~ (4X 1078) (P2,— P%) —C\(Ps:— Px) +(Bs:— Bx) =0. 
The first difference equation may be immediately 
solved for a numerical value of C;, 

10-*( P§:— Phe) 

P32) 


and with this value of C, inserted in the second 
difference equation, we may solve for the apparent 
mass of the unknown in water, By. 
Bx=Bg.—4X107*(P§,— Px) —Ci(Pai— Px). 

This experiment demands a minimum of three 
observations—Ps;, Pg, and Py. Our system has 
drift associated with electronic drift, drift in ambient 
conditions and (probably) other sources. In order 
to minimize the effect of this drift we use one of our 
loads as a drift monitor. This load is observed 
every second observation. For an experiment with 
only one unknown the observation format would 
contain five observations taken at uniform intervals, 
thus: Ps:—Ps:—Ps:i—Px—Ps:. The value _ of 
(Ps,—Ps2), required in the calculations shown above 
is obtained by subtracting the observed value of Ps. 
from the average of the first two Pg)’s. (P?si:— Px) 
is obtained in a similar manner. The three ob- 
served values of Ps; are averaged to give the effective 
value of Ps, during the experiment. Effective 
values of Ps. and Py used in calculations of C, and 
By, are obtained by applying the difference terms 
to this effective value. 
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The values of water density used in the calculation 
of apparent mass in water of the standards (from 
assigned values of mass and volume) are based upon 
an assumed density of 0.999973 g/cm’ at 4°C and at 1 
atm, and corrected for departure from | atm of the effec- 
tive values of observed pressure. Errors in the assumed 
value of water density exert only a second-order 
effect on the calculated density of the unknown, 
because our density knowledge is primarily obtained 
from the density standards used in the experiments. 
The water is merely the comparison medium. 

The diver system is currently being used in several 
density studies not pertinent to this report. To 
illustrate typical diver operation, some data from 
these studies will be presented here. Table 1 shows 
the observed raw data taken during measurements 
on two imperfect silicon crystals, X1 and X2. The 
standards used in this experiment were two perfect 
silicon erystals, SI and S2, whose assigned density 
value was 2.3290040 g/em*. The observed values of 
P, and the pressure differences were calculated from 
values of rising and falling pressures and velocities 
as previously explained. This table shows data 
taken during four independent experiments per- 
formed during January 1965, and the average value 
of a group of nonindependent similar experiments 
performed on the same crystals during the previous 
November. Between each experiment all crystals 
were removed from the hydrostatic system and dried 
and recleaned. 


TABLE 


Rising 


Pressure Ve- 


locity 


66, SOLO 2 66.7310 0 
2264 v4 64. 1564 

}. SOLO 922 66.7310 

5. OSO09 7 65. 0109 

». S110 - [Sde 66. 7410 
1246 Vi2 64. 0546 

> SLIO ()2 66. 7410 


> 8210 79 66. 7510 
2459 92 4.1859 

. S210 202 . 7510 
5. 1009 . 1766 5. 0309 
), 830) 72 1H. TH00 
. 1464 4. O764 
. S210 155 . 7410 


». YOO . 2056 >. S700 
3851 4 4.3151 

». 9400 912 S700 
5. 2201 {2 5. 1501 
». YAOO 183 iH. STOO 
5 2159 

). S800 


}. 9092 


}. 9092 
3559 


). 8992 


ints taken in November 
Silicon 


Table 2 shows the calculated values of effective 
P, and the assumed py associated with these values. 


Theapparent masses in water of the standards is based 
upon assigned true mass values of J/s;=2.0000912 g 
and Mg.=2.0004884 ¢. The effective values of 
(’, were obtained as previously described and they 
led to the values of apparent masses in water of the 
unknown cyrstals. 


TABLE 2 


Crystal 31 


Effective Pez l 66. 7692 64. 0863 
66. 7868 64. 1097 
66. 9094 64. 2377 
66. 9420 64. 2832 
67. 73bt 65. 0508 
Assumed pw 0. 999678 0. 9999660 


0. 999661 yO99667 


660 661 667 
661 62 HOS 
662 663 668 
667 668 673 
3 of standards 1. 1415717 
1. 1415717 
5716 
3434 5716 
3430 5711 
—(), QOO0Y9O2S 
9047 
YOH. 
Y1L09 


3 of unknowns 


As pointed out in the introduction, we may 
obtain values of Wy and Vy from a simultaneous 
solution of the air weighing and the water weighing 
equations. Alternatively we may enter the water 
weighing equation with the value of By obtained 
with the diver and a previously determined Wx 
and solve for the density of the unknown, px. 


p Pu 


where py is the assumed water density at pressure 7’ y. 
Very careful analysis of this point will show that the 


two methods are nearly identical. The values of 
px Calculated from the By’s determined with the 
diver and \/y,—2.0004774 ¢ and MW y.=2.0003478 ¢ 
were: 


il X1 


8289950 
2. 3289958 
2. 3289944 
2. 3289950 
2. 3289956 


6 January 
7 January 
Ss January 
9 January 


Nov. Avg 


BS290012 
B290018 
S2900017 
BS2O00L5 
3290019 
S2840951 


Average 2290016 


Pooled standard deviation 


ppm. 
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5. Use of the Diver System—Density 
Changes 


There is widespread scientific interest in the 
measurement of slight density changes and density 
differences as well as in the measurement of absolute 
density itself. The Cartesian diver system is partic- 
ularly well-suited to determination of changes and 
differences because of its ability to provide high 
precision data on samples of small size. In problems 
such as determination of density differences be- 
tween nominally identical materials but having 
slightly different physical properties (such as_ be- 
tween hard drawn and annealed wire), or determina- 
tion of small density changes when a material is 
exposed to a density-modifying process (such as UV 
irradiation), homogeneity of the sample is usually 
of prime importance. Homogeneity can frequently 
be verified or inhomogeneity evaluated in samples of 
a few hundred milligrams, but it is frequently im- 
practical or impossible to do so 1n samples of several 
grams. Some density modifying processes can pro- 
vide homogeneous density changes only in small 
samples, as, for example, processes depending upon 
x-ray beams that have limited diameters or pene- 
trating powers. 


byt f py 
t of ¢ 


Vst: ] of Interest 


Two measured values of By 


Two calculated values of px 


Average 


PX 


Change of px during modification 0. 000204 


The above data show the agreement in B, to be 
about 0.3 we, which is somewhat poorer than average 
diver operation. Nevertheless, although the change 
in density associated with the process under exami- 
nation amounted to only 5:10°, the system reported 
this change to agreement of about 7's percent. 

This experiment was performed upon a crystal 
of only 250mg. With larger samples, the reproduci- 
bility of By would remain about constant, hence the 
reproducibility of values of density would be much 
improved. 


6. General Comments 


Although our work was not primarily diiected to 
evaluation of diver performance as a function of 
diver characteristics, the following comments may 
be of interest: 


(A) We were never successful in working with 
Pyrex divers. Although the ones we tried out were 
always thoroughly annealed after fabrication, we 


In measuring density changes in an object sub- 
jected to a density modifying process, the diver 
system is largely free of the requirements for high 
accuracy standards. If the two standards are cut 
from the same piece of parent material from which 
the object of interest is cut, and to about the same 
geometry, then we are justified in assuming that the 
handbook value of the density of the material applies 
to the standards absolutely. 

The diver was recently called upon to measure the 
change in density of a single crystal sample of rutile 
weighing less than 260 mg, when the stoichiometry 
was altered by thermal reduction. The standards 
used were cut from adjacent positions in the same 
boule from which the object of interest was cut. 
All three crystals were carefully weighed and two 
independent measurements on the apparent mass in 
water of the unknown, By, were made based on the 
assumption that the density of the standards was 
4.249 g/cm’. 

The crystal of interest was partially reduced as 
called for in the experimental plan and was then re- 
weighed. 
upon By. 


Two additional observations were made 
The data of interest is shown below: 


Before modific 


After modification 


») 


0. 2586061 ¢ 0. 2585972 


0. 1977429 0. 197426 0. 1977389 0. 1977388 


4. 248973 4. 248952 4. 249169 4. 249162 


4. 248962 4. 249166 


never were able to achieve stability of much better 
than 10 pe. 

(B) The sensitivity of the sytem may be crudely 
defined as (—)AB/AP», that is the indicated change 
in apparent mass per unit change in equilibrium 
pressure. Defined thus, the sensitivity of the diver 
reported on was about 2 ug/g change in Py. We 
experimented upon highly compressible forms in 
which sensitivity was 10 ug/g, and also upon struc- 
tures with compressibilities so small that sensitivity 
approached 0.003 ug/g. We found that the re- 
producibility decreased as the divers had extremely 
small or extremely large values of sensitivity. For 
the extremely sensitive structures we experienced 
difficulty in obtaining hydrostatic balance, and for 
the relatively insensitive structures, zero drift was 
most annoying. 

(C) The system is quite sensitive to vibration and 
the mechanical shock associated with closing doors, 
etc. Spaepen also reported this vibration sensitivity 
and described his efforts to minimize it. Our experi- 


a/o 
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ment was conducted on a massive pier which had 
previously been used to support a high sensitivity 
balance used in the NBS mass calibration program. 
A vibration survey conducted prior to the experiment 
indicated that at 30 Hz, the vibrations were about 
610g, decreasing somewhat at night. Our best 
work was done at night and on weekends. 

(D) Many experimenters working with other 
buovant devices have reported difficulty in achieving 
smooth “lift-off of the buoyant structure from an 
alining cone at the bottom. We found that a blunt 
tip at the bottom of the quartz diver moved nicely 
put of a chromium-plated brass cone whose radius 
Was somewhat greater than the quartz. We experi- 
mented with cones of both Pyrex and sapphire with 
poor results. 

(E) In a private communication Spaepen informed 
us that a newly fabricated diver is rarely stable 
that a one or two months aging period is called for. 
This has also been our experience. We believe that 
annealing the quartz diver would probably reduce 
this aging period, 

(Ff) A major uncertainty in all hydrostatic weigh- 
ing is caused by the adherence of gas to the surfaces 
of the load. An air bubble large enough to displace 
1 ug of water is easily detecied visually and easily 
removed. An equal amount of gas dispersed among 
several hundred (or thousand?) bubbles would be 
impossible to detect (without @ priort knowledge of 
the density of the material being hydrostatically 
weighed). This could lead to serious systematic 
errors in the values of density caleulated from ob- 
served data. In order to randomize the surface 
trapped air, we always remove all loads (both 
standards and unknowns) between each hydrostatic 
experiment and thoroughly dry and reclean them. 
A stable volume of gas adhering to the surfaces of the 
diver body is less annoying, since its chief effect is to 
change the values of C, and G which are determined 
in each experiment. 

(G) In changing diver loads the portion of the 
tools used which come into contact with the loads 
are kept continuously immersed in water, never 
allowing them to come into contact with air. In this 
way we minimize the likelihood of transferring air 
from the tools to the loads. 

(H) Generally, a perfectly cylindrical diver body 
will be insufficiently compressible for practical diver 


operation. Our method of achieving the desired 
compressibility is to measure AG/AP on a newly con- 
structed diver in its initial evlindrical form. We then 
partially flatten the cylindrical walls by flame treat- 
ment in small increments, measuring AG/AP after 
each operation. This process is continued until the 
Figure 1 shows this 


desired sensitivity Is attained. 
partially flattened section. 
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Cryogenic Behavior of Selected Magnetic Materials 
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Institute for Materials Research, National Bureau of Standards, Boulder, Colo. 
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Measurements of a-c core loss and d-c magnetic properties including coercive force, 
residual induction, and hysteresis loss are reported at temperatures of 273, 195, 76, and 4 °K 
for a variety of commercial magnetic Fe-Si and Fe-48 percent Ni alloys. Of the materials 
tested the high silicon content (3 to 4.45 percent Si) Fe-Si alloys were found to be the least 
temperature sensitive. Hysteresis loss increases between 273 and 4 °K of up to 100 percent 
were measured for the iron-nickel materials 


1. Introduction 2. Materials Tested 


An increasing use of electric motors, transformers, Tables 1 and 2 list the various Fe-Si and Fe-Ni 
and other electrical equipment operating at cryo- | samples tested together with related pertinent 
genic temperatures has created the need for engi- | information concerning sheet thickness, type of 
neering data on the low-temperature magnetic | core, and chemical composition. The identification 
properties of ferromagnetic materials. To fill this | numbers given will be referred to throughout the 
need, a-c and d-c magnetic properties were measured | text. Materials 110 through 118 were kindly supplied 
at 4, 76, 195, and 273 °K for a variety of commercial | by United States Steel Corporation and materials 
Fe-Si and Fe-Ni alloys. 100 through 102 and 200 through 213 by Armco 

Steel Corporation. Materials 100 through 102 are 
identified by the trade name, “Oriented T,’’ 200 
through 203 by “48 Orthonik,” and 210 through 213 
; *Present address: Colorado Earth Science Laboratories, Inc., Boulder, Colo by “48 Ni.” 


TaBLe 1. Jron-silicon alloys tested at cryogenic temperatures 


Composition, percent 
Material AISI rype valance primarily tron) | Grain orientation | Thickness lype of core Identification 
Number product I 


umber 
Si Cc Mn 


‘hin grain-oriented Oriented. __- Tape-wound, B* 
nominally 3 percent Fe- Oriented _. 2 Tape-wound, 

Si alloy. Bd y 

43 . 03 

. 05 


3 
B 
Oriented. __. : Tape-wound, B 
noriented . . . 02 Laminated ” 


{ { 
43 Unoriented ; Laminated—l 

36 . 04 Unoriented _ . . Laminated—l 

, 04 Unoriented ; Laminated—l 

. 017 : Unoriented _ . ; Laminated—U, 

d I I 

if 

I { 

I 


Samm AZAR eOeon ese 


04 ‘oriented ; Laminated 
. 03 ‘norieuted Laminated 
. 021 noriented Laminated 
043 ‘noriented Laminated 


= 


boxed. **U =unboxed,. H R=hot rolled. y cold rolled. 


TaBLe 2. Jron-nickel alloys tested at cryogenic te mperatures } 4) Measuring Technique 


Nominally 48 percent Ni content i 
3.1. Equipment 


Material Thickness Type of core Identifica- A standard voltmeter-wattmeter method [1]! was 
t he . a bh e 
‘on number | used to determine a-c core loss values. The magnetic 

induction in the core is determined by measuring the 

et a : meres, . 2 secondary voltage with an a-e differential voltmeter 

FeNi Oriente = wound, B 2 . : x, ° 

FeNi Oriented jaa, © 2 calibrated to 0.5 percent accuracy. Induction levels 

FeNi Oriented Laminated, B : 

FeNi Unoriented Tape-wound, B 

FeNi Unoriented Tape-wound, B 

FeNi Unoriented Laminated, B 

FeNi Unoriented Laminated, B 


1 Figures in brackets indicate the literature references at the end of this paper. 
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values below that which resulted in 
greater 3rd harmonic distortion as de- 
a wave analyzer across the secondary 
windings. <A high input impedence electronic watt- 
meter was employed thereby making instrumental 
loss corrections negligible. Direct current hysteresis 
loops were plotted semiautomatically on an 2-y 
chart recorder. The magnetic induction was _ re- 
corded using an integrating d-c amplifier. The 
magnetizing force was calculated from the voltage 
drop produced by the magnetizing current flowing 
through a standard resistor. 


were limited to 
percent or 
termined by 


3.2. Samples 


cryogenic testing weigh 
approximately ‘x Ib and are considerably smaller 
than that specified by the ASTM standards [2]. 
This was necessary to prevent excessive liquid 
evaporation and the use of inconveniently large 
Dewars. All samples are in the shape of toroids. 


The samples used for 


The tape-wound and laminated samples vary slightly 
mean 
This 
distribution 
sample. 


in size. In all cases, however, the ratio of 
diameter to radial width is 8 to 1 or greater. 
ensures a reasonably uniform flux 


throughout the cross-sectional area of the 


Total core loss at 


Core loss value at 273 °K 


perature 


The mean diameter of the cores is approximately 
2 in. For most materials two samples of each 
were tested. 

The oriented iron-silicon material and all iron- 
nickel samples were placed in either anodized alu- 
minum or bakelite boxes prior to winding the second- 
ary and primary coils. These loose-fitting boxes 
prevent coil winding stresses from affecting the 
magnetic properties. Since the unoriented grades of 
Fe-Si are seldom boxed in normal use, these samples 
were prepared by wrapping the laminated toroid 
with 0.004-in. Teflon tape to provide electrical 
insulation between windings and core material. 
The primary and secondary windings were then 
placed directly atop the Teflon tape. To determine 
the effect of winding stresses on the magnetic proper- 
ties, four materials, 111, 112, 114, and 117, were also 
tested in boxed form. The large increase in core 
loss for the unboxed cores compared to the boxed 
cores at low temperature can be attributed to the 
detrimental stress effects produced by the differential 
contraction between the iron-core material and the 
copper windings. Because of the wide use of the 
unoriented grade Fe-Si in components where little or 
no emphasis is given to eliminating winding or 
clamping stresses, core loss data for the unboxed as 
well as selected boxed materials are shown in table 3. 


K and the ratio of core loss measured at te mperature T to that measured at T 973 °K 


is given in units of watts/pound 


Peak induc 





Material 


Fre- 
quen- 


rer 
pera 
ture 


ik induction kG 


10 


Fre- | Tem- Peak induc 
Material | quen- 
cy 


0. 0108 
28 

1. 69 

1. 93 


With the exception of sample 11S all materials 
were stress-relief annealed by the fabricators follow- 
ing their normal annealing procedures. Sample 118 
Was tested in both unannealed, as fabricated, and 
annealed conditions. Both unannealed and an- 
nealed samples were boxed to prevent the extraneous 
effects of winding stresses. The temperature be- 
havior of these two samples will be discussed later 
in relation to theories of coercive force. 

Internal heating was, in some instances, great 
enough to raise the sample temperature above that 
of the environmental bath. The temperature rise 
was determined using a Au-Co versus copper thermo- 
couple. Core loss values presented in table 3 have 
been interpolated where necessary to give loss values 
at the bath temperatures listed. 


4. Alternating Current Loss Data 


Alternating current core loss for each sample was 
measured at 4, 76, 195, and 273 °K for various peak 
inductions from 2 to 18 kg.? All materials were test- 
ed at 60 c/s (hertz). Where the type of material 
warranted it, measurements were also made at 400 
and 1000 ¢/s. At 60 ¢/s the upper induction level is 
limited by waveform distortion. At 1000 ¢ sit is lim- 
ited by the breakdown voltage of the magnet wire 
insulation. In the interest of saving space, only four 
graphs (figs. 1-4) showing core loss versus tempera- 
ture are presented. The bulk of the data are pre- 
sented in table 3. 

From the data given in figures 1 to 4 and table 3 
it is noted that the high silicon content (3 to 4.45% 
Si) Fe-Si alloys are the least sensitive to temperature 
variations. This is primarily a result of these alloys 
having a lesser temperature dependence of electrical 
resistivity resulting in a smaller increase in eddy 
current loss. Although there has been some criticism 
in the past of separating total core loss into hysteresis 
loss and eddy current loss components, we find that 
it is helpful, at least qualitatively, to speak in terms 
of a separate hysteresis and eddy current loss. 


rhe following relationships may be used to obtain equivalent value 
MKSA system of units 
1 kilogauss 0.1 tesla 
1 oersted 79.58 ampere turns/meter 
1 watt/pound= 2.204 watt/kilogram 
linch 0.02510 meter 
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loss at 60 c's for 


Material 102. 


Total core arious peak inductions 


For instance, upon examining r the data for material 
102 in table 3, it is noted that at 60 ¢/s and for low 
peak inductions there is a loss increase of approxi- 
mately S percent upon cooling from 273 to 4 °K. 
At 400 ¢/s and 1000 c/s however, the corresponding 
loss increase is approximately 20 percent. This 
difference may be explained by assuming the greater 
part of the total loss at 60 ¢/s to be due to a hysteresis 
loss while at 400 and 1000 ¢/s an increasing per- 
centage of the total loss is due to eddy current losses 
which, to a first approximation, may be expected 
to vary inversely as the material resistivity. From 
the resistance data given in figure 5 it is seen that for 
material 102, 2,/R2-3=0.82. This gives an expected 
eddy current increase of approximately 22 
percent which is in good agreement with the total 
core loss data at 1000 ¢/s. By similar arguments 
involving separate hysteresis and eddy current loss 
components a satisfactory comparison of the tem- 
perature dependence of core loss can be made among 
materials 100, 101, and 102 which are similar except 
for tape thickness. 

As a result of this study it may be said that, in 
general, a-c will increase upon cooling 


from 273 to 4 °K. The amount of increase varies 


loss 


core loss 
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Total core loss at 60 c/s for 


Vaterial 114 


FIGURE 2. arious peak nductions: 


inbored. 


greatly from one material to another, 
being from nearly no increase to almost 
for the materials tested. The 
the above statement is material 118. As seen in 
table 3 this sample, in the unannealed condition, 
showed a loss decrease of approximately 4 percent 
at 60 c/s. This decrease with decreasing tempera- 
tures is also evident in the d-c hysteresis loss measure- 
ments and will be discussed further in the section 
on d-c properties. 


the range 
200 percent 
one exception to 


5. Direct Current Properties 


Direct current hysteresis loops at peak inductions 
of 5, 10, and 15 kG were measured for the same 
samples used for a-c measurements. Values of 
the coercive force, H,, and the ratio of residual 
induction to peak induction, B,/B,,, were measured 
directly from the loops and should be accurate to 
within approximately 1 percent. The d-c hy- 
steresis loss, P,, was determined from planimeter 
measurements of the loop area and should be 
accurate to about 1 to 2 percent except for induc- 
tion levels which required a change of scale in the 
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FIGURE 3. oss at 60 c's for various peak inductions 


Material 114 boxed 


H-axis. 


In the latter cases the loss values may be 
in error up to approximately 5 percent. 

Inasmuch as the winding stresses are known to | 
influence the hysteresis loss temperature dependence, 
only the boxed samples will be considered here. 
The effects observed are therefore of purely magnetic 


origin. The results for the Fe-Si and Fe-Ni alloys 
are presented separately. 


TABLE 4. Coercive force (H, 


5kG 


ergicm/cycle 
537 


290 
163 
661 
541 
621 
386 
118 unannealed 3 .6 517 
118 annealed 3 5 281 


), ratio of residual induction to peak induction (B,/ Bm), 





+ 


y Peak Induction = 1/4 kG 
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Material 213. 


5.1. Fe-Si Alloys 


The temperature dependence of the coercive 
force, residual induction, and hysteresis loss is 
shown graphically in figures 6 through 10 with the 
results normalized to 1.00 at 273 °K. For com- 
parison purposes the absolute values of /7,, B,/B,,, 
and P, at 273 °K are tabulated in table 4. 


and hysteresis loss (P},) at 273 


wd 


10kG 15kG 


Pr 


erg/cm/cycle 

1956 

9S 

507 
1937 
1603 
2053 
1163 
1664 

821 
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iterials listed in tables 1 and 2 
117; (3) 116 5 


4) 113 and 100, 101, 102 114; (6) 112; (7) 110; (8) 200 


a. Hysteresis Loss 


The magnetization process in ferromagnets is 
commonly described in terms of domain wall motion 
and domain rotation processes. By measuring the 
d-c hysteresis loops at peak inductions of 5, 10, and 
15 kG it was intended to proceed from an induction 
level (5 kG) where comparatively simple domain 
wall motion processes predominate to an induction 
level (15 kG) where domain rotations and the more 
complex closure domains [3] become increasingly 
important. It is of course realized that this some- 
what arbitrary choice of peak induction levels will 
distinguish between the different magnetization 
processes to varying degrees for the various materials 
tested. ‘To illustrate this it would be well to con- 
sider materials 100 and 102. The hysteresis loops 
at a peak induction of 15 kG are shown for these 
two materials in figures 11 and 12, respectively. It 
is evident from the loops that the magnetization of 
material 102 has proceeded primarily by simple 
domain wall motion. On the other hand, the curve 
for material 100 that the magnetization 
process is well into the region of domain rotations 
and the associated closure domain wall movements. 
This influence of magnetization process on the tem- 
perature dependence of the hysteresis loss is shown 


suggests 
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Grain-oriented materia . 101; 102, Materials identified in table 1 


in figure 7. Here it is noted that at 15 kG, material 
100 exhibits a temperature independent hysteresis 
This corresponds to Chen’s statement [3] that 
at the higher temperatures the decrease in energy 
loss associated with irreversible domain wall move- 
ment may be completely compensated by an increase 
in closure domain loss. At 15 kG material 102 still 
shows an increased hysteresis loss at 4 °K. The 
increase, however, is only about 5 percent compared 
to 10 percent at 10 kG and approximately 13 percent 
at 5 kG. Thus it indicates that the magnetization 
of material 102 is just starting to enter the region of 
domain rotation at 15 kG. It would be expected 
that at some higher peak induction material 102 
would exhibit a temperature behavior comparable to 
the behavior of material 100 at 15 kG. This in facet 
is the case at 18 kG where the hysteresis loss for 
material 102 was measured to be approximately 
2 percent lower at 4 °K compared to the 273 °K 
value. Ignoring material 118 for the moment it is 
found that a similar situation exists for the unoriented 
alloys as shown in figure 9. While there is a hys- 
teresis loss increase of approximately 4 percent at 
76 °K for both 5 and 10 kG the hysteresis loss 
decreases by approximately 2 percent for a peak 
induction of 15 kG. 


loss. 
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b. Coercive Force 


Upon comparison of the curves for the tempera- 
ture dependence of coercive force and hysteresis loss 
it is found for materials 100-102 that for peak indue- 


tions of 5 and 10 kG the temperature behavior is 
nearily identical. Thus the hysteresis loss behavior 
is determined almost entirely by the temperature 
dependence of the coercive force. At 15 kG (18S kG 
for material 102), there is no such correspondence. 
This is due to the effect of closure domain loss as 
discussed previously. It should be observed, how- 
ever, that the increase in coercive force at 4 °K is 
very nearly the same for both the 10 and 15 kG 
levels of induction. 


c. Residual Induction 


As discussed by Chen the residual induction may 
be expected to increase with decreasing temperature. 
The primary reason for this is the temperature 
dependence of the magnetocrystalline anisotropy, K, 
in the expression for the wall energy. Inasmuch as 
this model is concerned with the nucleation of new 
reverse domains, the expression for wall energy 
would be expected to give an indication of the 
temperature dependence of the residual induction 
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Boxed unoriented material 111 s 112 
< 11S annealed. Materials identified in table 


only at induction levels approaching saturation where 
all the original domain walls have been eliminated. 
Indeed, for the samples shown in figure 10, with the 
exception of 118, the residual induction is essentially 
temperature independent at the lower peak induc- 
tions of 5 and 10 kG. 


d. Material 118 


As can be seen from the curves for 7,, B,, and P,, 
material 118 exhibits a behavior different from the 
other iron-silicon alloys tested. For this reason it is 
presented separately. We have no complete ex- 
planation for the temperature behavior of this 
material but feel it is probably associated with the 
vanishing magnetostriction constant at low temper- 
atures for a silicon content of 4.45 percent [4]. 

As discussed by Chen the current theories indicate 
the coercive force to be a function of » lor K IT. or 
both, where \, is the magnetostriction constant, /, 
is the saturation magnetization, and Ais the magneto- 
crystalline anisotropy. Whereas Chen considered 
the energy term associated with magnetostriction to 
be negligible compared to the crystal anisotropy 
term, this is not the case for 4.45 percent silicon 
iron. In the case of the unannealed, “as processed’”’ 
sample, the material is undoubtedly severely stressed. 
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Boxed unoriented material 111; A 112 114; Y 117; @ 118 unannealed; X 
118 annealed. Materials identified in table 1 
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Figure 11. Hysteresis loop at a peak induction of 15 kG. 


Material 100, 


Taking the internal stress to be of the order of the | The decrease in 7, shown in figure 8 is thus attributed 
yield stress (4.9 10° dyn/em?) and Aioo=2.010~, | to the decrease in the magnetostriction constant 
the energy associated with magnetostriction is | upon cooling below 273 °K. 

9.8 10* ergs/em* compared to a crystal anisotropy If this interpretation is correct the removal of 
energy of 30> 10* ergs/em* at room temperature. | internal stresses by annealing would be expected to 
Clearly the magnetostriction energy is not negligible. | result in a temperature dependence of coercive force 
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Ficure 12. Hrysteresis loop ata peak induction of 15 k@. 
Material 102, 


determined by the crystal anisotropy energy only. 
As shown for the annealed sample this result was 
partially realized; the coercive force at 15 kG in- 
creased slightly with decreasing temperature. It 
should be noted that the coercive force at 5 kG and 
10 kG still shows a small decrease at the lower tem- 
peratures. The decrease is, however, much less 
than forthe unannealed sample. It is possible that 
annealing the material at 1425 °F for 24 hr does not 
completely remove the internal stresses, thus leav- 
ing a small dependence on the magnetostriction 
energy. Annealing at higher temperatures was not 
attempted since it was feared that recrystallization 
might occur thereby introducing a new variable of 
crystal orientation [5]. 

The behavior of the residual induction for both 
the unannealed and annealed material is shown in 
figure 10. The unannealed sample shows a greater 
temperature dependence than the annealed material. 
For a material whose magnetostriction decreases 
with decreasing temperature this would not be 
expected to occur on the basis of the wall energy 
equation mentioned previously. Barring the possi- 
bility of recrystallization occurring during the an- 
nealing treatment one must conclude that it is an 
oversimplification to treat the residual induction 
in terms of the wall nucleation process only. 
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Figure 13. Hysteresis loops. 


Material 200. 


5.2. Fe-Ni Alloys 


The iron-nickel alloys tested are less well-behaved 
than the iron-silicon alloys discussed in the previous 
section. Their behavior is not completely under- 
stood at present. In view of their potential use at 
cryogenic temperatures however, and the large 
changes in the magnetic properties which were 


found to occur, some results are presented here 


without explanation. The curves are 
primarily for design and application use. 

Hysteresis loops for a variety of materials are 
presented in figures 13 through 20. The curves 
shown in figure 17 are for a material commercially 
supplied as 0.002 in. 48 Ni. However, metallurgical 
investigation indicated a very small grain size which 
is not considered to be typical for this material. 
The curves are presented in view of the strange 
behavior evident at 76 °K. It is not known whether 
such a temperature behavior would be the rule or 
the exception for commercially supplied material 
of this type. A laboratory prepared sample, of 
0.002 in. 48 Ni, supplied by Armco Steel, was also 
tested. The results are shown in figure 18 where 
it is seen that the kinked loop behavior is absent. 


intended 


6. Summary 


The a-c and d-c magnetic properties of selected 
Fe-Si and Fe-Ni alloys have been measured at cryo- 
genic temperatures. The results are presented in 
tables and graphs to provide engineering design data. 
Of the Fe-Si alloys, the thin grain-oriented material 
would be best suited for cryogenic use. In addition 
to its superior room temperature properties, it 
exhibits a considerably lower loss increase at cryo- 
genic temperatures than do the thicker laminated, 
low silicon content alloys. The Fe-Ni alloys are in 











To 25 


Oersteds 


B - KILOGAUSS 





FiGurRE 14. Hysteresis loops. 
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FIGURE 17. Hysteresis loops. 


Material 210; commercial preparation 
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Figure 18. Hysteresis loops. 


Material 210; laboratory preparation. 
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FIGURE 19. Hysteresis loops. 


Material 211. 
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Figure 20. Hysteresis loops. 
Material 213. 
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general more temperature sensitive than the Fe-Si 
alloys. Hysteresis loss increases on the order of 
100 percent have been measured for some materials 
at 4 °K. The large decrease in residual induction 
at low temperatures may make this material less 
desirable for certain applications. 
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Observations of dislocations and surface features in corun- 
dum crystals by electron transmission microscopy, D. J. 
Barber and Nancy J. Tighe, J. Res. NBS 69A (Phys. and 
Chem.), No. 3 (May-June 1965), 70 cents. 

Crystals of corundum, grown by the flame-fusion method, 
have been chemically thinned and examined by electron 
transmission microscopy. Some specimens were etch-pitted 
after thinning, and the relationship between etch pits and 
dislocations was investigated. Within subgrains there is a 
close agreement between dislocation and etch pit densities; 
this agreement breaks down along many boundaries. The 
etching behavior of crystals containing Cr,O;, TiO. and 
other impurities has been examined. Eteh-tunnelling along 
grown-in dislocations is commonly observed in the impurity- 
doped materials; it is particularly pronounced in ruby. In 
heat-treated crystals containing TiO:, local differences in 
solid solubility can lead to the formation of etch hillocks. 
Thermally etched surfaces have also been characterized, and 
gas atmospheres are shown to affect the surface topology. 
Evaporation experiments in the electron microscope are 
reported. 


Phase velocities and attenuation distances in the ionosphere, 
D. R. Croley, Jr., and B. 8. Tanenbaum, “Radio Science J. 
Res. NBS/USNC-URSI 69D, No. 6, 819-825 (June 1965) 
$1.00. 

Plane wave phase velocities and attenuation distances in an 
infinite, homogeneous, partly ionized gas composed of an 
interacting mixture of neutral molecules, ions, and electrons 
in the presence of a uniform magnetie field are determined 
for three sets of plasma conditions corresponding to the 80, 
100, and 300 kilometer levels of the ionosphere. Propaga- 
tion in the parallel direction and the perpendicular direction 
relative to the magnetic field in the frequency range from 
approximately 10-5 ¢/s to 10° ¢/s is considered. — 


International comparison of atomic frequency standards via 


VLF radio signals, A. 
Radio Science J. 
(July 1965), $1.00. 
A study was made of data obtained over an 18-month period 
(July 1961 to December 1962, inclusive) on the comparison 
of atomic frequeney standards located in seven laboratories 
in the United States, Kurope, and Canada, using the VLF 
signals of GBR (16 ke/s), Rugby, England and NBA (18 ke/s), 
Balboa, Canal Zone. Each laboratory observes the accumu- 
lated difference in phase over a 24-hour period (the same for 
all laboratories, or nearly its own standard 
(either laboratory or constructed) and the 
received VLF signal. A statistical analysis was designed to 
separate the observations at each laboratory into three com- 
ponents: (a) long-term 


H. Morgan, EF. L. 


Crow, B. E. Blair, 
Res. 


NBS/USNC-URSI 69D, No. 7 


] 


so) between 


commercially 


mean differences among the atomic 
standards; (b) estimates of the standard deviatioas, &;, at each 
receiving station; and (c) estimates of the transmitter stand- 
ard Each a; includes receiver fluctuations, 
propagation effects peculiar to the path, and measurement 
uncertainties; 7 includes the transmitter 
propagation effects common to all paths. 


deviations, 7. 


fluctuations and 
The study shows that @; at each receiver varied from a low 
or 0.39107! units of fractional frequency (that is, 0.39 
parts in 10") (GBR data) at LSRH to a high of 1.97 10-" 
(GBR data) at NRC with an average for all stations of 1.01 

10-!° measured against GBR and 0.99 10-!° when measured 
against NBA. Also, the average 7 for GBR is 1.26 107-12 


and for NBA is 0.68 10-". Finally, it is shown that: (1) 
the means of the frequencies of the seven individual atomic 
standards laboratories agreed with the grand mean of these 
seven standards to within +2 parts in 10'° for the 18-month 
period, and (2) the laboratory-type standards agreed with 
their grand mean to within +1 part in 10”. 


Control of WWV and WWVH standard frequency broadcasts 
by VLF and LF signals, B. E. Blair and A. . Morgan, 
Radio Science J. Res. NBS/USNC-URSI 69D, No. 7, 
(July 1965), $1.00. 

Since 1961 the NBS VLF and LF signals have improved the 
calibration and frequency control of the WWYV (Maryland) 
HF broadcasts. Similarly, better control of the WWVH 
(Hawaii) HF broadcasts was achieved in early 1963 by 
monitoring the NBS VLF broadeasts in terms of the WWVH 
control oscillator. In mid 1963 WWVL (20 ke/s) and WWVB 
(60 ke/s) were relocated from two sites near Boulder, Colo- 
rado, to a single site near Ft. Collins, Colorado, and the 
transmitter power for both broadcasts was increased several 
fold. These higher powered broadcasts resulted in more 
precise control of both HF broadcasts. Through the VLF 
and LF signals the 24-hour average frequency values of 
WWY are related to the United States Frequency Standard 
(USFS) within a few parts in 10". 

This paper describes the NBS low frequency broadcasts, the 
method of using them to control and calibrate the HF 
broadcasts, and gives an analysis of the precision of frequency 
control obtained at WWV over a 21-month period. An ap- 
pendix discusses the short-term phase stabilities and diurnal 
phase shifts observed in the low-frequency signals at WWV 
and WWVH, and examines the accuracy-limiting effects of 
propagation path characteristics and background noise levels 
in such received signals. 

Theory of coil antennas, T. Padhi, Radio Science J. Res. 
NBS/USNC-URSI 69D, No. 7, (July 1965), $1.00. 

In this paper, a method is presented by which the distribution 
of current on some structurally simple coil or multiturn loop 
antennas may be obtained. The input admittances of un- 
shielded and shielded coils are determined and their operation 
as receiving elements is considered. 


Their construction, maintenance, and char- 


Hamer, NBS Mono. 84 (Jan. 15, 1965), 


Standard cells. 
acteristics, W. J. 
35 cents. 

This Monograph contains information on the construction, 
maintenance, and characteristics of standard cells. The 
effects of temperature, pressure, electric current, light, shock, 
and vibration on standard cells are discussed. A history of 
the realization and maintenance of the unit of electromotive 
force is also included. A record of international comparisons 
of the unit of electromotive force is presented as well as 
information on the constancy of the National Reference Group 
of Standard Cells. 


Effects of thermal shrinkage on built-up roofing, W. C. Cullen, 
NBS Mono. 89 (Mar. 4, 1955), 10 cents. 
The effects of thermal movement of bituminous built-up roof 
membranes are related to some common failures of built-up 
roofing observed in field exposures. A procedure for deter- 
mining the amount of thermal movement of built-up mem- 
branes is described. Data are given for the thermal move- 
ment of various bitumens and reinforcing felts used in the 
construction of built-up roofs and for the composite mem- 
branes over a temperature range of +30 to —30 °F. The 
data obtained are related to field experience and suggestions 
| are presented which will prove beneficial in reducing the 
| incidence of built-up failure due to thermal movement. 
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Calibration of liquid-in-glass thermometers, J. F. Swindells, 
NBS Mono. 90 (Feb. 12, 1965), 25 cents. Supersedes NBS 
Circ. 600. 

This Monograph, which supersedes Circular 600, contains 
information of general to both manufacturers and 
users of liquid-in-glass thermometers, as well as those who 
wish to calibrate thermometers or submit them to the Na- 
tional Bureau of Standards for calibration. Instructions are 
provided for applicants requesting calibration services, and 
the techniques and equipment used in tne calibration pro- 
cedures are described. Important elements of thermometer 
design are discussed, and factors affecting the use of common 
types of liquid-in-glass thermometers are included together 
with tables of tolerances and reasonably attainable accuracies. 
The calculation of corrections the temperature of the 
emergent stem is given in detail for various types of ther- 
mometers and conditions of use. 


interest 


tor 


Automatic indexing: A state-of-the-art report, \I. E. 
NBS Mono..91 (Mar. 20, 1965), $1.50 

A state-of-the-art survey of automatic indexing systems and 
experiments has been conducted by the Research Information 
Center and Advisory Service on Information Processing, 
Information Technology Division, Institute for Applied Tech- 
nology, National Bureau of Standards. Consideration is first 
given to indexes compiled by or with the aid of machines, 
including citation indexes. Automatic derivative indexing is 
exemplified by key-word-in-context (K WIC) and other word- 
in-context techniques. Advantages, disadvantages, and pos- 
sibilities for modification and improvement are discussed. 
Experiments in automatic assignment indexing are sum- 
marized. Related research efforts in such areas as automatic 
classification and categorization, computer use of thesauri, 
statistical association techniques, and linguistic data process- 
ing are described. A major question is that of evaluation, 
particularly in view of evidence of human _ inter-indexer 
inconsistency. It is concluded that indexes based on words 
extracted from text are practical for many purposes today, 
and that automatic assignment indexing and classification 
experiments show promise for future progress. 


Stevens, 


The examination of weighing equipment, M. W. 
R. W. Smith, VBS Handb. 
sedes H3?. 

This Handbook is primarily directed to the presentation of a 
coordinated series of step-by-step Examination Procedure 
Outlines (EPO’s) for weighing equipment, recommended for 
adoption and use, as minimal requirements, by weights and 
measures agencies, commercial service agencies, weighing 
equipment owners and operators, and manufacturers of 
weighing devices. Supporting information embraces discus- 
sions on related topics such as reference and field standards 
report forms, tolerances, weighing principles, and elements of 
scale construction and performance. Extensive weights and 
measures tables—basic, interrelation of units, and equiv- 
alents—and a list of references for further study are appended. 
An alphabetical index is supplied. 


Jensen and 


94 (Mar. 1, 1965), $3.00. Super- 


Standard frequency and time services of the National Bureau 
of Standards, NBS Misc. Publ. 236 (Jan. 15, 1965), 15 cents. 
Supersedes editions of Dec. 1, 1960 and July 1, 1961. 

Detailed descriptions are given of eight technical services 
provided by the National Bureau of Standards radio stations 
WWV, WWVH, WWVB, and WWVL. These services are: 
1. Standard radio frequencies; 2. Standard audio frequencies; 
3. Standard musical pitch; 4. Standard time intervals; 5. Time 
signals; 6. UT2 corrections; 7. Radio propagation forecasts; 
and 8. Geophysical alerts. In order to provide users with the 
best possible services, occasional changes in the broadcasting 
schedules are required. This publication shows the schedules 
in effect on January 1, 1965. Annual revisions will be made. 
Advance notices of changes occurring between revisions will be 
sent to regular users of these services upon request. 

Standard Reference Materials: Sources of information, J. L. 
Hague, T. W. Mears, and R. E. Michaelis, NBS Misc. Publ. 
260-4 (Feb. 1965), 20 cents. 

This annotated listing of of Standard Reference 
Materials is prepared with references to all types of materials 


sources 


used to standardize analytical, physico-chemical, and engi- 
neering methods. The range of Standard Reference Materials 
included runs from high-purity substances and carefully 
analyzed metals, alloys, and rocks to materials of ind#finitely 
known composition for standardizing a single phenomenon. 


Standard Reference Materials: Accuracy of solution x-ray 
spectrometric analysis of copper-base alloys, R. Alvarez and 
R. Flitsch, NBS Misc. Publ. 260-5 (Mar. 15, 1965), 25 cents. 
X-ray fluorescence spectrometry has been investigated as an 
independent method for the analysis of certain major con- 
stituents in NBS standard reference materials. The deter- 
mination of copper and zine in cartridge and aluminum brass 
samples which were to be certified as standard materials, 
was selected for study. After a preliminary analysis, a solu- 
tion of the unknown and a closely-matched synthetic stand- 
ard were compared without the use of an internal standard or 
added control element. To achieve reliable results, close con- 
trol of the type and concentration of acid was found neces- 
sary, and, for cells employing films as windows, the position 
of the film was critical. The temperature rise of the cell 
during irradiation was 0.1 °C, which did not contribute signifi- 
cantly to error. Typical coefficients of variation obtained for 
individual results in a group of 8 runs were 0.3°% for copper 
and 0.4°;, for zine. The agreement between chemical and 
x-ray results was within 0.3¢% of the amount present for both 
elements. 


Optical scintillation; a survey of the literature, J. R. Meyer- 
Arendt and C. B. Emmanuel, NBS Tech. Note 225 (Apr. 4, 
1965), 70 cents. 

In this Technical Note, main emphasis is placed on providing 
the reader with an exhaustive survey of the literature, 
covering the field of effects of atmospheric refraction on the 
propagation of electromagnetic radiation at optical fre- 
quencies. 

One may distinguish systematic, regular, or normal refraction 
on the one side from random refraction on the other. The 
former can be predicted theoretically, using various types 
of atmospheric models; the latter requires analysis by sta- 
tistical methods. 

Numerous observational, experimental, and theoretical as- 
pects of random refraction, that is, of scintillation in its 
widest sense, are discussed. The following topics are dealt 
with: refraction in plane and in spherically stratified media, 
refractive index variations, radio refraction, scintillation as a 
function of aperture size, zenith distance, site location, dis- 
persion, and meteorological conditions, the frequency spec- 
trum of scintillation, terrestrial scintillation, image distortion 
and contrast reduction, refraction and diffraction theories 
of scintillation, autocorrelation analyses, radio star scintilla- 
tion, and coherence problems. Questions of atmospheric 
scattering, absorption, and depolarization are excluded. 
Finally, a brief review is given concerning newer experimental 
methods for the observation, recording, and analysis of optical 
scintillation, including suggestions as to what further theo- 
retical and experimental efforts should be undertaken. 


The Joule-Thomson process in cryogenic refrigeration sys- 


tems, J. W. Dean and D. B. Mann, NBS Tech. 
1965), 30 cents. 

A comprehensive analysis of the Joule-Thomson process as 
applied to eryogenic refrigeration systems is presented. The 
descriptions of the process already in the literature are usually 
for specific applications. In contrast, performance charac- 
teristics are presented here for helium, para-hydrogen, and 
nitrogen operations over a large range of process parameters. 


227 (Feb. 14, 


A tabulation of Airy functions, H. T. Dougherty and M. E. 
Johnson, NBS Tech. Note 228 (Sept. 18, 1964), 20 cents. 
A tabulation is provided for Wait’s formulation of the Airy 
function and its first derivative. The argument covers the 
range of real values from —6.0 to +6.0 in intervals of 0.1. 
As an aid for their application the tabulations are also graphed. 


Table of attenuation as a function of vane angle for rotary- 
vane attenuators (A 40 log; cos 6), W. Larson, NBS 
Tech. 229 (Jan. 1965), $1.25. 

The table of attenuation as a function of vane angle gives 
the attenuation in decibels produced by an angular displace- 
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ment of the vane of a rotary-vane attenuator assuming ideal 
behavior of the attenuator. The attenuations for vane 
angles corresponding to the function A 40 logy cos @ 
are given in the table, where A is the attenuation in decibels 
and @ is the vane angle in degrees. The attenuations in 
decibels are given for increments of 0.001° for the values of 6 
from zero to 89.999° to the sixth decimal place and for the 
values of 6 from 89,000° to 89,999° to the fifth decimal place. 
The table was designed for determining the attenuation in 
decibels of a rotary-vane attenuation utilizing a scale marked 
in degrees. Also, the table can be used to determine correct 
decibel markings for a rotary-vane attenuator with a direct 
reading scale regardless of gear ratio. The table can be used 
as an aid in analyzing the calibration data of a rotary-vane 
attenuator for characteristics which include inherent error. 
resettability, and gearing eccentricity. 


Practical aspects of the use of ac-de transfer instruments, 
E. S. Williams, NBS Tech. Note 257 (Mar. 9, 1965) 15 cents. 
Electrothermic transfer instruments may be used to make 
direct measurements of ac-de difference or frequency influence. 
With a somewhat similar procedure an ‘‘unknown”’ a-c volt- 
age or current may be accurately measured by determining 
the difference between it and a preset and accurately measured 
d-c equivalent. Test circuits and procedures are described 
and data taking and calculations are illustrated. 


Summary of current research on archival microfilm, C. 8S. 
MecCamy and C. I. Pope, NBS Tech. Note 261 (Apr. 16, 
1965), 25 cents. 

The discovery of aging blemishes on microfilm in storage has 
prompted a study of the statistical relationship of incidence 
to conditions of preparation and storage, the chemical and 
physical nature of the blemishes, and measures to be taken 
to prevent their formation. This paper is a report of prog- 
ress in that study with recommendations based on current 
information. With few exceptions, reports of aging blemishes 
refer to negative silver-gelatin type microfilm in roll form. 
Extremely small amounts of information loss have been re- 
ported. In some large collections of microfilm, no blemishes 
have been found. No evidence of biological attack has been 
found. Apparently, image silver is ionized, migrates, and is 
reduced to colloidal silver. The reaction often proceeds in 
the manner described by Liesegang for other silver-gelatin 
systems. In laboratory tests, blemishes formed in the pres- 
ence of minute amounts of peroxides, certain imbibed chemi- 
cals, and solid contaminants. Recommendations include: 
avoiding excessive densities on film, careful avoidance of 
physical damage to the image layer, the use of a small con- 
centration of iodide ion in the fixer, thorough washing, uni- 
form drying, storage at low temperature and humidity in 
sealed cans, careful avoidance of air-borne reactants, increased 
use of positive copies for archival storage, and careful periodic 
inspection of record films. 


A calorimetric determination of the enthalpy of graphite from 
1200 to 2600° K, E. D. West and SS. Ishihara (Proc. 3rd 
Symp. Thermal Physical Properties, Purdue Univ., Lafayette, 
Ind., Mar. 1965), Book, Advances in Thermophysical 
Properties at Extreme Temperatures and Pressures, pp. 146-151 
(American Mechanical New York, 
N.Y... 1965). 

Measurements are reported of the increase in the enthalpy 
of high-purity graphite between room temperature and 1200 
to 2600° K. A new apparatus is described briefly. Methods 
of operation are developed which eliminate some systematic 
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and random errors and are especially advantageous at high 
temperatures. The enthalpy data fit the equation 


H7— Hyg 1: 


28.9004 T—1.045 * 10-4 T? 


16,126.2 logy 7'/313.15—8907.3 
where H is in absolute joules per gram atomic weight and 7 


in degrees Kelvin. 


Impedance of sealed nickel-cadmium dry cells, R. 
and H. J. DeWane, Electrochem. Tech. 3, No. 
(Jan.-Feb. 1965). 


J. Brodd 
1-2, 12-17 





The impedance of various sizes of nickel-cadmium dry cells 
was measured as a function of frequency between 50 eps and 
50 kpes. The effects of aging, discharge, and temperature 
on the impedance of the cells were determined. The tempera- 
ture coefficient of electromotive force was determined for some 
cells. The discharge characteristics of AA-size nickel- 
cadmium and Leclanche cells were compared on two standard 
tests. It was found that the impedance of nickel-cadmium 
cells increased during discharge, decreased about 2 to 3 
percent per °C increase in temperature, increased slightly 
during storage, was less for larger cells, and that for cells of 
equal rated capacity, the sintered type had a smaller im- 
pedance than the pressed-powder type. 


Resins and technics used in constructing dentures, (;. C. 
Paffenbarger, J. B. Woelfel, and W. T. Sweeney, Dental 
Clinics of North America, 51-262 (W. B. Saunders Co., 
Philadelphia, Pa., 1965). 

Four technic dentures and usually eight clinical dentures were 
made with eight different types of acrylic resins, a polystyrene, 
a vinyl-acrylic copolymer, an epoxy resin, and hard rubber. 
The dentures were processed with the equipment recom- 
mended by the manufacturer. 

Thick dentures had less shrinkage on processing and less 
dimensional change during use, when wetted, dried in air, and 
heated in water. Therefore, the denture should be made as 
thick as comfort and function will permit. 

None of the dimensional changes of dentures that occurred 
during processing or in use could be detected clinically. The 
conventional acrylic resins processed with the usual dental 
technic of compression molding produced dentures that were 
just as stable in dimension as those produced by special resins 
and applicances. 


The polarographic analysis of titanium (IV) EDTA complex: 
Application to paint pigments, H. W. Berger and B. C. Cadoff, 
Offic. Digest Paint Tech. 37, No. 480, 28-34 (Jan. 1965 

The Ti*4-EDTA complex, buffered at pH 4.7, has been 
found to give polarographic waves suitable for the analysis of 
TiO, in paint pigments. A linear relationship between 
diffusion current and concentration of Ti*4 for the range 2.8 » 
10-5 to 8.4x10-3M has been observed. The pigments 
analyzed have been either in the dry form or extracted from 
whole paint. The method gives results in good agreement 
with the more complicated and time-consuming standard wet 
chemical method. 


A high-speed (milliseconds) method for the simultaneous 
measurement of enthalpy, specific heat, and resistivity of 
electrical conductors at high temperatures, A. Cezairliyan 
(Proc. 3rd Symp. Thermophysical Properties, 1965 Book, 
Advances in Thermophysical Properties at Extreme Tempera- 
tures and pressures, Pp. 253-262 (American Society of Me- 
chanical Engineers, New York, N.Y., 1965). 

A method for the simultaneous measurement of enthalpy, 
specific heat, and electrical resistivity of electrical conductors 
at high temperatures (1000 °K to the melting point of the 
sample) in experiments of short duration (milliseconds) is 
described. The rapid direct joule heating is achieved by a 
single energy pulse obtained from a battery. The average 
heating rate was approximately 6600 °K s-'. The sample is 
tubular in form with a small hole in the middle to approximate 
blackbody conditions, thus allowing the measurement of 
transient true temperatures by photoelectric methods. 
Results on molybdenum (up to approximately 1600 °K) are 
obtained and compared with available literature values. 


Application of cooled air detector in an optical system, 8S. 
Abramowitz, A. M. Bass, and A. E. Ledford, Jr., Appl. 
Opt. 4, No. 2, p. 255 (Feb. 1965). 

A new method for cooling photoconductive detectors is 
described in which the detector is mounted on a copper rod 
attached to a refrigerant. An ellipsoid which is mounted in 
a common vacuum chamber with the detector is 
focus the emergent beam onto the detector. 


used to 


On the relative roles of free volume and activation energy 
in the viscosity of liquids, P. B. Macedo and T. A. 


Litovitz, 
J. Chem. Phys. 42, No. 1, 245 


256 (Jan. 1, 1965). 
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The deficiencies in the rate theory approach of Eyring and 
the free volume theory of Cohen and Turnbull for liquid 
viscosities have been reconsidered. The rate theory has been 
reformulated using the Cohen and Turnbull expression for 
the probability of finding a vacant site and yields an equation 
of the form, 

] Ag exp [Ey RT y Vo vr. 
This result has been applied to many liquids ranging from 
fused silica to liquid argon and including polyatomic van der 
Waals as well as hydrogen bonded liquids. 
Consistent fits to both temperature and pressure dependence 
of viscosity were obtained. A significant conclusion is that 
the free volume of liquids is considerably higher than that 
suggested by the WLF and Cohen and Turnbull equations. 


Zeeman filter, D. G. Dessler and W. G. Schweitzer, Jr., J. Opt. 
Soc. Am. 55, No. 3 284 288 (Mar. 1965). 

A description is given of the Zeeman Filter, a very narrow 
optical band-pass filter. The filtering action is achieved by 
Zeeman splitting a resonant absorption line in a magnetic 
field. With the 2537 A-resonance absorption line of Hg! 
such a filter can produce a filtered line as narrow as .005 em°™!. 


Triggered fuse for load protection, J. H. 
Instr. 35, No. 12, 1715 (Dec. 1964 

This is a circuit designed to add positive protection for 
sensitive loads connected to high capacity power supplies. 
It is especially useful for protecting klystron microwave 
oscillator tubes. It is reliable and uses a minimum of com- 
ponents, hence it is economical and easy to construct. 


Rogers, Rev. Sci. 


Safety in the use of liquid hydrogen, D. B. Chelton, Book, 
Technology and Uses of Liquid Hydrogen, ed. R. B. Scott, W. 
H. Denton, and C. M. Nicholls, ch. 10, pp. 359-378 (Pergamon 
Press, Inc., New York, N.Y., 1964). 

Safety in systems using liquid hydrogen is a subject of 
considerable present-day importance. Reducing or eliminat- 
ing the potential safety hazards involved requires a basic 
understanding of the physical and chemical properties of 
hydrogen and the application of known safety technology. 
These factors are reviewed and recommendations are given 
for safety criteria. Comparisons are made with propane and 
methane to assist in establishing the relative hazard. Con- 
sideration is given in the use of liquid hydrogen on a large 
scale and on a laboratory scale. 


The performance of lenses made from inhomogeneous glasses, 
FW 


W. Rosberry, Appl. Opt. 4, No. 1, 21-24 (Jan. 1965). 
The performance of lenses made from glass elements contain- 
ing defects such as bubbles, streaks, striae, ete., has been 
measured and an attempt made to compare the effect of the 
various defects. Four defective and a good or control lens, 
were used. The performance in the image plane was de- 
termined by four different evaluation procedures. A sum- 
mary evaluation showed that the lens containing stria, al- 
though the most difficult defect to detect, produced the lowest 
quality image. 


Stress-induced Martensite transformation in 18Cr-8 Ni Steel, 
R. P. Reed and C. J. Guntner, Trans. Met. Soc. AIME 230, 
1713-1720 (Dec. 1964). 

A commercial 18Cr—SNi iron alloy (AISI 304L) was examined 
in tension at 300°, 76°, 20°, and 4°K. Continuous 
strain recordings were made, X-ray analyses at periodic stress 
(strain) intervals were obtained, and the magnetic measure- 
ments were taken. From this data the percentage of marten- 
sitic products [bee (a) and hep (e)] was computed as a function 
of stress (strain). It was found that up to 15 pet e phase 
forms at low temperatures. The amount of e formed increases 
to a maximum at about 5 pet strain, then decreases. This 
decrease indicates the additional transformation of ¢« to a’. 
The total amount of « and @’ was suppressed at constant stress 
(strain) at 4°K as compared to 76°K. Itis proposed that the 
suppression of « and a’ is associated with the decreased mo- 
bility of extended dislocations at very low temperatures. 
The yield strength decreased as the temperature was de- 


stress- 


Metallized polymer film guard rings for dielectric measure- 
ments, A. H. Scott, Mater. Res. Std. Tech. Note 5, No. 

76-78 (Feb. 1965). 
Dielectric measurements have been made using a metallized 
polymer film as a guard electrode over a tinfoil guarded 
electrode on three disk specimens. Comparisons with meas- 
urements using the regular tinfoil guard-ring indicate that 
errors of 1.59% or more in capacitance and 3X 10-4 or more in 
dissipation factor can be expected under certain conditions. 


All-glass sorption vacuum trap, F. B. Haller, Rev. Sci. Instr. 
35, No. 10, 1356-1357 (Oct. 1964). 

An ultra high vacuum sorption trap consisting entirely of 
glass is described. The trapping medium is a commercially 
available 96°) silica porous glass. Advantages offered are 
freedom from dust and capability of mounting in any position. 


Comparisons of tearing-strength 
Carter, J. Am. Leather 
1965). 

Slit-tear, double-hole stitch-tear, single-hole stitch-tear, and 
tongue-tear tests were made on different types of leather to 
determine their relative merits. Analysis of the data shows 
that the slit-tear and the tongue-tear tests are more sensitive 
than the double-hole stitch-tear and the single-hole stiteh- 
tear tests. The double-hole stitch-tear test was found to be 
the least sensitive. The tearing processes in the slit-tear 
and tongue tests are similar, and differ considerably from the 
tearing processes in the double-hole and the single-hole 
stitch-tear tesults show the tearing load to be un- 
related to the direction of tear. 


tests for leather, T. J. 
Chem. Assoc. LX, No. 1, 4-14 (Jan. 


tests. 


Thermal insulation storage, transport and transfer of liquid 
hydrogen, hk. B. Jacobs, Book, Technology and Uses of Liquid 
Hydrogen, ed. R. B. Scott, W. H. Denton, and C. M. Nicholls, 
ch. 4, pp. 106-148 (Pergamon Press, Inc., New York, N.Y., 
1964). 

This manuscript is a chapter in a book entitled ‘“‘The Tech- 
nology and Uses of Liquid Hydrogen” to be published by 
Pergamon Press. It presents information useful to the de- 
signers and users of equipment for the storage, transport, and 
transfer of liquid hydrogen. The topics discussed are, in 
general, confined to those which are unique to cryogenic 
equipment; those concepts and techniques which are used in 
connection with equipment that operates at normal tem- 
peratures are not discussed. The chapter is divided into 
three sections: Thermal Insulation, Containers, and Transfer 
Systems. The discussion of thermal insulation includes 
information for non-insulated equipment, for non-vacuum 
insulated equipment, and for vacuum insulated equipment. 
examples of other areas that are discussed include instru- 
mentation, cooldown requirements, pressurization, valves, 
pumps, and fluid mechanical design. 


The response of microwave refractometer cavities to atmos- 
pheric variations, R. O. Gilmer, R. E. MeGavin, and B. R. 
Bean, Proc. 1964 World Conf Radio Me teorology and 11th 
Weather Radar Conf., pp. 244-247 (1964). 

The microwave refractometer has been used extensively in the 
study of the refractive index structure of the atmosphere. 
Little work has been done evaluating the sampling cavity of 
the refractometer as an accurate measuring probe of the 
short-term variations of the atmosphere. Wind tunnel, water 
flow tunnel, and free atmospheric tests of the microwave 
sampling cavities currently in use have yielded information on 
aspect and velocity sensitivity, and on the spacial resolution 
of the cavities. 

The tests revealed that the air flow through the sampling 
cavity undergoes a transition from laminar to turbulent type 
flow at wind velocities in excess of 5 mph. The pressure in 
the cavity decreases as the wind velocity increases. The 
indicated error due to the pressure drop at 15°C is as high as 
1 N unit at 135 mph for a cavity oriented into the wind. This 
error appears to be a function of the aspect of the cavity 
relative to the direction of the wind, becoming a minimum at 
an angle of 45 


pressed below room temperature and then increased rapidly | Flow tunnel tests gave direct visual indications of the flushing 


near 4°K. 


time versus the aspect angle of the sampling cavities. A 
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characteristic flushing length of approximately 
complete flushing was measured at 0 
fect at 90° 


Characteristics of dosimeter films processed in phenidone- 


thiosulfate monobaths, M. 
1, 1-9 (Jan.—Feb. 1965). 
Characteristics of commercial x-ray dosimetry films exposed 
to photons of widely different energies were 
processing in phenidone-thiosulfate monobaths and in x-ray 
and phenidone developers. Monobath processing has a 
number of advantages for routine dosimetry applications. 
There is little dependence of response on processing time and 
temperature and, because of lower contrast, the useful ex- 
posure range is extended. But difficulties due to rate de- 
pendence are more likely, and film sensitivity is reduced. 
One of the main features of the sensitometric behavior and the 
grain structure obtained with monobath processing is the 
interplay of the characteristics of ‘“‘physical development,” 
prevailing at low density levels, and those of ‘‘chemical 
development” at high densities. In contrast to their appear- 
ance for conventional development, the silver aggregates 
corresponding to fog produced by the monobath and the 
aggregates corresponding to relatively low densities above fog 
are large, smooth, and compact. The corresponding photo- 
metric equivalent is high. For higher exposures and densities, 
the differences in photometric equivalent and in the structure 
of the image silver for conventional and monobath processing 
tend to disappear. This points toward a change in the 
predominant development mechanism with density level. 


Ehrlich, J. Phot. Sci. Engr. 9, No. 


The calibration of the National Bureau of Standards’ tritiated- 
toluene standard of radioactivity, S. B. Garfinkel, W. B. 
Mann, R. W. Medlock, and O. Yura, Intern. J. Appl. Radia- 
tion Isotopes 16, 27-33 (1965). 
The difficulties encountered, and their solution, when eali- 
brating a radioactivity standard of tritiated toluene by com- 
parison with the NBS tritiated-water standard by 
liquid-scintillation counting, are described. 
seals of commercial vials, giving rise to uncontrolled oxygen 
quenching, were eliminated by modifying the seals. White 
plastic caps gave higher count rates than black plastic caps. 
A new and highly precise method of dilution was used. 
The ratio of the activities of the two standards was deter- 
mined with a standard error of less than 0.1 percent. 


means of 


Leaks in the 


Conductance of glass immersed in molte 
J. Electrochem. Soc. 112, No. 2, 208-210 (Feb. 1965). 

The electrical conductance of high SiO. glass Was measured 
when immersed in a number of molten salts and salt mixtures 
between 800° and 900°C. Sodium and lithium ions diffuse 
into the glass, both from their own chlorides and from low con- 
centrations in KCl. The glass comes to equilibrium with 
each melt composition, the time of equilibration depending 
on the particular ions in the melt. The activation energy of 
sodium and lithium ion conduction in glass is 75-84 kJ (18 
20 keal), with potassium it is variable between 105 kJ (25 
keal) and 167 kJ (40 keal). The conductance method is not 
sufficiently precise to detect silver and cuprous ion diffusion 
into glass. 


n salts, K. H. 


Stern, 


Other NBS Publications 


J. Research NBS 69A 

1965), 70 cents. 

Electronic arreaeure and magnetic Beene rties of the neptunyl 
ion. J. Kisenstein and M. H. Pryce. 

Heats of eee in Pioadae oxide by differential 
thermal analysis. E. M. Levin and C. L. MeDaniel. 

Phase relations between iridium and the sesquioxides in air. 
S. J. Schneider, J. L. Waring, and R. E. Tressler. 

Phase equilibrium relationships in the system Gd:0;-TiOs. 
J. L. Waring and 8. J. Schneider. 

Acid-base behavior in 50-percent aqueous methanol: 
dynamies of the dissociation of protonated tris (hydroxy- 
methyl) aminomethane and nature of the solvent effect. 
M. Woodhead, M, Paabo, R. A. Robinson, and R. G. 
Bates. 

Observations of dislocations and surface features in corundum 


(Phys. and Chem.), No. 3 (May-June 


thermo- 


2.5 feet for 
as compared to 5.75 


compared for 


crystals by electron transmission microscopy. ). J. Barber 
and N. J. Tighe. (See above abstracts.) 

Low-frequency dielectric properties of liquid boric 
K. H. Stern. 

Mass spectrometric study of the effects of nitric oxide and 
other additives on “‘purified’’ active nitrogen. J. T. 
Herron. 

Synthesis and ring structure of 
galacto-heptulose. E. J. McDonald. 


oxide. 


7-acetamido-7-deoxy-.L- 


J. Research NBS 69B (Math. and Math. Phys. 

(Jan._June 1965), 75 cents. 

Lectures on matroids. W. T. Tutte. 

Menger’s theorem for matroids. W. T. Tutte. 

Single-element extensions of matroids. H. H. Crapo. 

Minimum partition of a matroid into independent 
J. Edmonds. 

Lehman’s switching game and a theorem of 
Williams. Edmonds. 

On the connection between the properties of oriented linear 
graphs and anlyses of lumped pysical systems. H. M. 
Trent. 

Character subgroups of 
Newman. 

A note on multipliers of difference sets. 

Modification of Edmonds’ maximum 
C. Witzgall and C. T. Zahn, Jr. 

On measurable sets and functions. A. J 

Characteristic formulation for nonperiodic solutions of 
Mathieu’s equation. T. Tamir and H. C. Wang. 

On the surface duality of linear graphs. J. Edmonds. 

Maximum matching and a polyhedron wlth 0, 1-vertices. 
J. Edmonds. 


Radio Sci. J. Res. NBS/ USNC-URSI 69D, No. 6 (June 1965). 

$1.00. 

Cyclotron harmonic 
ford. 

Propagation of waves across a magnetoplasma-vacuum bound- 
ary. R. L. Gallawa. 

Phase velocities and attenuation distances in the ionosphere. 
D. R. Croley, Jr., and B. 8. Tanenbaum. 
abstracts). 

lonospheric effects of electrostatic fields generated in the outer 
magnetosphere. G. C. Reid. 

An experimental study of plasma sheath effects on antennas. 

+. Tyras, P. C. Bargeliotes, J. M. Hamm, and R. R. Schell. 

An approach to improve re-entry communications by suitable 
orientations of antenna and static magnetic field. S. N. 
Samaddar. 

Concerning the 


, No. 1 and 2 


subsets. 


Tutte and Nash- 


F-groups. M. I. Knopp and M. 
R. A. Brauldi. 
matching algorithm. 


. Goldman. 


waves in warm plasmas. F. W. Craw- 


(See above 


mechanism of reflection of electromagnetic 
waves from an inhomogeneous lossy plasma. J. R. Wait. 


Selection of Papers Presented at 1964 World Conference on 
tadio Meteorology, Boulder, Colo., Sept. 14-18, 1964 


Angels in focus. JD. Atlas. 
An investigation of clear air 


stratification 
elevated instruments. D. R 


. Hay and K. 

On inferring the refractive-index structure 
from electromagnetic scattering 
Smith, Jr. 


with radar 
Naito. 

of the troposphere 
experiments. P. cs 


and 


The biexponential nature of troposphe + , gaseous absorption 
of radio waves. E. J. Dutton and B. Bean. 

Complete scattering parameters of ee haeanet hydrome- 
ter in the A 0.1 to 410 ecm range. D. Deirmendijian. 


Radio Sci. J. Res. NBS/USNC-URSI 69D, No. 7 (July 1965) 

$1.00. 

International een. agg of atomic frequency standards via 
VLF radio signals. H. Morgan, E. L. and B. E. 
Blair. (See above Ft acts). 

Control of WWV and WWVH standard freque a y broadcasts 
by VLF and LF signals. B. E. Blair and A. H. 


(See above abstracts). 


Crow, 


Morgan. 


Measurements of the total electron content and the equivale nt 


slab rg ae of the 


midlatitude 
Bhonsle, 


ionosphere R. 
and O. K. 


. da Rosa, Garriott. 
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D-region absorption at 10 and 15 Me/s during the total solar 


eclipse of July 20, 1963. G. M 
and J. M. Watts. 

Effect of the eclipse of 20 July 
gating over short paths. J. H 
A comparison of radar auroral r 

wave theory R. L 
Electromagnetic 
D. J 
Influence of an inhomogeneous ground on the propagation of 
VLF radio waves in the t waveguide. 
J. R. Wait. 
Aspects of the terrestrial ELF 
source or its antipode. L. G 
Analvsis of linear arrays focused in the Fresnel 
P. P. Lombardini, R. Doviak, and J. Goldhirsh. 
Theorv of coil antennas. T. Padhi. (See above abstracts). 
Small prolate spheroidal antenna in a dissipative medium. 
R. H. Williams, R. D. Kelly, and W. T. Cowan. 
Calculated curve for groundwave propagation over inhomo- 
geneous earth 1 topographical features 


K. Furutsu 


Lerfald, J. K. Hargreaves, 
1963 on VLF signals propa- 
Crary and D.E. Schneible. 
flection data with acoustic 
Leadabrand 
properties of a plasma 
Jacavanco 


covered antenna. 


earth-ionosphere 


noise spectrum when near the 
Abraham, Jr. 


region. 


pronoul Ce 


Radio Sci. J. Res. NBS/USNC-URSI 69D, No. 8 (Aug. 1955 
$1.00 


“Ultra Low Frequency Electromagnetic Fields’”’ Papers 


Preface to “Ultra 
Papers. W. H. Campbell and 38. 

Solar wind and its 
C. P. Sonett. 


Schumann 


Eleectromagnetie Fields’’ 
Matsushita. 


with the 


Low Freque I CV 


interactio! magnetosphere. 


resonances. J. Galejs. 
Karth-ionsophere cavity res 
ELF radio waves. J. R 


Resonances of the 


mances and the propagation of 
Wait 
earth-ionosphere cavity 
Kngland M. J. Rycroft. 
Experimental results on the dynamics of the F 
Becker, R. Riister, and J. Klostermeyer. 
Regular oscillations near 1 c/s observed atrmiddle and low 
latitudes. L. Tepley. 
Preliminary results of a micropulsation experiment at 
conjugate points. R. E. Gendrin and V. A. Troitskaya. 
Some characteristics of geomagnetic pulsations at frequencies 
W. H. Campbell and E. C. Stiltner. 
Propagation of hy lromagnetic waves in the magnetosphere. 
M. Sugiura 


Ionospherie perturbation 


observed at 
Cambridge, 


region. W. 


near 1 ¢/s 


the roles played by the ionosphere 
tic pulsations S. Matsushita. 
Effects of induced earth currents on low-frequency electro- 
magnetic oscillations. A. T. Price. 
Equatorial effects. R. Hutton. 


In geomagne 


Interpretation of early magnetic 


transients caused by high- 


altitude uclear detonations. S. L. Kahalas and P. 


Newman. 


Abstracts of ULF Conference papers not 


issuc 


published in this 


A note on the application of pulst 


compression techniques to 
ionospheric sounding dD. C 


Coll and J. R. Storey. 

Comments on a paper “Measurement of the phase velocity of 
VLF propagation in the earth ionosphere waveguide” 
by F. K. Steele and C. J. Chilton. H. F. Bates. 

Reply to H. F. Bates’ ‘Comments on ‘Measurements of the 
phase velocity of VLF propagation in the earth ionosphere 
waveguide’ by F. K. Steele and C. J. Chilton.” F. Kk. 
Steele and C. J. Chilton 


Report of the 49th National Conference on 
Measures 1964, NBS Mise. Publ. 263 (Feb. 1, 

Technical highlights of the National 
Annual 1964, NBS Mise. 
SL.00. 

Thermal properties of aqueous uni-univalent 
V. B. Parker, NSRDS-NS 2 (Apr. 1, 1965), 45 cents. 

Pipe nipples: Brass, copper, steel, and wrought-iron, CS5—65 
{ Jan. 15. 1965 10 cents. Supe rsedes CS5 16. 


Weights and 
1965) $1.00. 
Bureau of Standards 

Publ. 264 (Dec. 1964 


teport 


elect roly tes, 


Builders’ template hinges, CS9-65 (Jan. 15, 
(Supersedes CS9—-33). 

Aluminum alloy chain link fencing, CS269—65 (Mar. 1, 
10 cents. 

Ponderosa pine windows, sash, and screens (using single glass 
and insulating glass), CS163-64 (Mar. 17, 1964) 
Supersedes CS163—-59 and CS193-53. 

Wood double-hung window units, CS190—64, (Mar. 17. 
10 cents. Supersedes CS190-59. 

Wood awning window units, CS204—64 
10 cents. Supersedes CS204—59 

Wood casement window units, CS205-64 
10 cents. Supersedes CS205-59. 

Wood horizontal-sliding window units 
CS264-—64 (Mar. 17, 1964), 10 cents 

Wood horizontal-sliding window units 
operating sash), CS265-64 (Mar. 17, 

Wood single-hung window units, 
10-cents. 

Quarterly radio noise data September, October, November 
1963, W. Q. Crichlow, R. T. Disney, and M. A. Jenkins 
NBS Tech. Note 18-20 (Oct. 23, 1964), 50 cents. 

Quarterly radio noise data December, January, February, 
1963-64, W. Q. Crichlow, R. T. Disney, and M. A. Jenkins, 
NBS Tech. Note 18-21 (Jan. 25, 1965), 50 cents. 

Bibliography of fading on microwave line of sight tropospheric 
propagation paths and associated subjects, 1. T. Dougherty, 
NBS Tech. Note 302 (Aug. 31, 1964), 60 cents. 

Bibliography on tropospheric propagation of radio waves, W. 
Nupen, NBS Tech. Note 304 Apr. 1, 1965), $2.00. 

Atlas of Fourier coefficients of diurnal variation of foF». 
Part II. Distribution of amplitude and phase, W. B. 
Jones and R. M. Gallet, NBS Tech. Note 305 (Feb. 14, 
1965), 60 cents 

Resonance effects of electrostatic oscillations in the 
osphere, J. A. Fejer and W. Calvert, J. Geophys. 
5049-5062 (Dee. 1964 

Observatory report on JILA, L. M. Branscomb, Astron. 
J. 69, 700-703 Nov. 1964 Nonlinear ambipolar dif- 
fusion of an isothermal plasma across a magnetic field, 
Kk. R. Mosburg and K. B. Persson, Phys. Fluids 7, 1829 
1833 (Nov. 1964 

Laboratory studies of helium ion loss processes of interest 
in the ionosphere, KE. ; Ferguson, F. C Fehsenfeld, 
D. B. Dunkin, A. L. Schmeltekopf, and H. I. Schiff, 
Planetary Space Sei. 12, 1169-1171 (Dee. 1964 

Interaction of maxima in the absorption of soft X-ray, J. W. 
Cooper, Phys. Rev. Letters 18, No. 25, 762-764 (Dee. Zi, 
1964 

Complex Lorentz group with a real metric. I. 
ture, A. O. Barut, J. 
1964 

Climatology of elevated super-refractive layers arising from 
atmospheric subsidence, B. A. Cahoon and L. P. 
Proc. 1964 World Conf. Radio Meteorology and 
Weather Radar Conf., pp. 52-55 (1964). 

Annual report from the National Bureau of Standards, C. E. 
Moore, Astron, J. 68, No. 9, 1-3 Nov. 1964 

Nuclear magnetic moment of Pr'*! from the hyperfine struc- 
ture of doubly ionized praseodymium, J. Reader and J. 
Sugar, Phys. Rev. 137, No. 4B, B784—B789 (Feb. 22, 1965). 

Brightness temperature of the atmosphere using a biexpo- 
nential model in the 6-45 GlIz/s frequency range, EK. J. 
Dutton, B. R. Bean, and Ek. R. Westwater, Proc. 1964 
World Conf. Radio Meteorology and llth Weather Radar 
Conf., pp. 212-215 (1964 

Natural variations in the abundance ratio and the atomic 
weight of copper, W. R. Shields, 8. 8. Goldich, KE. L. 
Garner, and T. J. Murphy, J. Geophys. Res. 70, No. 2 
179-491 (Jan. 15, 1965). 

Optical heterodyne detection of the forward stimulated 
Brillouin seattering, D. A. Jennings and H. Takuma, Appl. 
Phys. Letters 5, 241-242 (Dec. 1964 

The uses of oblique ionograms in frequency utilization, K. 
Davies, Telecommun. J. 31, 281-287 (Oct. 1964 

Conductance of glass immersed in molten salts, K. H. Stern, 
J. Electrochem. Soc. 112, No. 2, 208-210 (Feb. 1965). 

Reversal of the diurnal phase variations of GBR (16 kilo- 
cycles per second) observed over a path of 720 kilometers, 


1965), 10 cents. 


1965), 


15 cents. 
1964), 
March 17, 1964), 
(Mar. 17, 1964), 
(all sash operating), 
one or more 


1964), 10 cents. 
CS266-64 (Mar. 17, 


non- 


1964), 


ion- 


Res. 69, 


Group struc 
Math. Phys. 5, 1652-1656 (Nov. 


tiggs 


Lith 
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. D. Crombie and H. L. 
5023-5027 (Dee. 1, 1964). 
properties of ilmenite-hematite 
C. F. Jefferson and R. G. West, 
No. 11, 544-548 (Nov. 1964). 

Thermodynamics of hydrogen solubility in cryogenic solvents 
at high pressures, M. Orentlicher and J. M. Prausnitz, 
Book, Chem. Eng. Sei. 19, 775-782 (1964). 

Liquid air, R. B. Scott, Encyclopaedic Dictionary of Physies, 
p. 1 (Pergamon Press, Inc., New York, N.Y., 1964). 

Properties of normal and para hydrogen, R. B. Stewart and 
H. M. Roder, Book, Technology and Uses of Liquid 
Hydrogen, ed. R. B. Seott, W. H. Denton, and C. M. 
Nicholls, ch. 11, pp. 379-404 (Pergamon Press, Ine., New 
York, N.Y., 1964). 

High resolution, low energy electron spectrometer, } oo 
Simpson, Rev. Sci. Instr. 35, No. 12, 1698-1704 (Dee. 1964). 

The calculation of clectron energy distribution functions in 
the ionosphere, L. R. Megill and J. H. Cahn, J. Geophys. 
tes. 69, No. 23, 5041-5048 (Dee. 1, 1964). 

Solar flares with associated active dark flaments and their 
relation to 2800 Me/s radio bursts, M. W. Haurwitz, 
Astrophys. J. 140, No. 3, 1236-1246 (Oct. 1, 1964). 

Mechanism of inactivation of bacteriophages by metals, N. 
Yamamoto, C. W. Hiatt, and W. Haller, Biochim, Biophys. 
Acta 91, 257-261 (1964). 

Electron-impact studies of aromatic 
Azulene and naphthalene, R. J. Van 
Wacks, J. Chem. Phys. 41, No. 10, 
1964). 

Karly strength, flow and dimensional changes obtained on 
amalgam perpared with a standardized mechanical technic, 
H. J. Caul, W. 8S. Crowell, W. D. Kimmel, and G. C. 
Paffenbarger, J. Am. Dental Assoc. 69, 742-748 (Dee. 
1964). 

Photochemical changes in thin layer chromatograms of poly- 
evelic, aromatic hydrocarbons, M. N. Inseoe, Anal. Chem. 
36, No. 13, 2505-2506 (Dec. 1964). 

Surface integral form for three-body 
mann equation, M. 8. Green, Phys. Rev. 
A905-A910 (Nov. 16, 1964). 

Nucleation characteristies of static liquid nitrogen and liquid 
hydrogen, J. Hord, R. B. Jacobs, C. C. Robinson, and L. L. 
Sparks, J. Eng. Power 86, 485-495 (Oct. 1964). 

Degree of metamerism, I. Nimeroff and J. A. Yurow, J. 
Soc. Am. 55, No. 2, 185-190 (Feb. 1965). 

Food for billions, A. T. MePherson, 
forum “Food in the Future: Concepts for 
Chicago, Ill., Oet. 8-9, 1964, pp. 
Industries Supply Assoc., 1964). 

Is life possible on other planets? Ww R. 
and Children 2, No. 4, 7 (Dee. 1964). 

Investigation of plasma boundaries with 
surface waves, H. W. Wassink and A. J. 
Phys. 35, No. 10, 2795-2800 (Oct. 1964). 

Comments on a paper by R. L. Closs, “‘An experimental 
investigation of back seattering or radio waves from the 
equatorial electrojet,”” R. Cohen and D. T. Farley, Jr., 
Proc. Phys. Soc. 84, Pt. 4, No. 540, 619-620 (J. W. Arrow- 
smith Ltd., London, England 1964). 

Changes in living, ‘g FE. Eberhard, Proce. 
vironment of Change, Sterling Forest, 
June 14-17, 1964, pp. 62-67 (1964). 

Liquefied handling of; R. B. Scott, Eneylopaedic 
Dictionary of Physics, p. 1 (Pergamon Press, Ine., New 
York, N.Y. 1964). 

Low-temperature apparatus, constructional materials for, 
R. B. Seott, Eneyclopaedic Dictionary of Physics, pp. 1-2 
(Pergamon Press, Inec., New York, N.Y. 1964). 

Electron-impact studies of aromatic hydrocarbons. IT Naph- 
thacene, chrysene, tripenylene, and pyrene, M. FE. Wacks, 
J. Chem. Phys. 41, No. 6, 1661-1666 (Sept. 15, 1964). 

Formation of NF from NF, by photolysis, D. E. Mann and 
J. J. Comford, Spectrochim. Acta 21, 197-198 (Pergamon 
Press Ltd. Northern Ireland 1965). 

The influence of metastable oxygen molecules on ozone and 
airglow, H. I. Schiff and L. R. Megill, J. Geophys. Res. 
69, No. 23, 5120-5121 (Dee. 1, 1964). 


tath, J. Geophys. Res. 69, No. 


solid solutions, 
J. Am. Ceram. Soc. 47, 


hydrocarbons. IIT. 
Brunt and M. E. 
2395-3199 (Nov. 15, 


Boltz- 
No. 4A, 


in the 
136, 


collision 


Opt. 


tomorrow's Proc. 


Planning,” 
51-738 (Dairy and Food 


Bozman, Science 


electromagnetic 


I’stin, J. Appl. 


Conf. The En- 
Tuxedo, N.Y., 


ZAses, 


Shape of the magnetospheric boundary under 
pressure, R. J. Slutz and J. R. Winkelman, J. 
tes. 69, No. 23, 4933-4948 (Dee. 1, 1964 

Standards—the foundation of industrial growth, | 
Stiehler, Industrial Times VI, No. 22 
Nov. 15, 1964). 

Hyperfine structure of the B?=>* state of CN, H. E. Radford, 
Phys. Rev. 136, No. 6A, A1571—A1575 (Dee. 14, 1964). 

Research for architecture, J. P. Eberhard, Res. News II, 
No. 1, 10-19 (Jan. 1965). 

Spectral properties of plants, D. M. 
J. C. Schleter, and V. 

11-20 (Jan. 1965). 

Blood:bone equilibrium — in 
MacGregor and W. E. 
359-361 (Jan. 23, 1965). 

Carbon monoxide on niobium, R. Klein and J. W 
Surface Science 2, 167-176 (North-Holland Publ. 
Amsterdam, The Netherlands, 1964). 

Fracture in corrosive media, H. L. Logan, Book, Fracture of 
ngineering Materials, ch. 8, pp. 143-161 American 
Society for Metals, Metals Park, Ohio, Nov. 1964). 

Gas-phase radiolysis of n-pentane. A study of the decom- 
positions of the parent ion and neutral excited pentane 
molecule, P. Ausloos and 8. G. Lias, J. Chem. Phys. 41, 
No. 12, 3962-3970 (Dec. 15, 1964). 

Moderation in instrumentation, J. K. Taylor, Science Teacher 
32, No. 3, 18-19 (Mar. 1965). 

Pressure-indicator-paste patterns in duplicate dentures made 
by different processing technies for the 
J. B. Woelfel and G. C. 
Assoc. 70, 339-343 (Feb. 

Radioisotope dilution and its application to the radioassay of 
Cel, L. A. Currie, G. M. Franee, III, and H. L. St : 
Intern. J. 
1965). 

teference buffer solutions for pH measurements in 50% 
methanol. Dissociation constants of acetic acid and di- 
hydrogen phosphate ion from 10 to 40°, M. Paabo, R. A. 
Robinson, and R. G. Bates, J. Am. Chem. Soc. 87, 415-418 
(1965). 

Spectra of matrix-isolated NiF, and NiCL, D. E. 
M. E. Jacox, and J. D. MeKinley, J. Chem. Phys 
3, 902-905 (Feb. 1965). 

The structures of anhydrous sodium trimetaphosphate, 
Na,P;QOy, and the monohydrate, NasP,0,-H,O, H. M. Ondik, 
Acta. Cryst. 18, pt 2, 226-232 (Feb. 1965). 

Vacuum — ultraviolet photochemistry. VIII. Propylene, 
D. A. Becker, H. Ikabe, and J. R. MeNesby, J. Phys. 
Chem. 69, No. 2, 538-542 (Feb. 1965). 

A magnetic amplifier for use with diode logic, E. W. 
IEEE Trans. Commun. Electron. 83, No. 75. 
(Nov. 1964). 

A study of the variation of the surface electrical resistance of 
lead iodide films with relative humidity at room tempera- 
ture, F. E. Jones (Proce. 1963 Intern. Symp. Humidity and 
Moisture, Washington, D.C.), Book, Humidity and Mois- 
ture. Measurement and Control in Science and Industry. 
1, Principles and Methods of Measuring Humidity in 

Sec. II], Paper 37, pp. 361-365 (Reinhold Publ. 
Corp., New York, N.Y., 1965). 

An investigation of psychrometric measurement techniques 
in air-conditioning calorimetry, J. C. Davis and P. R. 
Achenbach (Proc. 1963 Intern. Symp. Humidity and Mois- 
ture, Washington, D.C.), Book, Humidity and Moisture. 
Measurement and Control in Seience and Industry. 2 
Applications, Sec. IV, Paper 37, pp. 303-313 
Publ. Corp., New York, N.Y., 1965 

Calculation of the temperature of a flat-plate wet surface 
under adiabatic conditions with respect to the Lewis 
relation, T. Kusuda (Proce. 1963 Intern. Symp. Humidity 
and Moisture, Washington, D.C.), Book, Humidity and 
Moisture. Measurement and Control in > and 
Industry. 1, Principles and Methods of Measuring 
Humidity in Gases, See. I, Paper 2, pp. 16-32 (Reinhold 
Publ. Corp., New York, N.Y., 1965). 

Clinical behavior of o-ethoxybenzoie acid-eugenol-zine oxide 
cements, S. Civjan and G. M. Brauer, J. Dental Res. 44, 
No. 1, 80-83 (Jan.—Feb. 1965). 
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Gates, H. J. Keegan, 
R. Weidner, Appl. Opt. 4, No. 1, 


calcium homoeostasis, J. 
grown, Nature 205, No. 4969, 


Little, 
Co., 


same 
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1965). 


path hits, 
Dental 


inberg, 


Appl. Radiation Isotopes 16, No. 1, 1-8 (Jan. 
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Control of relative humidity and temperature in rubber 
laboratory of National Bureau of Standards, F. L. Roth 
and R. D. Stiehler (Proc. 1963 Intern. Symp. Humidity 
and Moisture. Washington, D.C.), Book, Humidity and 
Moisture. Measurement and Control in Science and 
Industry. 2, Applications, Sec. III, Paper 35, pp. 287-292 
Reinhold Publ. Corp., New York, N.Y., 1965). 

Measurement of radiation dose distributions with photo- 
chromic materials, W. E. McLaughlin and L. Chalkley, 
Radiology 84, No. 1, 124-125 (Jan. 1965). 

Soil corrosion, W. J. Schwerdtfeger, Encyclopedia of Electro- 
chemistry, pp. 1068-1071 (Reinhold Publ. Corp., New 
York, N.Y., Mar. 1965). 

A theorem on the automorphs of a skew-symmetric matrix, 
M. Newman, Mich. Math. J. 12, 61-63 (1965). 

On the location of the zeros of some infrapolynomials with 
prescribed coefficients, J. L. Walsh and O. Shisha, Pacific 
J. Math 14, No. 3, 1103-1109 (1964). 

On the structure of infrapolynomials with prescribed co- 
efficients, O. Shisha, Pacific J. Math 14, No. 3, 1039-1051 
(1964). 

Continuous-absorption hygrometry with a pneumatic bridge 
utilizing critical flow, W. A. Wildhack, T. A. Perls, C. W. 
Kissinger, and J. W. Hayes (Proce. 1963 Intern. Symp. 
Humidity and Moisture, Washington, D.C.), Book, Hu- 
miditvy and Moisture. Measurement and Control in 
Science and Industry. 1, Principles amd Methods of 
Measuring Humidity in Gases, Sec. VI, Paper 53, pp. 
552-570 (Reinhold Publ. Corp., New York, N.Y., 1965). 

Absolute measurement of W for polonium-210 alpha particles 
in nitrogen, argon, and argon-methane mixture, Z. Bay 
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